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Über ein ebenes unstetiges Temperaturspannungsproblem 
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Zur Auffindung der Wärmespannungen in einem ungleichmäßig erwärmten Körper selzte man bisher die 
Differenzierbarkeit der Wärmeverteilungsfunktion voraus. Diese Vorausselzung wird in der vorliegenden Arbeit 
fallen gelassen, und als Beispiel werden die Wärmespannungen in einem unendlich langen Zylinder bei Vor- 
gabe einer unstetigen Temperaturverieilung ermittell. 


To determine the thermal stresses in an unequally heated body, the heat distribution was, till now, sup- 
posed to have a derivative. In this paper, this hypothesis is abandoned, and as example, the thermal stresses 
in a cylinder of infinite length are determined, if a discontinuous heat distribution is assumed. 


Jusqu’a present on supposait que, pour deter miner les tensions ihermales dans un corps inegalement 

. chauffe, on powvait differentier la fonction de la distribution de la chaleur. Cette hypothese est'abandonnee 

dans cet article, ei pour evemple les iensions thermales dans un cylindre de longueur infinie sont deter- 
minees ü la supposition d’une distribution inconstante de la chaleur. ; 


UToöbt OnpeneluTb TePMHYecKue HAIPAIKCHHA B HEPABHOMEPHO HaTpeBaeMoM Tee NO CHX 
nop Nperuonaranach AuddepeHumpyeMOcTb PyHEIHH TEIIOBOTO pacnpezenenns. B HacToameh 
paboTe 3Ta IIpenuochLIka OTOPAachIBaeTcH, & B KAyecTBe MpMMepa BEIABIAWOTCH TEePMHYeCKHE 
HAIPpS:KeHHA B ÖECKOHEeYHO AIIHHHOM IHJIHHNPE IIPH 3a1&HHOM HePaBHOMEePHOM PACHpenelIeHHH 
TeMHEePaTyPB. 


Die in dieser Arbeit verwendeten Bezeichnungen sind: 
E der Elastizitätsmodul. 


Em 
een h E 
art) der Schubmodul 


m das reziproke Poissonsche Querkontraktionsverhältnis. 

&j die Verformungskomponenten einer nicht Kirchhoffschen Formänderung, der 
kein Verschiebungsvektor entspricht. Wahlweise wird auch gesetzt: &]] — &, Egg = & 
Egg = &z, E13 — Eaı — Yay» Eas — Esa — Yya » Esı — E13 — Var» 

0; die Spannungen des im Gefolge der obigen Verzerrungen im Körper vorhandenen 
Spannungszustandes (oder wie es heißt: Selbst- oder Eigenspannungszustandes). 
Wahlweise wird auch gesetzt: 0,1 =, =8, (og =,y=8, 0,3 =, =8= 071 = (g 
= Toy» Oag = Toy — 033 > Oz, = (13 = Top 

w, die Komponenten eines Verschiebungsvektors, der vom unverzerrten Anfangszustand 
aus gerechnet wird. 

e;; der vom undeformierten Anfangszustand angefangen und aus dem Verschiebungsvektor 
abgeleitete Verzerrungszustand. (Kirchhoffscher oder verträglicher Verzerrungs- 
zustand). Wahlweise auch: e,ı = & , &a2 = & , &3 = & » E12 = xy USW. 

6; =0füri #7, ml für =]. 


Problemstellung 
Die bisherigen Rechnungen, die Wärmespannungen in ungleichmäßig erwärmten Kör- 
pern zu finden, fanden statt unter der Voraussetzung der Differenzierbarkeit der Wärme- 
verteilungsfunktion. Der auf F. Neumann [1] zurückgehende Ansatz von W. Voigt/[2] 


[| m % \ , mEa oT 
Im ad Ta man, 
worin T die Differenz zwischen Körper- und Außentemperatur und & die Wärmeausdehnungs- 
zahl sind und die übrigen Größen die übliche Bedeutung haben, setzt differenzierbare Tempe- 
raturverteilung voraus. Die gleiche Voraussetzung gilt auch bei späteren Betrachtungen. 
Unter Beschränkung auf die Gleichgewichtsaufgabe, also auf den Fall ö = 0 sollen nun die 
Wärmespannungen in einem unendlich langen Zylinder bei Vorgabe einer unstetigen Tempe- 
raturverteilung ermittelt werden. Bei Inanspruchnahme der Verträglichkeitsbedingungen 
für den ebenen Verzerrungsfall wird eine inhomogene Bipotentialgleichung zum Ausgang 
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der Untersuchung. Für die Randwerte gilt dabei die Voraussetzung der Spannungsfreiheit 
des Randes. 
Zunächst ist dann eine Poissonsche Differentialgleichung zu betrachten. Die 
Kenntnis eines partikulären Integrales dieser inhomogenen Differentialgleichung verschafft 
die Randwerte zur Bestimmung der Lösung der homogenen biharmonischen Differential- 
gleichung. Das Problem ist dann nach Bekanntwerden der Spannungsfunktion, aus dieser die 
Spannungen auszurechnen. Auf die Ausrechnung der Verschiebung wurde verzichtet. Die 
Spannungen selbst wurden nur für die innere Umgebung des Randes ermittelt. u 


Einleitung 

Zu bestimmen sind die Eigenspannungen im Inneren eines von äußeren eingeprägten 
Spannungen freien, unendlich langen, kreiszylindrischen Körpers, der ungleichmäßig er- 
wärmt ist. 

Die Temperaturverteilung soll sich zunächst noch mit den Koordinaten des Kreisquer- 
schnittes stetig ändern. Von der Koordinate längs der Achsenrichtung des Zylinders soll die 
Temperatur nicht abhängen. 

Über diesen Gegenstand sind schon zahlreiche Untersuchungen gemacht worden, 
darunter auch solche mit zeitabhängiger Temperaturverteilung. In allen Fällen beschäftigte 
män sich mit der Herstellung der Ausdrücke für die Spannungen und die Verschiebungen. Die 
Temperaturverteilung wurde natürlich als stetige Funktion der Koordinaten vorausgesetzt. 

Außer dem vollen Kreiszylinder wurde auch der hohle Kreiszylinder betrachtet. 

Franz Neumann[3] untersuchte einen zusammengesetzten Körper. Der Quer- 
schnitt des Körpers ist rechteckig. Die Frage wird bei ihm als gewöhnliches Spannungs- 
problem behandelt. Bei spannungsfreier Oberfläche und dem Fehlen von Massenkräften 
sollte nach dem Eindeutigkeitssatz von Kirchhoff£f[4] auch das Innere eines Kontinuums 
spannungsfrei sein. Beim Beweis dieser Behauptung ist vorausgesetzt, daß die Verzerrungs- 
komponenten sich eindeutig aus einem stetigen Verschiebungsvektor ableiten lassen. Weiter 
wird vorausgesetzt, daß zwischen dem Verzerrungstensor und dem Spannungstensor eine 
eindeutige Zuordnung besteht. Zudem werden die Randbedingungen jeweils am unverformten 
Körper befriedigt. 

Die Verletzung einer dieser Voraussetzungen macht die Existenz mehrerer Lösungen 
des Elastizitätsproblems möglich. Geben wir die Voraussetzung der eindeutigen Bezogenheit 
von Verzerrungs- und Spannungstensor zueinander auf, so ist ein Selbstspannungszustand 
möglich. Dieser kann bedingt sein durch eine Unstetigkeit des Verschiebungsvektors oder 
dadurch, daß die Verzerrungen nicht eindeutig aus einem Verschiebungsvektor abgeleitet 
werden können. 

Die Unstetigkeit im Verschiebungsvektor kann hervorgebracht werden durch eine 
starre Verschiebung zweier Ebenen gegeneinander, wobei nach Vornahme der Versetzung diese 
Flächen wieder fest miteinander verbunden worden sind. Die Kräfte und Momente der 
solcherart erzeugbaren Distorsionen lassen sich aus sechs Elementardistorsionen herleiten, 
indem man die resultierenden Dynamen an den Schnittflächen in Komponenten zerlegt. Die 
erzeugten Spannungszutände sind im allgemeinen räumlich. 

Sind die Eigenspannungen nicht durch Unstetigkeiten der Verschiebungen, sondern 
durch Verzerrungen ohne zugehörige äußere Spannungen, aus denen ein Spannungstensor 
bestimmt werden könnte, aufgebracht, so läßt sich die Aufgabe durch Definition eines ver- 
träglichen Verzerrungstensors lösen. | 

Das Auftreten nicht axialsymmetrischer Eigenspannungsverteilungen macht zudem noch 
eine Annahme über die Eigenschaften der Selbstspannungsursachen in Abhängigkeit von der 
Koordinatenrichtung der Zylinderachse notwendig. 

Die Unabhängigkeit der Eigenspannungsursachen von dieser Koordinate rechtfertigt 
nicht die Annahme einer gleichen Tatsache für die Eigenspannungen selbst. 

Die Gewährleistung des ebenen Verzerrungszustandes ist bei endlicher Zylinderlänge 
nur gegeben bei der zusätzlichen Forderung, daß eine bleibende Dilatation in Richtung der 
Zylinderachse nicht erfolgt. Allen solchen nur von den Zylinderkoordinaten abhängigen Ver- 
formungen ist gemeinsam, daß man keine mechanische Ursache dafür kennt und daß dem- 
zufolge eine eindeutige Zuordnung zwischen dem Verschiebungsvektor und den Verzerrungs- 
größen nicht stattfindet. Diese nicht von eingeprägten Kräften oder Spannungen verursachten 
Verformungen können durch elektrische, magnetische, thermische oder andere physikalische 
und chemische Ursachen hervorgebracht werden. Solche temporären oder auch bleibenden 
Verzerrungen werden nicht ausschließlich Dilatationen wie im Falle einer reinen Dichte- 


. sein müssen, sondern sie können auch gestaltändernden Einflüssen ihre Entstehung 
verdanken. 
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Nach notwendigen Festsetzungen über die Werte der Elastizitätskonstanten beim Auf- 
treten innerer inkompatibiler Verformungen, die ihre Nachbarschaft in einen Spannungs- 
zustand versetzen, kann man das Hookesche Gesetz in folgender Form ansetzen 


Se | Ta) | s ) | 
e; = Dr er — 2m, == 30. Gt 2G8; —) ENT 27 - 
Berücksichtigt man Massenkräfte nicht, so lauten die Gleichgewichtsbedingungen 
el Ber 
0%; 


Die zugehörigen Oberflächenbedingungen ausgedrückt durch die Eigenspannungen lauten 
0,508 (N, %;)=0. 


Bereehnung der Eigenspannungen aus bekannten Anfangsverformungen 


Die Tatsache, daß Eigenspannungen einem Verzerrungstensor nicht eindeutig zugeord- 
net werden können, hat häufig zur Folge, daß man sie direkt aus den Gleichungen des 
Hookeschen Gesetzes berechnet, wenn die Gesetze der Anfangsverformungen bekannt sind 
oder wenn sinnvolle Annahmen über diese gemacht werden können. 

Kehrt man das Hookesche Gesetz mit Verwendung der Abkürzungen 9 und % für 
die Summen e;; und &;; um, so ergibt sich 


2 : 
= @ ee 9) Hey eyl du t=1,2,3. 
Setzt man dies in die Gleichungen der Gleichgewichtsbedingungen re 03,9 == 1 2,5 
ein, so erhält man 0%; 
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Das heißt, daß ein Eigenspannungszustand immer entstanden gedacht werden kann durch 
nichtkirchhoffsche Verformungen und elastische Verzerrungen in einem ursprünglich unver- 
formten Körper. 

Ordnet man nun den auseinem Verschiebungsvektor mit den Komponenten u; ableit- 
baren Verzerrungen e;;einen Spannungstensor über die Beziehungen 


0%; 2 


zu, so erhält man die Relation 
00;; | 2 0» 085; =) 5 
Er ae -. er 
2 ep rue i 


nicht über i summieren ! 
Der kompatible Verzerrungstensor e;; ist einem Spannungstensor zugeordnet, der in 
einem anfänglich spannungslos gedachten Körper durch die Volumkräfte 
2 0% 08; = . 
ee =D 
de — 20%, 2 0% 0%, I . 
aufgebracht zu denken ist. 


Be. 
Bei reiner Dilatation würde sich daraus wegen &; —=(0, &; =. (nicht summieren) 
ergeben 
0 3—2+m Em 9% 
de ur d— = 
en 3(m — 2) 3 (m — 2) 0%; 


ist also eine Potentialfunktion der Volumkräfte. 
Daraus ist zu ersehen, daß die Volumkräfte im Falle von reinen Dichteänderungen als 


Gradienten des Potentials ® ausgedrückt werden können und daß an der 


Em 
3 (m — 2) 
Berandung des Körpers eine fiktive Spannung vom Betrage dieses Potentials anzubringen ist. 
Nach der Ausrechnung der Spannungskomponenten aus den Divergenzbedingungen der Eigen- 


; Em 

spannungen muß die Größe 3m 3) ö davon abgezogen werden, d.h. aber nur von den 
m — 

a ( P 2 00,5 : d Ei w 

Normalspannungen. Die verschwindenden Divergenzen werden ja nur von den Eigen 


0%; 
9* 


ge 


Elastizit all 
bereitstellen kann, zu errechnen. Die 
Einführung einer Hilfsgröße in Anspruch ge 
Verzerrungstensor der Deformation aus dem unv 
“auch mit Verwendung der für diesen Tensor gültigen atibilit: 
nach den Eigenspannungen als definierte mathematische Aufg en; 
d.h. daraus eine Randwertaufgabe der Elastizitätstheorie zu machen. 


' S i he u ” 
Berechnung der Spannungen aus den Kompatibilitätsbezieh 
Die zu Anfang in einem Körper vorhandenen vorübergehenden oder 


nichtkirchhoffschen Verzerrungen genügen wie gesagt den Kompatibilitätsbe« ing u 


Fe 
Ihnen ist kein Verschiebungsvektor so zugeordnet, daß seine Komponenten sich unabhängig 
vom Integrationsweg durch Linienintegrale 


ou; 2 es 
u = (w)+ | de; = 2 


darstellen lassen. 

ou; 
0% Ei, | 
Integrationsweg sein. Sie müssen also gleichfalls vollständige Differentiale sein, d.h. es müssen 
die Gleichungen gelten: 


Gelten diese Gleichungen für die w,, so müssen die gleichfalls unabhängig vom 


ou; a: 2) 0? U; 2 ö 
O0, \0zyo O0 : 
Damit dies zutrifft, dürfen die folgenden Bedingungsgleichungen nicht verletzt werden: 
o2e, Pe, ey N o2e, Es ae Oeyz I — En 
oa oa Roy 0:02 0% oe ..0y 03 PT 


Das sind die Kompatibilitätsbedingungen. 
Nun genügen aber weder die Anfangsverformungen noch die Spannungen diesen Glei- 
chungen. Die Kirchhoffschen Verzerrungskomponenten aus dem unverformten An- 
fangszustand und die ihnen zugeordneten gedachten Spannungen genügen ihnen. 


Damit hat man aber, wie wir weiter unten an einem Beispiel zeigen, im Falle des zwei- 2. 
dimensionalen Elastizitätsproblems eine definierte Randwertaufgabe vor sich, sofern die 
Anfangsverformungen bekannt sind. == Bi 

Man hat, um das einzusehen, nur die Gleichungen 'r“ 


e; = Re Gi + 2 Ge; > (0 + =) 6 L 


in die Kompatibilitätsbedingungen einzusetzen und erhält, da die e,;,; diesen genügen, u.a. 

die folgenden Gleichungen: 
e1: I Bye DM (5 | 6 ta A De we 

2G\loy 0 m+1\oa: ' 92 020) OyE | dm may’ \ 
1 ( 20, 1477025 o ( Or OT, 2) O?e, f) Oyy, Or MR 
or m+10y% 9a\ a 02 2, 0% + Ze + m) 


Man erhält also ein System von sechs unabbhängigen Differentialgleichungen, welche die 
gesuchten Spannungen in Beziehung setzen zu den bekannten Anfangsdeformatiönen oder 
auch umgekehrt, bei bekannten Eigenspannungen Relationen bilden mit den gesuchten verur- 
sachenden Deformationen. In dieser Allgemeinheit wird das Problem der. Auffindung der 
Spannungen aus den Anfangsbedingungen, von Sonderfällen abgesehen, wenig Lösungsaus- 
sichten haben. Man wird sich daher, um ein Urteil über den physikalischen Sachverhalt zu 
gewinnen, auf durchsichtige ebene Probleme beschränken. 


Das zweidimensionale Selbstspannungsproblem 


. Wir betrachten nunmehr zur konkreten Durchführung einer Frage aus dem Gebiete der 
Eigenspannungen einen Körper von kreiszylindrischem Querschnitt, dessen Länge im Ver- 
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gleich zu seinem Durchmesser groß ist. Die Zylinderachse machen wir zur z-Achse des einzu- 
führenden Koordinatensystems. 

Die Temperaturverteilung, der wir uns den Körper ausgesetzt denken, sei nur von den 
Polarkoordinaten r und @ des Querschnittes abhängig. 

Wenn der Körper ungleichmäßig erwärmt wird, ist er in seiner Beer Ausdehnungs- 
möglichkeit behindert. Es werden auf diese Weise in ihm Spannungen auftreten, ohne daß 
äußere Spannungen dabei mitwirken. 

Nach den Voraussetzungen über die Temperaturverteilung hängt diese von z nicht ab. 
Damit sind auch die durch die ungleichförmige Temperaturverteilung wachgerufenen Span- 
nungen mit Ausnahme derjenigen in der Umgebung der Stirnfläche des Zylinders unabhängig 
von ihrer Lage in der z-Richtung. 

Das wird bedeuten, daß die Zylinderquerschnitte in sich eben bleiben müssen. Ver- 
formungen dieses Charakters nennt man ebene Verzerrungszustände. Aus Symmetriegründen 
verschwinden hierbei die Schubspannungen 7,, und r,,. Außerdem wollen wir die Massen- 
kräfte, da sie den Selbstspannungen einfach überlagert werden können, nicht ansetzen. 

Damit schreiben sich die Gleichgewichtsbedingungen mit ya =0, Ya —=0 und 
e, —consk. 

00; 


in. Se ae 92 
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wenn hinzukommend beachtet wird, daß keine der Größen von z abhängt. 
Bei Spannungsfreiheit der Oberfläche lauten die zugehörigen Randbedingungen: 


0; c08 (n, 2,;) = 0 71,2% 
Die sechs Verzerrungskomponenten sind voraussetzungsgemäß ebenfalls unabhängig von 
z. Da zudem 7,, = (0, 7, = O0 und e, = const., bleibt nur die eine Kompatibilitätsbedingung 
o?e, ö8e, Ya 
OUT 0 ON 


übrig. Nun gilt 


ee] ( 8 
Baal, 201-0, 4) 
oder 
0, = E(e, —&T) — — (0, + 0,) 
und dann 
1 } 8 
= Ru + (26aT Oj > 


denn eine reine Dichteänderung, wie es die Wärmeausdehnung ist, verursacht wenigstens 
primär keine Winkeländerungen der Volumelemente. 
Mit dem Ausdruck 


5, =E(,—aT) — — (0, 0,) 


können die Verzerrungen auch folgendermaßen geschrieben werden or, 3): 
4 s—0,\, | e m-+1 
Ber 2 - A (6 | m ) di J r lmau n) du. 


Setzt man diese drei Verzerrungsgrößen in die Kompatibilitätsgleichung ein, so erhält man die 
folgende partielle Differentialgleichung: 


0°C 027, 0°0, Ema _, 
a een Nee 
Nun befriedigt nach Air y der folgende Ansatz: 
OF o®F o®F 
A May ET day’ Buy] 


die Gleichgewichtsdivergenzen. Setzt man diese Ausdrücke in die Differentialgleichung der 
Spannungen ein, so erhält man die folgende inhomogene Bipotentialgleichung für F 
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‘kennt, da die Integrale der homogenen Gleichung sämtlich in der homogenen Bipotential- 


m—1 r 


mit der Bedingung verschwindender Spannungen auf dem Rande. 
Setzt man eine stetige Verteilung der Temperatur voraus, so lasser 


tionen W* = — ®und V* aus den beiden Differentialgleichungen | 
NN i . BISu En”: 


bestimmen. Dabei errechnet man zuerst aus der Poissonschen Gleichung W*, m 
Randbedingungen zur Bestimmung von V zu erhalten. Es genügt dabei, wenn man m 
partikuläres Integral der Gleichung : 


p2W* 2 uRann 
m — 1 


gleichung enthalten sind. Die Spannungen selbst stellen sich inden V und W wegen F< RE 
als der Spannungsfunktion in folgender Weise dar: | — 
02V. 02 W. oa 02lW. 02V 2W 


a a ag ran 
Hat man ein partikuläres Integral der inhomogenen Poissonschen Differentialgleichung 
gefunden, so ist die Funktion V durch Bestimmung der Lösung von VY?V = 0 mit den Rand- 
werten 


nl 


dev _doW dv dow Sm Be 


dd ey aeg a; 

zu bestimmen. = Y 

Setzt man voraus, daß die Ableitungen auf dem Rande eindeutig sind und sich stetig m 

ändern, so läßt sich mit Hilfe der stetigen Temperaturverteilung 7 eine solche Funktion z 
finden [5]. ag 

Dabei wurden bisher vorwiegend unendlich lange zylindrische Körper mit Kreisquer- N 


schnitt betrachtet. Für Sonderfälle stetiger Temperaturverteilungen wurde die Frage auch bei 
Hohlzylindern beantwortet. Bei einfach berandeten Bereichen, die durch Polynome die kon- 
forme Abbildung des Einheitskreises auf sich gestatten, kann das Eigenspannungsproblem 
gleichfalls für den langen Zylinder gelöst werden. 

Ferner gelingt es dann, die Frage nach den Eigenspannungen in zusammengesetzten 
Körpern aufzugreifen und für den Kreisquerschnitt bei konzentrischer Anordnung zu beant- 
worten [6]. 

Gibt man die Stetigkeit der Temperaturverteilungsfunktion auf, indem man sie stück- 
weise stetig oder stückweise konstant annimmt, so hat die so gebildete Funktion keinen defi- 
nierten Gradienten an den Sprungstellen mehr. 

Die Lösung der Bipotentialgleichung ist dann aus den Werten zu suchen, die das partiku- 
läre Integral der inhomogenen Differentialgleichung elliptische Typs und seine erste Ableitung k 
auf dem Rande annehmen. Diese erste Ableitung ist nun aber keine differenzierbare Funktion. 

Die zu erwartenden Lösungen dieses Randwertproblems werden demnach Diskontinui- 
täten in den zweiten Ableitungen haben. Die Frage ist nun, ob es solche Lösungen des Pro- 
blems gibt, deren Unstetigkeiten mit der Differentialgleichung vereinbar sind, beziehungs- 
weise ob durch geeignete Festsetzungen über die Air ysche Funktion verträgliche Ergeb- 
nisse liefernde Lösung gewonnen werden können. Nachdem man einheitliche konstante Größen 
für das Elastizitätsverhalten und die Wärmeausdehnung festgesetzt hat, soll also auch die für 
stetige Dichteänderungen als Folge ungleichmäßiger Temperaturverteilung aufgestellte Diffe- 
rentialgleichung bestehen bleiben. 

Die aus der Lösung dieser Differentialgleichung durch zweimalige partielle Differentiation 
zu gewinnenden Spannungen müssen so beschaffen sein, daß die Normalspannungen in Rich- 
tung der Normalen der Unstetigkeitsfläche und die Schubspannungen die Diskontinuitäts- 
fläche stetig durchsetzen. Eine Unstetigkeit dürfen nur die in der Tangentialebene dieser Fläche 
liegenden Normalspannungen beim Durchsetzen derselben erfahren, 


ET WER, EBEN OR WERE 


partielle zweite Differentialquotienten überhaupt nicht definiert sind, die Differenzierbarkeit 5: 
auf der Sprungfläche i in der Art aufhört, daß endliche Sprünge auftreten. 


Man wird sie e vielmehr ee müssen und sie dann der inhomogenen elliptischen Diffe- % 
| rentialgleichung i in der Form , 


je J ns 
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als Belegungsfunktion beigegeben. 
Die Lösungen des simultanen Gleichungssystems mit einer diskontinuierlichen Störungs- 
funktion für die Poisson sche Differentialgleichung ist in geschlossener Form nach einer auf 


 Gauchy zurückgehenden Form möglich [7]. 


Es wird sich danach ergeben, daß mit Ausnahme einiger singulärer Randlinien, an denen 


Das ebene Eigenspannungsproblem bei unstetiger Diehteänderung re 


Mit den oben gemachten Voraussetzungen und Einschränkungen gehen wir also von der u 
Spannungsfunktion E 


aus. 

Im konkreten Fall sei nun die Dichteänderung als die folgende unstetige Funktion vor- 
gelegt. Innerhalb eines konzentrischen Kreissektors betrage die Temperaturerhöhung T°, 
während sie im übrigen Teile des Kreises den Wert Null haben möge. In dieser Form ist das 
Problem der Bestimmung der zugehörigen Spannungsfunktion der Berechnung der Durch- 
biegung der allseitig eingespannten Kreisplatte mit unstetiger Lastverteilung sehr ähnlich. 

Harry Schmidt hat für den Fall unstetiger Lastverteilungen auf der längs ihres 
Umfangs eingespannten Kreisplatte ein Integrationsverfahren angewandt, das für die ver- 
schiedensten Arten unstetiger Lastverteilung seine außerordentliche Kraft erwiesen hat. Es 
ist von ihm selbst und seinen Mitarbeitern [7] wiederholt mit Erfolg angewendet worden. 

Um dieses Integrationsverfahren auf das obige Eigenspannungsproblem anwenden zu 
können, transformieren wir die Poissonsche Differentialgleichung und auch die homogene 
Gleichung vierter Ordnung sachgemäß erst auf Polarkoordinaten und erhalten, wenn wir fest- 
setzen 

Ders, - Syn tt denen, 
die folgenden partiellen Differentialgleichungen: 
ca lea Wo 7 
ar re rar 
a 10 1% 
0° 
or? ar r or 2 09% 
Dazu kommen nach Ermittlung eines partikulären Integrales der inhomogenen Differential- 
gleichung für W und V die Randbedingungen 
oVK(r, a - (9) 
r=r N— en 


Mr, 9) =Wlr, y) und as Er 
Die Temperaturverteilung hängt nur ab von o. Sie ist folgendermaßen definiert: 
| 0- für —n So<— m. 
Ir, 9) = T(p) = 1 für -QRSpPp<+tM 
VB für n<PO< N: 
Dafür kann mit Hilfe der Dirichletschen diskontinuierlichen Funktion der Ansatz ge- 


macht werden: H 1 | - ee+ ele—p.) 2 ' | 
(9) 9ari | 2 u | 2 | 5 


(©) (€) 


Adi sen de 


r 5 . an, H Ye h. 
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Das Integral ist dabei längs einer Parallelen zur imaginären Achse in positiver Richtung zu er- 
strecken. Ändert sich pim Intervall —x, +, so springt das erste Integral der Klammer an der 
Stelle 9 = — 9, vom Werte Null auf T°, um von da an konstant zu bleiben. Das zweite Inte- 
gral macht den gleichen Sprung an der Stelle 9 = 99, SO daß von = g,bis = rn die Funktion 
T(g) wieder verschwindet. 

Wegen der Gültigkeit des Superpositionsprinzips genügt es zunächst, ein partikuläres 
Integral der Poissonschen Differentialgleichung 


AW(r, g) = er? 
zu kennen. 
Ein solches Integral läßt sich angeben, es lautet 


Wr: a 

- 77 (+25), 2) 

Zur Herstellung der Lösung des Problems ist hierzu diejenige Lösung der homogenen Gleichung 
vierter Ordnung AAV = 0 zu überlagern, die den Bedingungen 


er? 


UK, o, )r=n = {W (rn, ®, 2)}r=n = na) @ a) — 9(P, ©) 


und 
a en Fe 
er Dr ea 
genügt. 
Dieses Randwertproblem hat für den Kreis die Lösung [8] 
+2 
= 7 —hda 
2 pP — pe al | er 
2 ZITENg = 2 r?+r3 — 2rr,cos(& — p) 


_.n 


Feten | 
Bi {r? +3 — 2Arr,cos (a — 9)}? $ 


Der Ausdruck kann für unsere Zwecke umgeformt werden und lautet 


+n 
nn | di 


Are, Jr +19 — Arrgcos (& — ) 


+7 
= (15 — r?)? | g da 
VERA T 


4arz 5 — !rr,cos (a — ) 


+n 


j g da 
{+ —2rr9cos (a — p)}? 


an (r2 I, r2)° 


An? 


Setzt man schließlich 
EmaT 
irn € u 
MET 


so ist durch 
1 | " eP02 " e-Pı2 
Ar p) er: Iril | x 2 irle P5 2) er Mr, P 2)} da — Ay, {er pP; z) E Wr, pP, 2)} de) 
© a” 
die Lösung der inhomogenen Bipotentialgleichung gegeben. Die Quadratur muß zuerst über 
den Anteil des Integranden erfolgen, der vom partikulären Integral der inhomogenen partiellen 
elliptischen Differentialgleichung herrührt. Denn dieser Teil bestimmt. die Randwerte für die 


Lösung der homogenen Bipotentialgleichung im Falle des Kreises nach der oben stehenden 
Poissonschen Formel. 


Die angekündigte Integration liefert 


y2 eP #902 dz | y? 
ae al Test al) ze 


© lo für oSTo. 
Bei der Bestimmung eines partikulären Integrals der ellipti j j i 
partikuläre ptischen inhomogenen Diff tial- 
gleichung gibt es keine Schwierigkeiten. Anders ist es bei der homogenen Biotonne 
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‘Die Randwerte des zu ihrer Bestimmung dienenden Integrals erfüllen nicht die Bedingungen 
viermaliger Differenzierbarkeit auf dem Rande. Diese ist aber Voraussetzung für die Existenz 
einer regulären Lösung der Bipotentialgleichung. Die Funktion gist nur einmal und nicht vier- 
mal dort stetig differenzierbar und daher gibt es keine reguläre Lösung für das Randwertpro- 
blem. Damit gibt es keine Möglichkeit, die Einzigkeit der Lösung nachzuweisen. 

Weder Af noch nn sind bis in den Rand hinein stetig. Nun zeigt man bei regulären 


Lösungen die Eindeutigkeit auf folgende Weise. Man geht von der Green schen Formel aus, 
indem man in ihr die Funktion «durch Aw ersetzt und hat dann 


; 08V ni 
[au Av — v Adu) do a (1u® — Da ds 
für den Kreis und den Rand. 


Setzt man die gleichfalls reguläre Differenz U zweier regulärer Lösungen dieses Problems 
in die obige Formel ein, nachdem man noch vu = v = U gemacht hat, so ergibt sich mit a 
Randbedingungen 


oU 
sofort die Gleichung 
J(AU% do =. 


Wegen der Regularität von AU folgt daraus AU =. 


Da auf dem Rande U = 0 besteht, muß die Funktion U identisch null sein. Nun ist unsere 
Lösung nicht regulär, also läßt sich der Eindeutigkeitsbeweis nicht führen. Die gewonnene 
Lösung ist also nicht die einzig mögliche. Nach dieser Zwischenbemerkung möge nun die Qua- 
dratur fortgeführt werden. 


Die anderen Quadraturen werden so vorgenommen, daß man die Reihenfolge der auf- 


‘einander folgenden Integrationen vertauscht. Nach Ausführung der Integration nach z sind 
dann noch die folgenden Integrale zu berechnen: 


[ cos {n(a + 9,)} da a [  cos{n(a — Pu) da 
{r? + 13 — 2rr, cos (a — p)YF {r? +r} — 2rr, cos (a — p)}* 
k=1,2,n=0,2. 


Die unteren Integrationsgrenzen sind dadurch festgelegt, daß die Integranden jeweils für 
— nSas— g,0der -asSasoy, verschwinden. 


Die Berechnung dieser Integrale ist sowohl im Komplexen als auch im reellen Gebiet in 
geschlossener Form möglich. 


Für das Reelle sei die Auswertung an dem Integral 


7 
| da 
Jr +r) — 2rr,cos{& — p) 
Po 
vorgeführt. 
Zu dem Zweck setzen wir 


WI): Ant er er 
cos (a — p) = a sin (« — 9) = Are 
a ae RL 
cotg — = u; da = FAyE 


Damit erhält man für das unbestimmte Integral 


da E 
ee ern BEE: = ); = m 


ro 
mereig Ge „colg > — 5 


— 
ro 
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Das bestimmte Integral erhält somit den Wert: 


2 er, ER rt u, |. 
ei 


Die Größe & ist hierin folgendermaßen definiert 
cc 0 für —n<o<— 

—r fir -_ H9<P<I. 

Neben dem reellen Integrationsverfahren wurde abwechselnd auch die Integration im Kom- 


plexen herangezogen, wofür gleichfalls noch ein Beispiel gegeben werde. Zu dem Zweck nehmen 
wir das Integral 


S 


N—P 
Br E dy 
= a TE 2 
{? +13 — Arncos (a — p)}? Artnr a v) 
277, 
Be —(9+9,) 
Setzen wir 
$ r Fo 
z= eiW, kZı =» 20 — —, 
Yo 2 


so kann dieses Integral als Realteil von 
eıl2—P) 
% 2. dz 


KU REIN, (2 — 21) (2 — 23)} ? 


aufgefaßt werden. 


Durch Partialbruchzerlegung erhält man daraus 


e&lT—p) 
ee | a =), Ar I Zal2ı — 2) ie 2ı + = dz 
5 en na 
eo (P— 9.) (2 — 21)? 2 — 2%] (2 — 22)” 2 — 231(2ı — 22)? 
lm) 

FI 2ılıı —2) Zıt2 | 22l&ı = 2) 21+ 22 

f 2 2) 3 2 I 2, —| dz. 
we, en 2 — 2,) PR (2 — 2;) a 


Daraus folgen nacheinander 


el P) 
dz 1 letze) ' ; 
—i(P+9) (2 — 2)? LK — 21\e ie —Pı) emp) — 2, eit7—P) = 2, 
€ 0 | 
_ 008 (P — 90) — For + irssin (P+ Po) __ Na c0os p+ Tor + irasin 
r? +17; — 2rr,cos (p-+9,) r?+n+2rr,c059 
damit 
h eilr—p) 
N ia 2) de | 
y2 a! y2 3 @-2)| 
| 0) Bar 1) 
— 1 j rr,sin (9 + 9) Ps © 


(rt — rl +7? —2 170,008 (P+9) tr +2rr,cosg)' 
Das Integral 


E44 


da 
| {r? + 73 — 2rr, cos (a — po)? 


== 


» 
EN N N “u 
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läßt sich mit diesen Ergebnissen in der folgenden Form schreiben 


27 | rr,sin(p-+r) rr,sin @ 
(r? — rn) Ir? Hr — 2 rr, cos ae -) 
2(r+ 9%) rosin (p+ 2) 
re t ee 5 
de N (ae (+) — {r — r,c08 re, 
Be: arcig r,sin (P + 9) 
Vet — 2rrcos (pt) — {r— 10cos (p + go)} 
Ben on 
(rn) Vr? +2 — 2rrgcos (pP +2) — {rs — r cos ne 


Es 9%) 73 
Yr +13 — 2 rroc0s (P+ 90) — {ro — r cos (P + 0)} 
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Die Arcustangensfunktionen lassen sich zusammenfassen nach den für diese geltenden Addi- 


tionstheoremen. Danach ergibt sich 


arctg rsin (@+ 9%) <<) 
Yr®+ 13 — 2rrcos (P+ 90) — {n— rcos (p-+ u} 
a arctg | i rg sın (9 Ss 2) ) 
Yr® + rg —2rrgcos (P+ 9) — {r— 1, cos (P + 9o)} 
== — arctg | 7 Kur r ° ctg es T 2) - — arctg (e- — ctg Be & e) 5 
0 
Das obige Integral erhält damit den Wert 
2 | rr,sin(pP-+9) rr,sin@ \ 
(r?— 2)? \r2 + r2 a cos (9 + 9) u +72 +2 rr,c0so| 
2(r® + 5) [rn — Y a tor P+nl 
(r? — »2)3 are In, - Er r Se Ban BE Ing - +r „cotg 2 ! 
2 (r® + ro) 
Tea 


Dabei ist & 
er für — Q<p<R 
| Is für — app Sn: 


Weiter erhält man auf ähnliche Weise 


ln 9) 
„inne de | 
a 28 ER 
| ( 0) eilP+ 90) (e 2.) | 
1 f rr,sin (9 + 9) rr,sin® | 


= Mr Hr —2rncos(Pp+9) Frrnt2rn cos oJ” 


dann erhalten wir noch die beiden Integrale 


Ulm Oo) 
„lines + FH 22) | 
er SE A | 
BEIELD RE (sin (p +) 
ern) Yr? +12 — Arr, cos (pP +) — {r—rgcos (P +r)} 
er arte 7, sin (P + 9) 
Vr + — 2rr,cos (P+ 9) — {r — 70008 (P + Po} 


Bor 
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und 
Fa?) | 
rer ae 
er) rs 
— 2(n + r?) H rsin.(o+ 7) 
= are | ——— —an ee E 
=) Yr® +73 — 2rr,cos (p+n) — {n —rcos(p +7 
erregen a = . 
Yr? +3 — 2rrocos (P + 90) — {rn — r cos (P + Po)} 


Der Integrationsweg ist längs des Einheitskreises in der z-Ebene zu erstrecken. Um die aus- 
. . . .. . . > j Yo .. T 
tretenden Funktionen in dieser Ebene eindeutig zu erklären, ist diese von 2, = „ über A 


f) 
längs der negativ reellen Achse aufgeschlitzt. Der Wert der Größen cund chängt davon ab, wie 
oft der Integrationsweg den Schnitt überschreitet. Der Wert der Arcustangensfunktion ist 

2 7 47 4 
durch den veränderlichen Winkel bestimmt und ändert sich zwischen — 67 und + >. Die 


Werte der Größen eund c seien in Abhängigkeit von pin der folgenden Tabelle gegeben 


M n für — rn <o<— 
en A 
Er a > BE — Tue 0 <p<r— 9: 
7 MÄÖP<T 0 für m<P<n. 
Schlitz von = nach — Schlitz von n nach — x. 
0 


Alle auftretenden Integrale können nun auf die eine oder andere Art ausgerechnet werden, so 
daß man die gesuchte Spannungsfunktion in der folgenden Form erhält; mit den Erklärungen 
für den Arcustangens: 


En) Tfg), 


dabei ist 
SE 9) =f im Sektor Iu. II 
2 
=f—- {1- c00s2(@+ 9) im Sektor IIIu. IV 


4 
2 (r2 — r?)? — Art 
A 


= im Sektor Vu.VI 
r r* 
et — 005 2(P +9) +7 (10052 (pP — 9) 
DD. ehe 
Ä 1 ur: | 2 (r ER 4r 


er Bild. 1. 16 7 r 
Einteilung des Integrationsgebietes 9 (r7 u Fa)° em,! 
47? -| Sr 
| 


r& 
: — 0052 (P — gu). 


“cos 2 (+ oo) 2 


16 


Darin ist f eine Funktion, die im Anhang dargestellt ist. Die Art der Einteilung des Kreis- 
querschnittes lehnt sich an die hergebrachte [7] an und geht aus dem Bild eindeutig hervor. 
Die Spannungen ergeben sich aus der Spannungsfunktion durch Differenzieren. Man erhält: 


en 02F a 


m — = —— 
"2 9r2 r or’ Or’ Or 


r 09 


Der durch J dargestellte Teilder Lösung wird im Punkte r — 0 sinnlos, und bei seiner Herleitung 
wurde diese Stelle dann ja auch ausgeschlossen. Wegen der fortgesetzten Division durch r 
wurde stillschweigend r £ 0 vorausgesetzt. Der Wert für fan der Stelle r = 0 kann entweder 


- 
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durch eine Grenzwertbetrachtung, ein Verfahren, dessen man sich im Hinblick auf die späteren 
Differentiationen nicht bedienen kann, oder aber auf dem Wege von Reihenentwicklungen der 
Integranten ermittelt werden, deren Konvergenz in der nahen Umgebung von r = 0 gutist. 
Die betreffende Untersuchung wurde nicht gemacht, da außerhalb dieser Stelle der ge- 
schlossene Ausdruck zur Ermittlung der Spannungen zur Verfügung steht. Zunächst wurde 
durch entsprechende partielle Differentiation der Spannungsfunktion die Schubspannung be- 
stimmt. Dafür erhält man wegen i 


EmaT 0) IE s 

9% aan; und u 7% 

S 1 £ 3 no Hleor 
unter Verwendung des im Anhang dargestellten Differentialquotienten — andre den fol- 
genden Ausdruck er 

BBELEE a 1 in den Sektoren Iu. II 
m—1il 9r\r pl 

a Eu. Zr a in den Sektoren Illu. IV 
m—ıl 9ar\r 
> 13n —2rr+rt 

Snap in 2lpr 2) 

2 4 r 

= 2 BET al 14 - ST in den Sektoren Vu. VI 
m—1ı|l ar\r 


4a} u: 
+57 5in2(P+ 9) — z sin 2 (9 — 9) 


13n —2rnr-+r,. uch 
en: Samt (sin 2(p+ 9) —sin 29-9)! 
Setzt man darin r = r,, so verschwinden, wie verlangt, die Schubspannungen. In den Punkten 
r=ny,P=—-9,undr=r,P=+ ist die Funktion nicht definiert. An den Grenzflächen 
9 = — 9, und 9 =+ 9, verhält sich Schubspannung mit Ausnahmen stetig. Die Ausnahme- 


stellen liegen für die hergestellten Funktionen auf den Erzeugenden des Zylinders durch die 
Randpunkte r,, —9, und 7, +99. Die Verhältnisse auf der Achse r = 0 können damit nicht 
beurteilt werden. Es würde die Untersuchung auf diese Punkte auch nicht weiter ausgedehnt. 


Der für die Schubspannung r gefundene Ausdruck wurde auf sein Verhalten in den Stellen 
fg —P, und r,, +9, untersucht. Dazu wurden in den Funktionen 5 + r? — 2 rır cos (9 — 9,) 
und 3 -+r? — 2 r,rcos (p + 9,) abhängigen Gliedern diese beiden Argumente in der Um- 
gebung der kritischen Punkte durch I ersetzt. Die anderen jeweils interessierenden Größen 
lassen sich dann durch / und einen Hilfswinkel y, der durch die Beziehung: = rn +1? — 
2r,lcos y definiert ist, darstellen. Es zeigt sich dann daraus, daß r auf den fraglichen Linien 
kein Wert entspricht, da die Grenzübergänge jeweils Werte erzeugen, die von der Richtung 


2 
abhängen. Mit den genannten Ausnahmen läßt sich nun auch die Normalspannung o, = = 
definieren. Die Werte für 02f/or?ersehe man im Anhang. T 

2 
Op = Era. im I. und II. Abschnitt 
$  _ m— 19 
EmaTy0°f 


De im III.und IV. Abschnitt 
», m—1 lor? 


Ar, 1 yd — Zr ; 
Er 2 (1 — cos 2 (9 ai P0)) = (8 — re: (COS 2 (o -- 0) 
ige ee im V. und VI. Abschnitt 
HU — Y 


= b3 i 
ES (1 — cos2(p+9)) + hd cos 2 (p — 9,)) 


1 Ar 1 n— 3 
zE ze ca n — 0052 (@- ) 1 (2 = ” c0s 2 (op — 0) i 


134 Hieke, Über ein ebenes unstetiges Temperaturspannungsproblem  pa.34 Nr 4 Apri/Ma1 1054 
ee ss 2 1 o%f 

Die Normalspannung 0, läßt sich mit den im Anhang aufgeführten Werten für 73 dg% und 
m in folgender Weise darstellen: 


Ge et je = a an im I. und II. Abschnitt 
"T m—1 Ir 99% ' rorl 
2 nn May im III. und IV. Abschnitt 
ml Wer 
1 
— 0052 (pP +90) —z (1 — 005 2(P-+ 9) 
8 (rn — r?)? — 16 r* 1 4(r, + r*) 
RN n cos2(P+9o)+ Ak: — am cos2 (pP -+9)) 
0, >= en En im V. und vr. Abschnitt 


m—1\r 99% 7 or 

1 1 \ 
— 008 2(p-+ 90) + cos 2 (pP) — z (1 — 005 2(P+ 90)) + Z(1-e0s2(P- Po)) 
8 (rg — r2)? — 16%, 


16 r# = 
4 4 
nr (cos 2 (P + 96) — ©05 2 (pP — P))) 


Diese Funktion ist im Nullpunkt nicht erklärt. Desgleichen nimmt sie auf den Erzeugenden 
des umschließenden Zylindermantels durch die Punkte r,, —9@, und 79, +9, keinen Wert an. 
Damit ist für die innere Umgebung der Zylinderfläche mit den genannten Ausnahmen für diese 
selbst die Aufgabe, die Eigenspannungen für die geschilderte unstetige Dichteverteilung auf- 
zufinden, gelöst. Einschränkend wurde noch verlangt, daß die Dichteänderung bei endlicher 
Zylinderlänge nur in Richtungen senkrecht zur Achse des Zylinders stattfindet. Ein Einzig- 
keitsbeweis läßt sich nicht führen, da die Funktion nicht mehr als einmal differenzierbar ist 
Parallel zur Achse, also in z-Richtung, wirkt die Spannung 


2 
,=b&,—aT)+— (+). 


cos 2(P + 9) — c0s 2(p — 9o)) 


Zusammenfassung 


Die Herstellung der Ausdrücke für die Spannungen durch gehörige Differentiation der 
Lösung der Randwertaufgabe der Bipotentialgleichung, also der Air yschen Spannungs- 
funktion ergibt einen stetigen Durchgang der Normalspannung in Richtung der Flächennor- 
male der Diskontinuitätsfläche durch diese Trennfläche. Die Schubspannung durchsetzt die 
Trennfläche gleichfalls ohne Sprung. Allein 
die Normalspannung parallel zur Tangential- 
ebene zur Trennfläche erfährt eine Diskon- 
tinuität. Die Unterbrechung der Stetigkeit 
besteht in einem endlichen Sprung. 

Dabei wurde das Verhalten der Spannun- 
gen im Nullpunkt nicht untersucht, da für 
diesen die Air ysche Funktion nicht gesucht 
wurde. 

Eine zweite Einschränkung ist zu machen 
bei den Punkten der zwei Mantellinien des 
Zylinders durch die Punkte r,— 9, und 
"g +99. Auf diesen Linien existieren für die 
Ableitungen keine Werte. Die Air ysche 
Funktion wurde durch die Hauptwerte der 
Arcustangensfunktionen erklärt. Diese so defi- 
nierte Funktion erfüllt somit die an sie zu 

Bild 2. Scheibe aus dem unendlich langen Zylinder stellenden Forderungen mit den genannten 

Ausnahmen. Sie kann daher als Lösung des 

Differentialgleichungssystems bei verschwindenden Randwerten gelten. Diese nichtreguläre 

Lösung ist also mit Ausnahme der Singularitäten auf den Erzeugenden die Antwort auf die 

Frage nach den Eigenspannungen des Wärmespannungsproblems. Läßt man die unstetige 

Temperaturverteilung dadurch in die gleichmäßige übergehen, daß man @g,=r wählt, so ver- 
schwindet der Spannungszustand, wie es sein muß. 
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Bild 3. Schubspannung in der Trennfläche 
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Bild 5. Normalspannung Oy längs der auf Druck beanspruchten Seite von AB 
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Die gleiche Erledigung hätte die Frage gefunden, die Spannungen zu finden bei einem 
gleichförmig erwärmten Körper unter der Voraussetzung, daß dieser aus stetig anschließenden 
Sektoren mit verschiedener Wärmedehnung zusammengesetzt wurde. 


Zum Schluß sind die allgemeinen Ausdrücke des konzentrischen Kreisausschnittes für den 
speziellen Fall zurechtgemacht, daß die Dichte sich längs einer parallel zu einem Durchmesser 
und der Zylinderachse verlaufenden Diskontinuitätsfläche unstetig ändert. Die Bedingung für 


eine solche Fläche lautet 9, = 5 . Der solcherart erzeugte Spannungszustand wurde für einige 


Durchmesser numerisch gegeben für zz z =. 


Bild2 stellt eine Scheibe dar aus dem unendlich lang vorausgesetzten Zylinder, in 
der das fragliche Stück der Diskontinuitätsfläche durch Schraffur hervorgehoben wurde. Die 


Schubspannung längs des Radius = — Z also in der einen Hälfte der Trennfläche, sind 


in Bild 3 gezeichnet. Als Abszissenmaß wurde das Verhältnis r/r, gewählt, wobei r, der Kreis- 
radius der Grundfläche und r ein innerer Punkt auf dieser sind. Als Ordinaten ist die Schub- 
spannung 7 eingezeichnet, die längs dieses Teiles der Trennfläche auftritt. Die Schubspan- 
- nung geht stetig durch die trennende Ebene hindurch. Die Normalspannung o, wird in der 
Trennfläche null, während die Normalspannung o, vom Wert + 1/2 auf der einen Seite auf 
den Wert — 1/2 auf der anderen Seite springt. Beide Spannungen sind nicht eingezeichnet. 

In Bild 4 und 5 sind die Spannungen o,und o,, die längs der Linie AB des Bildes 2 auf- 
treten, gegen die Abszisse r/r, aufgetragen. Längs dieser Linien verschwindet die Schub- 
spannung r. Die Normalspannungen sind zugleich Hauptspannungen. 


An den Stellnr =r,„p = — z undr=r,9=- 5 hat keine der Größen einen de- 


finierten Wert, sondern hängt jeweils ab von der Richtung, in welcher man in die betref- 
fenden Punkte einläuft. Ist y der Einlaufswinkel gegen die Trennlinie gemessen, so lautet 
z.B. der Ausdruck für die Normalspannung o, dort: 


; 
Go 571 2sin?ysin 2y + 2arctg (cot y)}. 


| Anhang 
Es gilt I(r,9;y) mit v=g, und y= — 9, Wobei 
R=2(y)=r?+r?— 2rr,cos (+ Y), 
R?(a)=r?+r5+2rr,c0s op. 
gesetzt ist. 
Dann wird mit der unten angegebenen Festsetzung für 7’ 


F= IT, 9; 9) — Ir, 9; — 9) 


A-07 gMrio au) el SP 


und entsprechend 


Se; or r or 
usw. 
Anhangl 
ee 1 2 (nn — r)? — Art nn —r 7 
Ir, 0:9) = 16z j: r2 = — (0829 + ) larctg > cot nn 


En n—r P-+Y y2 —y2 
— arctg & 3% cot — ) = Er @ — nr) G — 9) cos 2(p+Y) 


+ sin (+ Y)+sin (#+2%)+sin2(P-+ y)In 


+ Rp) Zrart sin (P+Y)+ (+ 7)sin2(P+Yy) — (ri + 107) sin 3 (P+ y)} 
+ BR>(a) {2 77° sin (9) — (+ P) sin 2(P+Y)— ri tnr)sin gt | 
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Anhang 2: 


Ri - en 2 - = 


ı er rr—2rr? + r)sin2(P + Yy) (areie 1, Fr eot zz aretg nn " cot Ze) 


— 2ngcot PIP 21. cot PI* 
+(8r5—4r{n —r2}?)sin 2(p+ y) +% | = 


(mr + neo (nennen 


2 
+(4rt — 2318 — 2rör? + ricos2(p-+ y) 3 (3 — r?)(R2(n)— R2(p))+ (475 + {27 (13 — 72)? 
— 47%} cos2 (p+ Y))(Rea) ro{(rg + 72) cos (Pp+ m) —2ror} — Rt(p)ro{(r$ + 7?) cos (P + Y) 


—2nr)+2(—-3r+2rRr? + rt) (nr — y)sin2(@-+y) 


a Br+ 2 Het 2) (os (Pt Y)+ cos @+2Y) 
—- —3n+2r7°+r)sin2(P+ y) (R(y)2rn,sin(p+ Y)+ R(r) 2 ryr sin p) 


— r(rg — r?)? sin 2 (9 +Yy) (R*(p) 2 rofrg — 7} sin (+ Y)+ RA(n) 27, {r} — r?} sin p) 


—2r+4r° 


— R(y)2nrt (nn — r?)sin(p+yY) (—2r,sin (P-+Y) — Ar sin2(P +) 
= N sin3(p + y)) — sin (P+y) Rp) 4 rort(r? — 3 rör 


+2 15c0s (e+y))(— 2ror sin ee rer IF in3(@+Y) 


en... u 0052 (P-+Y) 


A une cos 3 (+ 9)) + R?(w)(r} — (en ie pre Y)— 6(r, zu 


Tor 


) 


u 


—+ R%(n) Argrt sin p(r — 3rgr — 2 rg cos (5 a0 siin2(+y)+ sin(9-+2 n) 


+ Rn) Ar} — r?) = - a “ cos 2(P+Y) — 


6 
Ad nagr2r) 


2n+r° 


j — —n+3rt 
+2 R4(r) r,r%(r} — r?) sin @ ze sin 2 ( Fr En 


sin (39+ 29) 


ei 


— RA) 2 r,r(r3 + r2)cos + 2un(2®, ur -cos2(e+y)+ 3% „a 008 39 + 24) 


— Ar(r 3 r3)8 cos 2(p+ y) (R*(p){r — 1008 (P + y)}— Re(m){r + r.cos 2) 


10 
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Anhang3: 

1 or, 95%) _ 

r or % 
zu + o+m 
er (8 rt — A{r$ + rt} cos 2(p + %Y)) (aretg Pe ” eotg . — arcig I aa T cotg 5 ) 

en 

+(4 5 +[2r (13 — 2)? — 4 rd}cos 2(P+ Y)) ae 

(n-+r)+ (n—r)2cot?TF 5 = 

Zr? tr 


- 


It 


es rue E 2(p — a) (rn — r}) cos2 (P+ y) 
(r + 10)2+ (nr)? cot? I 


sin 2 rl sin 2(p+ y)2(Rr)ir — nocos (+ y) 
— R:(r) {r+ r,cos p}) er (3rt—2 Hr? — ri) (sin (P+Yy)+ sin (P+2Y)) 


Da 2 nrsin +) sin2(p +») 
ar — 27 


+ sin 3 (p+ m) Rn) (—2nsin pr) sin2(P+Y) 


ni 


- 


ar 2 in 3P+y))+2nf- 2ror — (rn + r?) cos pr? Rn ( 


4 nn — sin (3 +29) + Pan) Ra), Sin ++ en 
Anhang4: 
170 7190,90) 
72 99° 
- 07 N — 8{n—r?}2) cos 2 (p+ y) (arctg (= Sr T eot - 


2 


\\+ (87 — 4{rg — r?}?)sin 2 (P+ Y)(R:(a) {rd — r} — R&(y) {ri —r}) 
— A(r-r2)2(y — re y)+(Art+ {2 (ri r?)? Art cos2(p+y))Lkilpiin- rt} rorsin (P+ y) 
+ Rı(a){rg — 1} or sing) + 0 — 12)? (sin (P+Y) + sin (pP +2 Y)) 


— aretg | ° — ! cot 
N tr 


+4 — 1?)? sin 2(p+ y)In ao Ar — rrcos2(p+Y)(R’ly)Aror sin (p+ Y) 


= IE sin 9) — 2 pr (rg — r?)? sin 2(p+ RW r+ r?) cos (+ Y) — 2 rer} 
+ Rıla){(r$ + 7?) cos +2 ror}) — Rey) 2 rgr? (rd — r2){(rö +72 +2rr, cos (P+y)) cos (P-+Y) 


—4 a 2r,r sin (P+Y) — sin2(p-+y) + Br sin3(p+ ) 


[2 . 5 
— AR) ra r)sinp + y)(- 2 Peos(p+ 2"; Ar Apr + ET 0083 (ey) 
+ Ro) en) (Bin pr TEE WET in s(a+ 7) 


+ R%(r) 2 vor? (1? ee Pr Due 55in2(P+%) 


+ sin &p+ 29) +4 Ra) )rorsir sin (2 cos 2(p+%) 
ne -00s(3+2y))+ Rrla)r. (r -)( in 2p+y)— in g+2n) 
Tor 


— 


Anhang: 
ll, 9 Y) _ 
or? Fr 
= ae rt—4{rt— Ir} 60s2(p+y))(arctg I cotg vn —arctg Br cotg.? z a) 
27, ct FIR Dr cr IR A 


+2(8 7 —4r{rg + rttcos2(p+Yy))| - = ae TE 
j (n+r +{n—r) or (r,+r)?+ (rn, — r)?cot? a en 


A r,cot Br I( + r) — (nn — r) cot? PEN 
+ (4°+ 27° {(13— 72)? — 47} cos2(p+Y)) a Ten 


(ro Fer> Fin root? e4r pi 


A er let N u Er | 
ee +2 —r)(3r5+ 7%) cos2(p-+ Y) 
— 2 (r+3r5)sin 2(p-+y)In Pu +4r(rg— rt) sin 2(P+ Y) (2 R?(y) {r—ro cos(p+ y)} | 
— 2Re(a){r + nacosg}) — (rt +27) (in (Pr y)+ sin (@+2Y) 

— 27? (n —r?)sin 2(9+ Y) (R(y){2rorcos (Pp+ y)—r?+nn(1—2c0s:(p+Y))} 


+ Rıa){2rgr cos +1? — 1, (1 — 2cos? p)} 


127° +6r$ 


, Ar. 
+ in 


es 


sin 3 (p+ 9) 


f \ . . 6 
— Arte (2) (1 00s®(p4p)) Rsp)(— 2rrosin (p+y)— I sin 2 (P+Y) 
0 


ed ? 
+ BEE sin 3(g+ y))+ An {—2r0r+ 4 )eos (pt y)} Ri) —2rosin (P+v) 


Tor 
4r—2r) . Tor. 
a, sin2e+y)+ —— sin 3 (e+m) 
AR TgT 
6 6 4 „4 
— Artrt (3 — r?) (1— cos?) R®(r) in sin2(9+%)+ wer sin®9+2 ) 
0 0 


45 — 28 . 34 pi 
er roraf 2ror (r2 Hr?) cos a} R4 m sin 2 (p- Y)- — = sin (3 + „) 
0 f) 


DEE aan: 
Far sin (3p+2Y)) 
0 


127° +60 . 
+ Aral sin 2 (p+y)+ 
0 
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Fehlerabschätzungen zur Störungsrechnung für lineare N 
Eigenwertprobleme'beigewöhnlichen Differentialgleichungen” 
| Von J. Schröder in Hannover 


Ausgehend von einem an anderer Stelle!) bewiesenen Satz über die Störung des Spektrums von Ope- 
ratoren eines Hilbertschen Raumes werden Fehlerabschätzungen für lineare Eigenwertprobleme bei gewöhn- 
lichen Differentialgleichungen hergeleitet. An Beispielen wird die Güle der Abschälzungen geprüft um 
die Größe der zulässigen Störungen veranschaulicht. a 


The author has revently proved a theorem!) on the perturbation of ihe spectrum of operators in a Hübert 
space. By use of that iheorem, estimations of the error for linear eigenvalue problems of ordinary diffe- 
rential equations are deduced. The goodness of the estimations is evamined, and ihe amount of ihe per- 
missible perturbalions is demonstrated, by evamples. 


Prenant pour point de depart un th6orömet), prouve & un aulre endroit, "sur la perturbation du spectre 
d’operateurs d’un espace d’Hilbert, l’auteur deduit des estimations d’erreurs pour des probl&mes lin£aires 
de valeur-propre ü des &quations differentielles ordinaires. Le qualit£ des estimations est examine par quel- 
ques examples et le grandeur des permissable perturbations est monire. 


OrTupas.aach OT Ipeskze!) NOKA3aHHOH TeOpeMbI O BO3MYINEHHH CIISKTPA ONSPATOPOB THJIB- 
ÖeproBa IPOCTPAHCTBA, B PA6OTE BEIBONATCH OMCHKH OIMHOOK WIM JIHHEÄHBIX IIPOd.IeM COScT- 
BEHHBIX SHAYEHHÄ IIPOCTEHIX AuhhepeHnuansHsx ypasHeHunü. Ilposepaerca kayecTBo ONTeHOX 
H NOKA3EIBAETCH, HA IPUMepaX, YHCIO NOLYCTHMBIX BO3MYITeHuK. 


Es kommt vor, daß man eine gegebene Eigenwertaufgabe Ap=4p (Beispiel: 
— 9" +g(2) 9 = Ap für Funktionen p(z) mit 9(0) = p(1) = 0) nicht in geschlossener Form 
lösen kann, für ein ‚„‚benachbartes“ Problem A,p =Ap (hier etwa — eg" =/9) Lösungen A,, 
9, aber ohne Schwierigkeiten anzugeben vermag. Man kann dann versuchen, die gestörte 
Aufgabe A(e)p=Ap mit Ale)p—= A,p+teA,p und A,py=ApP—-Ap (- p’+teEeg(e)P 
=/_), welche für e= 1 in das Ausgangsproblem übergeht, durch Potenzreihenansatz in 

A=M)=h+Eehteiht..., 9=pe)=MtEep t®Mp+..-. 

zu lösen ?2). Die Möglichkeit einer solchen Entwicklung von A und 9 wurde zuerst ausführlich 
von F. Rellich?) allgemein für gewisse Klassen von abstrakten Operatoren A (e) unter- 
sucht. F. Rellich gab auch Fehlerabschätzungen für die nach dem n-ten Gliede ab- 
gebrochenen Reihen an. Andere Autoren führten diese Untersuchungen fort ?). Die Fehler- 
abschätzungen wurden vom Verfasser an anderer Stelle !) so weit verbessert, daß sie für die 
Anwendung in vielen Fällen brauchbar scheinen. Die dort erhaltenen Ergebnisse sollen hier 
zur Gewinnung von Aussagen für gewöhnliche Differentialgleichungen benutzt werden >). 
Um die schematische Anwendung der Abschätzungen möglichst einfach zu machen, ist im 
Abschnitt 3 der Rechengang in seiner am häufigsten vorkommenden Gestalt beschrieben. — 
Im folgenden bedeute e einen komplexen Parameter, den Störungsparameter. Setzt man in 
einem der betrachteten Störungsprobleme e=0, so entsteht das zugehörige ungestörte Problem. 


1. Die Ergebnisse für abstrakte Operatoren 


Es seien u, v, 9... die Elemente eines Hilbertschen Raumes 9. (u, v) bedeute das 


innere Produkt zweier Elemente u, v ; ||u]| = Y(w, %). Der folgende Satz wird hier ohne 
Beweis genannt. Die darin vorkommenden Zahlen a, y, o und Funktionen f,(e) erhält man 
unter Benutzung der in den Voraussetzungen erklärten Größen Q(&,n), ö und /, sowie der 
- Zahl A, = (A, 99 9,) aus den Gln. (5). 


Satz 1. Es seien A, und A, lineare Operatoren des Hilbertschen Raumes $ mit 
dem gleichen Definitionsbereich Y. 


*) Diese Arbeit stellt im wesentlichen einen Teil der Dissertation des Verfassers dar, Fehlerabschätzungen 
zur Störungsrechnung bei linearen Eigenwertproblemen, Diss. Hannover 1952, Ref. Prof. Dr. L. Collatz und 
Prof. Dr.-Ing. W. Quade. Einerster Teilerschien an anderer Stelle, s. 1). Gegenüber der unter 1) genannten 
Arbeit ist hier die Bezeichnungsweise geringfügig abgeändert. ; 


ı) J. Schröder, Fehlerabschätzungen zur Störungsrechnung bei linearen Eigenwertproblemen mit 
Operatoren eines Hilbertschen Raumes, Math. Nachr. 10 (1953), S. 113— 128. 

?2) Eine ausführliche Darstellung des Verfahrens der Störungsrechnung für Differentialgleichungen findet 
man bei L. Collatz, Eigenwertaufgaben mit technischen Anwendungen, Leipzig 1949, S, 378 ff. S, auch: 
R. Courant — D. Hilbert, Meth. d. math. Phys. I, 2. Aufl., Berlin 1931, 296—300; W. Meyer zur 
Gapellen, Meth. z. angen. Lösg. v. Eigenwertpr. m. Anwendg. auf Schwingungsprobl., Ann. Phys. (5. Folge) 
8 (1931), S. 297—352 und Genäh. Berechn. v. Eigenwert., Ing. Arch. 10 (1939), 167 — 174. 

°) Literaturangsben über die Störungstheorie findet man am unter !) a. O. 

4) Insbesondere B. v. Sz.-Nagy. 


5) Wir beschränken uns hier auf die Störung einfacher Eigenwerte, am unter !) a. O. wurden auch mehr- 
fache behandelt. 
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a) Der Operator 4A, sei hermitesch ®) und besitze ein in 9 vollständiges Orthonormal- 
system von Eigenelementen. 

b) Es gebe eine für nichtnegative Werte der reellen Veränderlichen &, n erklärte, ein- 
deutige und stetige Funktion Q(£, n), derart daß 


KENSOErn) für ES n<r......2..(a) 
und 
Kakan)=ad(E,n) N) 
für jede nichtnegative Konstante a ist und für alle Elemente # aus U 
IAul sQCllall Au) - - 2.2.2.0...) 
gilt. 


c) Der Operator A, besitze den einfachen Eigenwert A, und es gebe eine positive Kon- 
stante ö, derart daß im Intervall A, — 6<A<A,+ 6 außer A=4, kein Eigenwert von A, 
liegt. 9, sei ein zu A, gehöriges, normiertes Eigenelement; A, 9, = Ay 9 » ||Yol| = 1. 


d) Es sei 
AU (Aa + 5) « 
und 


RoS(a— +2 +e4,) fürs Aelı<og 


s\—1 
wenn Ä die beschränkte Fortsetzung des beschränkten Operators (4 — + a auf ganz 9 


bedeutet. 


Dann gilt: Für |e|< o besitzt der Operator A(e) = A,+ e 4, einen regulären Eigenwert 
Ale) =Ag +&/ı +89, +--- und ein zugehöriges reguläres?) Eigenelement $(e)= 
+E9+89+::-, für das (p(e),9) =1lund A,ole)= Am tEehpyı tr AP+ 


Die Störungsglieder A, und 9, (aus W) sind die eindeutigen Lösungen der Rekursionsformeln 
n—1 
An = (Aı Pn-ı; 90); (Ad — A) 9 = — Aıfn-ı+ 2 An Pı » (9; 9) = Nee 
Für |e|< e gelten die Fehlerabschätzungen 


ko Lern <dakliae, | Zen) she), 


ERGO -2 a e) <(ö+ AN): nr=L23, 


Zusatz:1. Unter den Voraussetzungen a) bis c) des Satzes 1 existiert für |e|< o der 
$ 1 
Operator (4 — + 7 .- 4) : 
2. Wenn auch A, hermetisch ist, gilt sogar 


A) Seh söYylelfn(e)- 


2. Anwendung auf Eigenwertprobleme bei gewöhnlichen Differentialgleiehungen 


a) In diesem Abschnitt werden Eigenwertaufgaben bei gewöhnlichen Differential- 
gleichungen als Probleme in einem Hilbertschen Raum gedeutet, und es wird gezeigt, daß 
für diese dann die Voraussetzungen des Satzes 1 erfüllt sind. Die daraus folgenden Ergebnisse 
findet man unter Vermeidung der Operatorenschreibweise zusammengefaßt. Wir nennen zwei 
Eigenwertaufgaben identisch, wenn sie die-gleichen Eigenwerte und -funktionen besitzen. 
x sei eine reelle Variable. Aussagen über Eigenschaften von Funktionen dieser Variablen 
(Stetigkeit, Integrierbarkeit usw.) sollen sich nur auf das Intervall 0 sz s1 beziehen. 


6) Ein (dicht definierter, linearer) Operator A heißt hermitesch, wenn für jedes Paar von Elementen u, ® 
aus seinem Definitionsbereich (Au, v) = (u, Av). Pe 
?) Ein Element o(e) beißt regulär für |e|</ go, wenn es von e unabhängige Elemente 99: ®1, Pa, - - - gibt, 


n 
derart ußo()= m te ++... fürle<o, d. nel | Ip(e) e2® %||=0. 
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u, v, 9, seien quadratisch integrable Funktionen der Veränderlichen &, die gegebenenfalls | 
auch noch von e abhängen können), . 
b) Wir betrachten Eigenwertprobleme der folgenden Form: 


M,(o) ke Mo) re 
vVo=20 „ela nm an 
mit 


m 3 2zm—1 EN 
M,(p) = Mı(p(a)) = iz 1 (sta) PM (a))9,  Mılp) = 1em(2) M(p)+ I 9,(2) 99 (a); 


U, = I lau O)+ Br)  m>0). 


Die darin vorkommenden Funktionen 9,, 9, q, und Konstanten &, P,, seien gegeben und 
haben folgende Eigenschaften: die Funktionen p,(x) sind reellwertig und besitzen stetige Ab- 
leitungen bis zur v-ten Ordnung, pm (2)> 0, 9(&) ist stetig und positiv, die Funktionen 9 (x) 
sind stückweise stetig, die Konstanten a,, und ß,, sind reell, bei festem u verschwinden nicht 
alle Konstanten a, und Ba»(vr = 0,1, - -,2m — 1). a 
Gesucht sind Zahlen A (Eigenwerte) und Vergleichsfunktionen ® (Eigenfunktionen), die 

der gegebenen Differentialgleichung (3a) genügen. Als Vergleichsfunktionen seien dabei die- 
jenigen Funktionen bezeichnet, die stetige Ableitungen bis zur (2m — 1)-ten Ordnung und eine 
stückweise stetige Qm-te Ableitung besitzen und den gegebenen Randbedingungen (3b) ge- 
nügen. 

: c) Wir benutzen im folgenden ferner die Voraussetzung, daß die Randbedingungen (3b) 
voneinander linear unabhängig sind und für jedes Paar von Vergleichsfunktionen %, v 


J Mut 500 = fu ifo) de EEE 


gilt ?). Diese Voraussetzung ist bei allen Aufgaben erfüllt, für die in den Tabellen I und II 
eine Abschätzungsfunktion angegeben ist. 

d) Alle Probleme (3), für welche auch die unter c) genannte Voraussetzung erfüllt ist, 
lassen sich als Eigenwertaufgaben (4,+ e4.)p=Apin einem Hilbertschen Raum mit 
einem Operator A, deuten, welcher der Voraussetzung a) des Satzes 1 genügt. Nur für einige 
spezielle Typen jedoch werden hier Abschätzungsfunktionen Q(&, n) angegeben, mit denen 
(1) und (2)gilt.. (Diese Problemtypen und die zugehörigen Funktionen Q sind in den Tabellen I 
und Il aufgeführt.) Es wurde aber versucht, die praktisch wichtigsten Aufgaben zu erfassen. 
Als ungestörte Differentialgleichungen wurden in der Hauptsache solche mit konstanten 
Koeffizienten — nämlich — ’ = 9 und p!Y =Ap — gewählt, da diese sich in geschlossener 
Form lösen lassen. Für die wichtigsten Randbedingungen sind die erhaltenen Abschätzungs- 
funktionen sehr einfach. 

e) Zur Ermittelung der Abschätzungsfunktionen Q(£&, n) aus den Tabellen I und II sei 
folgendes gesagt. Die in den Tabellen benutzten Bezeichnungen sind unten in Tabelle I er- 
klärt. In den mit (A), (B), (C) bezeichneten Zeilen werden durch die Angabe spezieller Formen 
für M,(p) und M,(p)und zusätzlicher Voraussetzungen über die in der Differentialgleichung 
auftretenden Funktionen drei (sich z. T. überschneidende) Typen von Differentialgleichungen 
(3a) charakterisiert. Liegt eine Differentialgleichung von einem der angegebenen Typen vor, 
so kann man die darunter stehende Funktion Q(£, n) verwenden, falls die gegebenen Randbe- 
dingungen eine der jeweils in der mit ‚Randbedingungen‘ überschriebenen Spalte genannten 
Formen haben und in den Fällen (A) und (C) außerdem die in der betreffenden Funktion ge- 
gebenenfalls vorkommenden Konstanten M, N bzw. K, Lden Tabellen zu entnehmen sind. Diese 
Konstanten erhält man jeweils unter den in der nächsten Spalte genannten zusätzlichen Voraus- 
setzungen (über die in den Randbedingungen auftretenden Konstanten) aus den beiden letzten 
Spalten. Ist der Faktor c, eines Ausdruckes gleich Null, welcher eine solche, noch zu be- 
stimmende Konstante enthält (d. h. ist q, (x) = 0), so setze man das ganze Produkt gleich Null. 


®) Im folgenden wird pim allgemeinen von z und gabhängen. Je nach Bedarf kennzeichnen wir sowohl die 
Abhängigkeit von & als auch die von e durch runde Klammern [9 (8), # ()]. Verwechselungen dürften kaum zu 
befürchten sein. pP). bedeutet immer die v-te Ableitung nach x, (0) bzw. (1) ist @ (2) an der Stelle x =0 
bzw. © = 1, entsprechend 9 (0), ” (1) usw. 

°) Mit dem Querstrich wird der Übergang zur konjugiert komplexen Funktion gekennzeichnet. Sind alle 
Funktionen 9, reellwertig und interessiert man sich auch nur für reelle Werte & (etwa e = 1), so braucht man auch 
nur reellwertige Funktionen v,v,@,...zu betrachten. Dann kann der Querstrich hier wie im folgenden weg- 
fallen. — Ist (4) erfüllt, so ist das ungestörte Problem (3) selbstadjungiert in dem etwa bei L. Collatz a.a.O 
S. 48 ff. erklärten Sinne. Die Bezeichnung ‚,‚Vergleichsfunktion“ findet man an der gleichen Stelle. 


N 


Dh a a a an 6 
j i Y KıOPM 


er 


- 6; i Sat 22: Die Aufgabe (3) habe für e= a den einfachen Eigenwert Ag» Po sei eine auge“ 


er 
 hörige, reellwertige, durch [ gleol? dz = 1 normierte Eigenfunktion. 
ö 


zugehörige Eigenfunktion rel für die J g9P(E)Y9, de=1. A(e) und o(e) lassen sich in Potenz- 


reihen nache 


4() A EAg Ehe re) le) = at Em+&®p+: 


(mit Vergleichsfunktionen Pn) entwickeln. Die teen: der Eigenfunktion gilt bei festeme 


bezüglich x im Sinne gleichmäßiger Konvergenz. Die Zahlen A, und Vergleichsfunktionen o,sind 
die Lösungen der folgenden Rekursionsformeln. Für je] <o gelten die unten angegebenen 
Fehlerabschätzungen. 


Rekursionsformeln n=1,23,3,.-.): ® 
| 1 
In = [ Mı(9n-) 90 de; 
n—1 
M, (Pr) — Ag pn = — Mı (Pn—ı) +9 hr per 


 Du(@n) =-0(w=1,2,::,2m), Jorında =0. 


Bezeichnungen:Es sei 


Q(£, n) eine bestimmte, für das a Problem in den Tabellen Iund II angegebene 
Abschätzungsfunktion (vgl. 2e); 


ö eine positive Konstante, derart daß im Intervall u — d<A<A, +6 außer A — As 


kein Eigenwert des ungestörten Problems liegt (man wähle ö möglichst groß); 


UN. A, „ln. u ee 
a FON + KOMITEE Fi 
ee  eluiye, = (at r)enıta 2 nn, Basel... °6); 
Hefe). HO-fO— 2 kt. U DEREN, DNA ADUS ES (5d); 


@(x, y) die Green sche Funktion des ungestörten Problems (d.h. des Differentialaus- 
drucks M,(g) zu den Randbedingungen (3b)), falls diese existiert !?). 


10%) Auf Grund der Ungleichung r, <y zo kann man zur Vereinfachung r3, durch » und y seinerseits durch 
& ersetzen, wodurch jedoch u. U. wesentlich an Genauigkeit verloren geht. 

11) Eine insbesondere bei Berechnung einer größeren Anzahl von Störungsgliedern genauere Werte liefernde 
Funktion ist 


oe (für eo #0) [r-"Z usr ale) 


12) Für die wichtigsten Probleme 2. und 4. Ordnung ist die Greensche Funktion bei L. Collatz 
a.a. 0: 8.425 f, we 


_ Dann besitzt die Aufgabe für |e| <e (o durch (5a) erklärt) einen Eigenwert A(e) und eine 
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Fehlerabschätzungen®) n=1, 2,3, -»): 


Abschätzung unter der jeweils links stehenden zusätz- 


zusätzliche Voraussetzungen lichen Voraussetzung 


<öalelfn (e) 


n 
Meyer, 
=0 


v 


ei 
| 


1 Inn 
V; alle) — BE dx < fn (e) 


2 Für jedes Paar von Vergleichsfunktionen u, v gilt 


3 006 KR 
SM, (u) ode = f uM,() de = öylelfn(e) 
[) 0 


n 
Ale) — 3 eh 
v—( 


ar, 
p(e) — Zen <\@(z, x) Vö + |Ao| fn (e) 


3 | alle Eigenwerte der ungestörten Aufgabe sind positiv 


—ı1 1 
4 ) = (ist kein Eigenwert der ungestörten Aufgabe b (e) - Zeo|< V S 9(y) I@(, y)l? dy(ö + ]Ao|) fn(e) 
v=0 0 


Tabelle L —Differentialgleichungen 2. Ordnung 


(4) M,(p) = — 9”, Mı (pP) = — 959” +919 + 909, zus. Vor.:g(@)=]1 
gen) Ma + YN)n + VMine+e.E 
Randbedingungen zusätzl. Vorauss. | M | N 
ab>0, cd>0 1! [) 
€ 
b=1 a>o— T 0 
ap (0) -bp()=0| d=1 
a 
cp (lW)+Ädp()=0 ee n 
a a 
A: am (+4 +5) 
ledl—cd[, b® 
A=ad+tbc-+bd A=+0 0 32 (a +ab +5) 
oder 
I: ad? —abcd +b?2c2 
6 . 342 
1 
D=d=0 0 = 
M=ap(l) —_ 1 0 
I 1 
a9 (0) =P() = 0 im 
(B) N) = —- (MP) +99; Mı @) =9%M,(P) + 919’ + 909, zus. Vor.: 9, (2) 0 
Randbedingungen Q(&n) 
a9’ (0) —.bo (0) = 0 it ab 0 ’ 
er t+dom=0 " ca>o nt u me+ te 
r*» 
pP (0) =ap(l) 
h t ß 
ap (0) 9° (0) = pl) 9° (1) N 


u u FF PD EEE BEE GE TE re an u 
Bezeichnungen in den Tabellen I, II 


% sind Konstanten, mit denen |g, («) | zo , 
A ist der Betrag des dem Betrage nach kleinsten Eigenwertes der ungestörten Aufgabe, 
x» =|/A. 


———— 


12) L.Gollatz weista.a. O. S.384 ff. darauf hin, daß A +Ad i i j 
‚0. S. ; : 0 en mit der Funktion @, gebildeten Ra y- 
leighschen Quotienten der gestörten Aufgabe für = 1 darstellt. "Ist auch diese Aufgabe De 


folgt daraus, falls A (l) der dem Betrage nach kleinste Ei i i 
a g einste Eigenwert dieser Aufgabe ist, [A,(l)| < |]A+A,|, also 
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Beweisdes Satzes 2. Essei 9 der Hilbertsche Raum der quadratisch integrablen 
1 
Funktionen %, v, 9, : : - mit deminneren Produkt (u, v) = ) guwödz. A,und A, seien die durch 
1 ie 
Au = Er M,(w), Aıu = & M,(u) in der Linearmannigfaltigkeit X der Vergleichsfunktionen 


definierten Operatoren. Das Problem (3)ist dann mit der Aufgabe (A, + eA,) 9 =/p (paus A) 
identisch. Wir prüfen zuerst, ob die Voraussetzungen des Satzes 1 erfüllt sind. 


Tabelle IIT—Differentialgleichungen 4 Ordnung: 


(©) MO=W, MO-u tar’ +nP’tnP Hp zus. Vo: geo)=1 
ge, Manta EM + alnE+nQK, Em +cE£ 
T 
RK + mE falls K #0 Verla falls L = 0 
mit Qı (K, 5 n) — 1 3 f ’ Qs (L, ER n) — 3 N 
nt & falls X=0 nt falls L = 0 
| Randbedingungen | zus. Vor. | K L 
1 | Periodizität — 0 0 
(0) = 0,09”) —be)=0 ZVe- eb Habe 
2 Per MO, La>ol FVac@a Fb) We ta) 
a,b,c,d>0 ee 1 17 de 
A=3(ad +be) +bd en a 
"=," W-IpgW=0| zo ae, le 
3 19(0)=0,cp’(N)+dp(l)=0 b(c + d) c+dY15 
a,b,,d 0 Na 0 0 
Fr i TER Ir 
_ lee (120 +b) (Adc+d) nicht berechnet 
(0) =0,ap" bp (l)=0 ' 4ye« + ab + 20 
4 (0) =0,cp’(M)+Ap(l)=0| °= Va ia +5) 
a,b,,d>0 in a 
A=Matetd tb(ctd -/ 360° + 1508 + 358° 
| 2 1/13 
N) 0 | ee 
WC 35/3 oder 0,119 
Rt ap’’’(0) +bo (0) = 0 Störungsprobleme mit derartigen ‚„symmetrischen‘‘ 
co" (0) — dep (0) = 0 b>0 Randbedingungen lassen sich, falls Hla)=m(l — a) 
vr] NE ist, auf Probleme mit Randbedingungen vom Typ 2 
a) ve DPON- und 3 zurückführen 
el) +4p()=0 e1 
EN b=d=0 nicht berechnet 0 
+90) =0 | 3 R 
Be An is en = ac+ta*rct=( 0 nicht berechnet 
6 a*9”’ (1) — b*o(l) =0 on nn 
EN FEAHN=t bd-+b*d*=0| nicht berechnet 0 
alle Koeffiz. > 0 


1. Der Operator A,ist hermitesch, denn es gilt infolge (4) für jedes Paar von Vergleichs- 
funktionen u, v 


(A,u, v) = ] M,(u)ü da = uH,0) da — (u, Ay). 
ö 0 


Außerdem besitzt das ungestörte Problem, d.h. der Operator A,, nach bekannten Sätzen 1) 
ein in 9 vollständiges Orthonormalsystem von Eigenfunktionen. Wie man ferner weiß, häufen 
sich die zugehörigen Eigenwerte im Endlichen nicht !*). Also gibt es sicher eine Konstante 6, 
derart daß im Intervall u — 6 <A<A, + 6 außer A— A, kein Eigenwert von A, liegt. Die 
Voraussetzungen a) und c) des Satzes 1 sind also erfüllt. 


14) S. etwa E. Kamke, Differentialgleichungen, Lösungsmethoden und Lösungen I, 4. Aufl., Leipzig 
1951, S. 208 £. 
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2. Es werde nun zunächst angenommen, daß auch die Voraussetzung b) mit der ange- 
gebenen Funktion Q(£, n)erfüllt ist, und unter dieser Annahme die Voraussetzung d) behandelt. 
Im folgenden sei Je|< e. Wie man aus den Tabellen abliest, gilt Q(&, n) < cgmn, wenn |gem| Scem- 
Daraus folgt, falls ggm # 0, 


c 
?m|1-+ € gem] Z Pm(l — lel com) > Pm(l — 0 02m) pm ee = m[1 — 7. 0, 
ur 
der Koeffizient der höchsten Ableitung im Differentialausdruck 
ö 
Mu) -Agutz gu+e M,(“) . . FR a a ya en et je nt (6) 


ist also von Null verschieden. Da außerdem nach dem Zusatz 1 zu Satz 1 die Reziproke des 


Operators (4 — + S+ € 4,) existiert, läßt sich die Greensche Funktion I'(z, y, e) des 


Differentialausdrucks (6) zu den Randbedingungen (3b) konstruieren. (4 — At —_ A) 


ist dann aber die Menge der stückweise stetigen Funktionen $’, denn einerseits ist jede 
Funktion v der Form 


0) (4 + 5 e Ay) u = „| — gut 2 gu = eM, (u) re 
(mit u aus A) stückweise stetig, und andererseits ist die Gleichung (7) bei beliebiger stückweise 
stetiger Funktion v durch die Vergleichsfunktion u = j T'(z, y, e) 9(y) v(y) dy erfüllt. Daraus 
folgt sofort A,AC< (4, — + 2) X, da jede Funklon v— A,u (mit u aus WA) stückweise 
stetig ist. 
Die Fortsetzung K des Operators (4, — + . auf ganz 9 ist durch 


1 
Ru=) (m 9. 0)o(y) uly) dy 
für alle quadratisch integralen Funktionen u erklärt, also ist jede Funktion v der Gestalt 
v = Ku stetig und damit für Je|< gauch KÄHCH’ = (4, — A +2 + cA,) 4. 


3. Die Herleitung der Abschätzungsfunktionen in den Tabellen kann hier nur an einem 
Beispiel erläutert werden. Wir benutzen dabei die in Tabelle Ierklärten Bezeichnungen. sei 
im folgenden eine beliebige Vergleichsfunktion. Von den erhaltenen Abschätzungsfunktionen 
8 n) zeigt man leicht, daß sie den Beziehungen (1) genügen. Wir beweisen hier nur, daß mit 


ihnen 
dıwl) Ss QAllell Lu) Free 
Beispiel: —Q’ +e[g(— P)+ q1P+ pP] = AP, p(0) — dp(0)=0, P(1)=0. 


Für b=#0 existiert die Greensche Funktion des ungestörten Problems @(x, y) = F +% 
1 

(für z<y). Damit ist dann u(«) = / @(z, y)M,(u(y)) dy. 
0 


gilt. 


Alle hier vorgeführten Abschätzungen gehen von der folgenden Ungleichung aus: 


dull Ss |gedowll + IgiwWll + Ilgowll S cz ||Asull + cı Il] + co Ill. 
Hierin ist nur noch |[w‘|| abzuschätzen, was auf verschiedene Weise geschehen kann. 
a) Ist b=+0, so hat man 


I? = [ 


insgesamt also (8) mit 


|2 


[a0 y) Dorn du < | | 1592 y)a2ay Au? = 4 lau, 


nt 1 
Oma c\ N 3) 08 


—— ah 


Tr Fer, H an #9 r“ 4 R 
5 Br > 2 A Bi 
Wuu)<— AR + + Il4sull Im! Let 
Ei lAsullell. 
Im Falle b<0O a man weiter ab 


Ni (0, ) Alu) in 


OR = < [ 60. May lau = v  Asull. 


| Für b =+ 0 erhält man wegen ||| jü Maul also insgesamt (8) mit 


en = (at 7" +08. 
ur) 1 b> 0, so gilt (8) wegen (9) mit 
En) = en+ a /nE+ 08 
4. Falls für jedes Paar w, v von Vergleichsfunktionen ? m, (u)ö de = af u M,(v) % gilt, ist 


auch A, hermitesch, dann ist also auch die in Zusatz 2 zu Satz 1 genannte Voraussetzung er- 
füllt. 

5. Aus dem bisherigen Beweis folgen mit Satz 1 und Zusatz 2 die PeBaUpSUNBen des 
Satzes 2 bis auf die im folgenden bewiesenen, daß die Reihe 9 +e.91 +89, +: gleich- 
mäßig gegen p(e) konvergiert und die rungen 3 und 4richtig sind. 

Existiert @(&, y), so hat man 


Wal = | ta Due) du </Tomw Marla. - 00) 


Ist A=0 h-facher Eigenwert der ungestörten Aufgabe zu den reellwertigen, orthonormalen 


Eigenfunktionen y (z)(v = 1,2, : - -, h) ((y’, y*)=ö,,) und @* (x, y) die zugehörige, reellwer- 
I 
tige, symmetrische, durch /[ G*(z, y)g(a) y(a)de=0 (=1,2,::-,h) festgelegte verall- 


ö 
gemeinerte Greensche Funktion ®), so gilt an Stelle von (10) 


Wal < 2 Ku wr)l Ir @)l--|] 0°C Wlntutı) au 


< 3 Ww(alllilt V} Iswle*@, YlraylAsull 


Infolge dieser Ungleichungen er wegen der Beschränktheit der Funktionen g9(x), @(x, y) 
G*(z, y), y’(x) sowie der im Satz 1 genannten Konvergenz der Reihen 9+ e91+ &ps+ : 
und Yast EA, Pa+ &AsPa+ : : gegen p(e) bzw. A,p(e)im Sinne der Metrik desHilbert- 
schen Raumes 9, konvergiert die Reihe + e91 + &Pa + : : : beifestem e auch gleichmäßig 
in x gegen o(e). 

Sind alle (ihrer Vielfachheit entsprechend oft gezählten) Eigenwerte #/(v=1,2,3, + -) 
der ungestörten Aufgabe positiv — 9’(z) seien zugehörige, orthonormale Eigenfunktionen, 
(9, p#) = 6,, — so erhält man Aiet Benutzung des Mercer schen Satzes !®) 


Warum) SEE. Earl pe = Ga 2) (An, W) 
se@ollaulllll. . . - - (m. 


Aus den Tnsleichungen (10), (11) und den im Satz 1 angegebenen EN 
folgen die Fehlerabschätzungen 3 und 4. 
Damit ist Satz 2 vollständig bewiesen. 


3. Rechengang und Ergebnis im einfachsten Fall !”) 


Im allgemeinen ist eine Eigenwertaufgabe nicht bereits in der Form eines Störungspro- 
blems (3) mit einem Parameter e gegeben, sondern man muß, um die hier genannten Ergebnisse 


anwenden zu können, erst versuchen, ein Störungsproblem (3) zu ermitteln, welches für e=1 


15) Vol. E. Kamke a.a. O. S. 190 ff. 
16) Vgl. R. Courant —D. Hilbert a.a. O0. 8. 117. 
. 17) Oft wird man von den Störungsgliedern nur A, berechnen und auch nur den Eigenwert abschätzen. Für 


“ diesen einfachsten, zugleich aber praktisch wichtigsten Fall wird hier der Rechengang beschrieben. 


TE 


stelle fest, ob A, ein einfacher Eigenwert ist !%) u 


1 Peg, - 
Au= | Mılpo) 90 de. a a 
3. Man ermittele die oben unter ‚Bezeichnungen ‘“ angegebenen Hilfsgrößen 
lich (5b), wobei man überall e= 1 setze, und stelle fest, ob g > list !%). (Zur Ern 1% 
Q(£&, n) hat man die Funktionen g,(z) abzuschätzen und gegebenenfalls weitere in Q vorl om- 
mende Zahlen zu berechnen.) AT 
4. Ist A, ein einfacher Eigenwert der ungestörten Aufgabe und o> 1, so besitzt die Auf- 
gabe (3) für e = l einen Eigenwert A, für den A — (A, +A,)| S da f(1). Ist für alle Vergleichs- 
funktionen u, ® ö P 2 . 


1 1 
J M;(u) öde — [u M,(v) da 
0 0 = 
(das wird oft der Fall sein und ist mittels partieller Integration festzustellen), so gilt sogar 
A (+ A) sy)". - 
4. Zwei Beispiele | 

Beispiel 1 (bei diesem Beispiel werden die Näherungen für die Eigenwerte mit den 
exakten Werten verglichen): 
(et? p) =i*p, O)=0, yU=0, a=%. I 


Als Störungsproblem, welches für & = lin diese Aufgabe übergeht, wählen wir 

M(9)+ eM,(p) = 49, p(0) = 0, p(1) = 0 mit M,(p) = — 9”, M,(p) = 9.(—-P")+q1P (12) 
et 1 ben 

und A=b4*, 8, =-b(a ta” 1, 1 = — Te: a+a-+ Du Es wurde — z? abge- 

spalten und bso bestimmt, daß c,ein Minimum wird. 

Die Aufgabe (12) ist von der unter 2b beschriebenen Form und speziell vom Typ (A) in 
Tabelle I. Wir können die dort angegebene Abschätzungsfunktion mit M=1und N = Over- = 
wenden. 

Eigenwerte und -funktionen des ungestörten Problems: 


” = 272,9 (2) = Yasinyvanv=1,2,3,--- 


2 
Störungsglied A, für A, = n?: = Ah — ı —= —0,002722 - - - 
‘ r 41 84 er 
Hilfsgrößen: = Am E&;n)=cen+cYn& 69=3n%, 
a& — 0,08620, yv —= 0,02685, r, = 0,000092, 0 — 5,01, 
(1) = 0,029387, k(1) = 0,094342, (1) = 0,029471, @(2, 2) = al — eo). i 


Der Maximalbetrag des Koeffizienten der zweiten bzw. ersten Ableitung im Differentialaus- 
druck M,(@) (c, bzw. c,) beträgt etwa 5 bzw. 10 % vom Betrage 1 des Koeffizienten der zweiten 
Ableitung in M,(p). Es handelt sich also um eine kleine Störung. Die Störungsrechnung er- 
scheint als geeignetes Mittel zur Lösung der Aufgabe. 

Fehlerabschätzungen für den Eigenwert (an der Stelle & = 1): Für jedes Paar von Ver- 


1 1 
gleichsfunktionen u,v gilt [ M,(u) öde = [wM,(v) da. Daher ist 
6 0 
A (A, + Anl <öyflt) = 0,0234. 
Das ergibt A* = 4149,02 + 9,86 (der genaue Wert ist A* — 4146,306 - - :). 
al at 
Et 


Ist dies nicht der Fall, so kommt man mit Hilfe der hier genannten Ergebnisse zu keiner Abschätzung. 


Für a = 10erhält man c, = ‚ 0=2,55, 4*—=1089,45+ 11,15, A*— 1086,72 .. j E 


18) 
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Benutzt man die Tatsache, daß A der kleinste Eigenwert der gestörten Aufgabe füre=1 
ist 19), so ergibt sich nach L. Collatz®2) A <A,+A, also 4139,16 <A* < 4149,02, 
)* — 4144,09 +4,93 [füra = 10: 1078,30 <A* < 1089,45, 4* = 1083,88 + 5,58]. 


Fehlerabschätzungen für die Eigenfunktion (e = 1, 9, = Y2 sinne): 


vl \p(&) — ol? da < 0,02%. 
0 . 
Da alle Eigenwerte A” positiv sind, gilt nach Abschätzung 3 ° 
\o(z) — 9, <Yx(1 — x) 0,1852 < 0,0926. 


Beispiel 2 (an diesem Beispiel wird veranschaulicht, wie sich der mit einer arts 
ten Funktion Q(£&, n) gebildete Wert o bei verschiedenartiger Störung verhält): 


— "tell Pe’) + Pf + pl =Ap, pl) = 0, P()=0. 


Be. 2 
Die Eigenwerte der ungestörten Aufgabe: A = De Tale, Wir 


unterscheiden auch die zugehörigen gestörten Eigenwerte und Zahlen &, y, o durch Anhän- 
gen eines oberen Index ». Die unten genannten Zahlen wurden mit Hilfe der Funktion 
Q=c;n + cıYnE& + 0% erhalten, r, wurde überall durch y ersetzt (vgl. !%)). 


Die folgende Tabelle gibt die Werte o” für drei spezielle Zahlentripel (c,, c,, c,) an 


Gain 7Oij. E60 Til mar Aa Su | ew>1) ung 
) 1 2 2(— 1) 
ui Absie 6e] 2atl.| ABI 42 — 3 . 
N ee: 2(— 1) Li = 
0 1 An | An ı 2 2y—1+ Aw +49—7 9 
1 1 Musa (v— 1)? 
For re ERS 2 
Eingezeichnet sind o!, A! (e), A? (s) und die Schranken 
Ao+8i + öyls|f(e) fü u = z £ 
Ist speziell g = — 1, q, =%= D, so erhält man #(e) = X (1 — e); für e = 1 entsteht 


ein Eigenwert (Null) unendlicher Vielfachheit & das Bild, in welchem o0,?°), Al(e), A2(e) 
und die Schranken A, + eA,+ öyelf(e) für 4, = u indeaschnet sind). 


192) Dies folgt hier aus den Sätzen der Störungstheorie, vgl. ®). In vielen Fällen wird man eine derartige 
Aussage aber auch auf Grund der physikalischen Bedeutung des gerade vorliegenden Problems machen können. 
20) Eine genauere, im Falle eines hermiteschen Störoperators gültige Formel (s. die unter *) genannte 


Diss.) führt auf o*=$ als untere Schranke für den Konvergenzradius von Al(e). 


Eingegangen am 11. März 1953. 


geschwindigkeit der Scheibe. Baer 
This paper deals with the turbulent boundary layer at the surface of a rotaling dis z 
the awis ofrotation. The thickness of ihe boundary layer and therefore also the moment of the rotat 
strongly on the ratio of the velocity of Ihe awial flow to the peripheral velocity ofthedisc. 
L’6coulement turbulent & un disc plat rotatoire, lequel s’6panche en möme temps dans la direc el 
tournant, est caleule, L’&paisseur de la couche entrainee par la rotation ei avec cela aussi le moment du di 
rotatoire dependent beaucoup de la proportion de la vitesse de l’&coulement parallele & l’axe de rolation ü 


vitesse peripherigue du disc. 2 we 
l u - ” > “ - Ka 2 h Bir 

B HacrToameä CTATbe BBIUHCAAETCH TYPÖY!IeEHTHOE TeyeHHe OKOMIO BPAIHAMINETOCH ILIOCKOTO 
ANCKa, ONHOBPEMEHHO HATEKAEeMOTO B HAIIPABJIEHHH OcH BpareHua. Tormmma c10A, 3aXBayaH- 
HOTO BPAINaATeJIBHBIM ABHIKeHNEM, & BMecTe C TeM HU MOMEHT BpäammeHuf, B 3HA4YHTENBHOH CTe 


IMeHH 3ABHCHAT OT OTHOMEHHA CKOPOCTH HATEKAHHA K OKPY:KHOH. cIOPOCTH AHCKaA. Ir? ,t 


1. Einleitung nr 7 3 

Die Strömungsvorgänge an einem angeströmten Körper, welcher eine Drehbewegung 
ausführt, spielen bei Strömungsmaschinen und Geschoßkörpern eine wichtige Rolle. Über die 
Bedeutung dieses Einflusses wurde bereits in einer Arbeit von H. Schlichting undE. 
Truckenbrodt [1] berichtet, die sich mit der Berechnung der laminaren Reibungsschicht 
an einer angeblasenen rotierenden Scheibe 
befaßt !). Die vorliegende Arbeit stellt die 
Erweiterung der in [1] behandelten Aufgabe 
auf den Fall turbulenter Strömung dar, der 
technisch von größerer Bedeutung ist. Der 
Sonderfall einer in ruhender Flüssigkeit um- 
laufenden Scheibe ist schon frühzeitig von 
Th. von Kärmän [2] und später von 
Ss. Goldstein, vgl. [3], S. 404 behandelt 
worden. 


2. Grundgleichungen 


Für die Berechnung der Strömung legen 
wir Zylinderkoordinaten nach Bild 1 zu- 
grunde. Die Scheibe liege in der Ebene 2 = 0; 
sie drehe sich um die z-Achse mit der Winkel- 
geschwindigkeit w. In der Ebene der Scheibe 
seien die Koordinaten rund @. Die Geschwin- 
digkeitskomponenten seien # inradialer Rich- 
tung und vin Umfangsrichtung. Entsprechend 
seien 7, und 7, die Wandschubspannungen in 
den beiden Richtungen. 

Nehmen wir an, daß sich die Scheibe in 
radialer Richtung bisins Unendliche erstreckt, 
dann gilt für. die potentialtheoretische Ge- 
schwindigkeitsverteilung der in Richtung der 
negativen 2-Achse angeströmten Scheibe (ro- 


Bild 1. Erläuterungsskizze, 


Koordinaten Geschwindigkeiten 2 & “ 
ee U tationssymmetrische Staupunktsströmung) 
azimutal 9 © VW 
axial 2 w, W D.= ar; Va; W= —2az ep 


*) Die vorliegende Arbeit stellt den ersten Teil einer von der Fakultät für Maschinenwesen der Technischen 
Hochschule Braunschweig angenommenen Habilitationsschrift mit dem Titel ‚Die Strömung um rotierende 
Drehkörper bei axialer Anströmung‘‘ dar. (Berichter Prof. Dr.H. Schliehting, Mitberichter Prof. Dr. H. 
Blenk). Der zweite Teil mit dem Titel „Ein Quadraturverfahren zur Berechnung der Reibungsschicht an axial 
angeströmten rotierenden Drehkörpern‘ ist inzwischen erschienen, [5]. 

Diese Untersuchungen wurden mit Unterstützung der Deutschen Forschungsgemeinschaft im Institut für 
Strömungsmechanik der T. H. Braunschweig durchgeführt. 

\) Inzwischen ist die exakte Lösung der in [1] behandelten Aufgabe von A. N. Tifford und Sheng 
To Chu angegeben worden [4]. Die Übereinstimmung mit der Näherungsläsung [1] ist recht zufriedenstellend. 
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Es bedeutet @eine Konstante, die wir für eine Kreisscheibe mit dem Radius R nach [1]folgen- 
dermaßen berechnen können: 

wir aUWE 

© = je Ja . . . . . . . . . . . “rs . . (2). 


Hierin sind W. die Anströmungsgeschwindigkeit und V, = wR die Umfangsgeschwindigkeit 
der Scheibe. 

Zur Lösung unserer Aufgabe benutzen wir das von Th. von Kärmän [2] erstmalig 
angegebene Impulsverfahren. Da wir die uns interessierenden Formeln bereits in [1] abgeleitet 
haben, erübrigt es sich, hier im einzelnen näher darauf einzugehen. Wir erhalten danach zwei 
Impulsgleichungen, und zwar eine für die radiale Richtung und eine andere für die Umfangs- 
richtung. Wir entnehmen aus [1], Gl. (15), die Impulsgleichung für die radiale 
Richtung wie folgt: 


8 ö 6 
d dU an 
L [uud 4 | U-mar+ wart unser VB). 


Entsprechend gilt für die Impulsgleichungfür die Umfangsrichtung 
nach [1], Gl. (17), : 
ö 


d z 
— 4 p2 EEE JO ET EEE 
dr a r z . (4). 


0 


Die obere Integrationsgrenze haben wir nicht mit 2 = oo sondern mit 2 = ö angegeben, wobei 
ö die Grenzschichtdicke ist. Dieses ist zulässig, da die Integranden für z >ö verschwinden. 
Um dieses einzusehen, seien die Randbedingungen für die Geschwindigkeitskomponenten noch 
einmal angeschrieben: 
Ne, Van ero; v_ß0| (5) 
Zoo: 4 =Dl-=ar:, v0 | : 


Um mit den beiden Impulsgleichungen (3) und (4) arbeiten zu können, müssen geeignete An- 
sätze sowohl für die Geschwindigkeitsverteilungen « (r, 2) und v(r, 2) als auch für die Wand- 
schubspannungen 7,(r) und r,(r) gemacht werden, worüber im nächsten Abschnitt berichtet 
wird. 


3. Näherungsansätze für die @eschwindigkeitsverteilungen und Wandsehubspannungen 
a)Geschwindigkeitsverteilungen 
Nehmen wir an, daß wir es mit affinen Geschwindigkeitsverteilungen zu tun haben, so 
können wir allgemein schreiben 


2 


u=15(3)=1f und = rgl3)= 100) = A» Re eG) 
mit den Randbedingungen 
a ER Menge et 


2 
Wir haben gesetzt t = 5° 


Die lineare Abhängigkeit in r folgt daraus, daß einmal die vorgegebene Geschwindigkeits- 
komponente in radialer Richtung U = ar und zum anderen die vorgegebene Umfangsge- 
schwindigkeit wr linear von rabhängen. Die Gl. (6) stimmt auch mit dem Ansatz von Th.von 
Kärmän [2]für den Fall der in ruhender Flüssigkeit sich drehenden Scheibe überein. 

Die Abhängigkeit der Geschwindigkeitskomponenten von 2 bestimmen wir unter Be- 
achtung der Verhältnisse an einer längsangeströmten ebenen Platte. Wird eine Platte mit der 
Geschwindigkeit U„ längsangeströmt, dann gilt für die Geschwindigkeitsverteilung in der 
Reibungsschicht das bekannte Potenzgesetz 


m 
Eher siet « (D- 


Hierin hat der Exponent mje nach der Reynoldsschen Zahl und der Wandrauhigkeit Werte 


in der Größenordnung 1/6 bis 1/10. Als brauchbarer Wert wird gewöhnlich m = 1/7 ange- 


nommen. 
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Wenden wir das Gesetz (7) sinngemäß auf die Umfangsgeschwindigkeit an, dann muß 
vonder ander Wand bereits herrschenden Geschwindigkeitskomponente v, = wr eine durch die 
Reibung bedingte Geschwindigkeitsabnahme erfolgen. Nehmen wir die „Anströmungsge- 
schwindigkeit der Platte“ mit U. = wr an, dann wird 

De EUER bau. areys], 
wr DI 

Für die Geschwindigkeitskomponente in radialer Richtung wollen wir annehmen, daß 
das Potenzgesetz (7)in Wandnähe (z— 0) gültigist. Außerdem sind noch die Randbedingungen 
(5) zu erfüllen. Mit einem zunächst unbekannten Faktor & können wir dann schreiben 


U 2" a\2z w A 
= -() == 2321. =lweoe 1) +atll"..... (9). 
Während (8) unmittelbar mit dem von Kärmänschen Ansatz übereinstimmt, gilt dieses 
für (9), wenn darin «a = 0 gesetzt wird. 


b) Wandschubspannungen 


Bei der Aufstellung der Formeln für die Wandschubspannungen gehen wir entsprechend 
dem Vorschlag von Th. von Kärmän [2] vor, indem wir zunächst die resultierende Wand- 


schubspannung ermitteln. 
Bezeichnen wir mit u, die resultierende Geschwindigkeit in Wandnähe, dann gilt für die 


Geschwindigkeitsverteilung für kleine Wandabstände z die bekannte Formel 


us 2\” 
vo) a 


ug v 


T 
a ES A hi] 
2 y: ( ) 


die sogenannte Schubspannungsgeschwindigkeit ist. ce = c (m) ist eine von m abhängige Kon- 
stante, deren Zahlenwerte man z.B. in [3], S. 366, findet. Für m = 1/7 ist ce = 8,74. Die die 
resultierende Wandschubspannung bestimmende Geschwindigkeit ergibt sich aus den beiden 
Geschwindigkeitskomponenten zu 


ya (0 on area ee rer (12), 


wobei für v und v die Werte aus (8) und (9) einzusetzen sind. Da wir uns auf wandnahe Punkte 
beschränken (z—0), können wir schreiben 


wobei 


u Ei a Marine . - 
mals un rn Rz = N ST! 
Setzen wir in (12) ein, dann wird 
Un _ ZUR 
re (1 + aa)ıl8 (4) 2 ee 


Wir fassen jetzt (10), (11) und (14) zusammen und erhalten für die resultierende Wandschub- 
spannung 


> 1 2m 2 
a + an Ten 11 DEE (15). 
Zur Abkürzung haben wir gesetzt 
2 
-. 
k= C rm CHIC Sanet Han ET N En ee u (16); 


und zwar gilt für m = -; »k—=0,022. 


Die Komponenten der Wandschubspannung in radialer Richtung und U i 
ergeben sich in einfacher Weise zu ö 8 mfangsrichtung 
U v— wr 


ye+tT u ee 7 
Sr Dust: To EEE REINER 
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Setzen wir die Beziehungen (13), (14) und (15) ein, dann finden wir schließlich 


“4 k nr 

oe (rar a(l +02)? ee a 
ee a 

ACH (erö) (1-+. 02) Eu ee enger 3 


Der Einfachheit halber haben wir den neuen Exponenten n eingeführt, der sich aus m folgender- 
» maßen berechnet: 


Bi 1 EN 1 
Für m = 7 ergibt sich n = ver 


Wie man sofort sieht, besteht zwischen den Komponenten der Wandschubspannung die 
Beziehung 


en a ei ee re er 


4. Rechenverfahren 


a) Grundlagen 
Die Ergebnisse des Abschnittes 3 setzen wir in die Impulsgleichungen (3) und (4) des 
Abschnittes 2 ein. 


Mit 3 —=t ergeben sich nach kurzer Zwischenrechnung die Impulsgleichung 


für dieradiale Richtung zu 


1—n 


2) Iso na+2 | +207-3r +90 =anura) rar...) 


unddie Impulsgleichuug für die Umfangsrichtung zu 


ö\r wr? —_ 
Bet 2) [mar=efer ME yon 
Diese beiden Gleichungen lassen sich noch erheblich vereinfachen, wenn wir die Substitution 
2 re 
2-7) EHRE LEN, 
vornehmen. Setzen wir in (22) und (23) ein, dann erhalten wir 
1—n 
kall +02) 2 B-Um —- Dlwun,o) 2. 2m anne a (25) 
und 
en 
Bilewa) Feld) 2. ec en a (26) 


wobei wir als neue Abkürzungen eingeführt haben 


1 
1 ur 2+n 
D, (a, w, n, = lle+ af—2 P +9) a | 


Ii+n i+n 
und fr er 22) 


1 
15+3n 
D,(a, w,n,c) = |, fydt | 


Die Größen & und ß stellen die beiden Unbekannten unseres Problems dar. 
Zwischen (25) und (26) besteht noch die einfache Beziehung 


FILE ED (AO,MO) ernennen). 
11 
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i i i ücksichti Ansätze für die Geschwindig- 
Die Funktionen Ö, und Ö, lassen sich unter Berücksichtigung der ; 
Beate E ji und g nach (9) und (8) elementar auswerten. Tabelle 1 gibt eine Zu 


sammenstellung dieser Auswertung. 
Tabelle 1: Konstanten zur Berechnung der Funktionen ®, und ©, nach (27) 
: rechnung dnrFunkmene man In 


allgemein m + ın= 4 
en DIE ee Te m. 
REM ERBE JE SERIE EEE Eee rn er Bee u 
1 
A= 0,408 
1 (1+m)(2-+m) 
[fd Awa+ Ba > 
3 | B= —— 0.466 
2-+m 
1 
( = — 0,207 
(1 + m) (1 + 2m) (3 + 2m) 
1 1 
= 0,266 
C»0°?+Dawa +Ba? D man 
1 
u 0,304 
B 3+2m 
1 2 m? 
= 0,0278 
in a “SarmdaHt2m x 
3m 
= 0,0681 
1 Gneeeer 
[fg d Ga +Havw ” 
$ > 0,0292 
HS gm) Hm) 
3tn SER 
1l+n „2 Um die Abhängigkeit von m zu ‚6 
2+tn erhalten, ist zu setzen 
2 K 3,60 
IF BET. 
1+m 
ie L 


Setzen wir die gefundenen Formeln in (27) ein, dann erhalten wir 


2 
8,=(1 + BJ BR)[2) HF (DE- AN a CK 

Nie: 
D, = AU” +GLa 


Die Größe a berechnen wir aus Gl. (28), in die wir die Beziehungen (27a) einsetzen. Wir 
erhalten 
a\? 


RE EhTE) ee 


(40) 


Die Formeln und Werte für R, Sund T gehen aus Tabelle 2 hervor. 
Die Lösung der quadratischen Gleichung für a liefert dann die gesuchte erste Unbekannte 


«= x + r(2) Ze A 


Die Bedeutung der Größen X, Y und Z, sowie ihre Zahlenwerte sind in Tabelle 2 mitgeteilt. 
Da a aufgrund des Geschwindigkeitsansatzes (9) positiv sein muß, ist die Wurzel mit positivem 
Vorzeichen zu versehen. 

Die zweite Unbekannte ß, die nach (24) unmittelbar mit der Grenzschichtdicke 
zusammenhängt, ermitteln wir aus (26), wenn wir für ®, den Ausdruck (27a) einsetzen: 


i en 

B= (er = | (hol) Kl 
a 

M+Na 


Die Bedeutung der neuen Größen j, M und N und ihre Zahlenwerte gehen wieder aus Tabelle 2 
hervor. 


Se 
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Tabelle 2: Konstanten zur Bereehnung der Größen a,a*, 
der Grenzschichtdicke ö, der Wandschubspannung und 
des Drebmomentenbeiwertes 


3 1 1 
allgem = _—_» Be 
gemein m 7 N ni 
DRE= AL-2HT 
R | CK + .. 0,0294 
= F 
8 CK +GL 0,0262 
fr IB ar 
MP kan 1,0557 
x N 0,0262 
R 2 
Y (2) +T 1,0559 
R 
Z = 0,0147 
M HL 0,1342 
N GL 0,3131 
P ee ion 
[M + N(fr - zZ)" (fr - zZ) ‚0861 
BR Asa, 
9 Zitn —. un | 0,0481 


Mit der Aufstellung der Formeln zur Berechnung der Größen a und ß ist das Grundsätz- 
liche unserer Aufgabe gelöst worden. Im Abschnitt 5 wollen wir einige Fälle für verschiedene 
Werte von a/w durchrechnen. Ähnlich wie beim laminaren Fall, vgl. [1], wollen wir dabei die 
beiden Fälle > a und a > w unterscheiden. 

Der Fall A (w > a) schließt dabei den Fall a —= 0, d. h. die rotierende Scheibe in ruhender 
Flüssigkeit, mit ein. Unsere Lösung ist hierfür mit derjenigen von Th. v. Kärmän [2] 


1 
identisch, wenn man m = 77 bzw.n = Z setzt. 


Der Fall B (a > w) liefert für ® = 0 den Fall der stillstehenden angeblasenen Scheibe 
(rotationssymmetrische Staupunktströmung). Während wir für den Fall A unmittelbar die 
oben abgeleiteten Formeln (29a) und (30) anwenden können, empfiehlt es sich für den Fall B 
noch eine Umformung vorzunehmen. Mit 


FR Ms „(2 a a en ° 
+ =—a und ß E ß oe (Sl) 
erhalten wir 
5 W\? 
= (+ r-z BREIT TR a 


und 


gt 18 (e SE (83). 


v ee, BI No. 


Für die Staupunktströmung (w = 0) erhalten wir außer aus (33) noch eine zweite einfachere Lösung, wenn wir nur 
mit der Impulsgleichung für die radiale Richtung rechnen. Diese Lösung ergibt sich mit v = 0 aus (3) zu 
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REN RE BE RE SE SE ae I 


Für die Geschwindigkeitsverteilung machen wir den Ansatz 


s-lö)-r 


Für die Wandschubspannung setzen wir analog den Überlegungen des Abschnittes 3b 
Tr k 


 — —— 


gar)? (aroin 
v 
Setzen wir uund 7,in die Impulsgleichung ein, dann erhalten wir entsprechend (22) 


1 1 
k 
„ [ra-mar, [arm -3ena= gm 
0 


dr arö\n 
v 


Führen wir die Integrationen aus, dann erhalten wir nach kurzer Zwischenrechnung die Lösung für die Grenz- 
schichtdicke der Staupunktströmung zu 


1 
ne ern 
Bi ( ) los.) 02 se (35). 


Nach dieser Zwischenbemerkung wenden wir uns jetzt der Ermittlung der anderen Grenz- 
schichtgrößen zu. 


b) WandschubspannungundDrehmoment 


Die Wandschubspannungen in beiden Richtungen r, und 7, berechnen wir 
nach (18) und (19). Setzen wir (30) ein, dann wird 


1—n n n 
7. ; 2(1+n) @ itn/wr?\ 14m 
Te Wa re a} 
1—n n n 
T, 2 2 (14) a ı4n/wr?\ i+n 


Unter der Annahme, daß die Reibungsschichtdicke klein im Vergleich zum Scheiben- 
radius R ist, können wir das Drehmoment einer einseitig benetzten Scheibe wie folgt 
berechnen: 


R 
Ms—ln jr ae Be 
) 


Die Abhängigkeit der Wandschubspannung 7, vom Radius r finden wir aus (37) zu 


2 
oe) _ an 
Te (R) riW 
Führen wir den dimensionslosen Momentenbeiwert 
M 
Cm a —, . . . . . . . . . . . . . . . (40) 


4 
> 08 Bö 


u 


ne), 


ein, dann ergibt sich nach Ausführen der Integration über r 
i+n 7(R) 
Ran el ee 


Zur Auswertung dieser Formel ist die Wandschubspannung nach (37) für r = R einzusetzen. 
Für große Werte von a/w nimmt der Ausdruck (41) eine besonders einfache Form an. 


Beachten wir, daß hierfür die Größe « = (Yr— Z)>1 wird, dann können wir schreiben 


aM=—4 


[00) 


1 n 
T(R) _ , ae a 
o(w R)? | v 

wobei die Formel und der Zahlenwert für Pin Tabelle 2 angegeben worden sind. 


(aa, 
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Wir setzen in (41) ein und erhalten für den Momentenbeiwert bei großen Werten von 


a/w den Ausdruck 
1 
| na 
I+n © 


na (om (vw <a ) RE (43). 


v 


Auf eine Besprechung der in diesem Abschnitt abgeleiteten Formeln wird im folgenden Abschnitt 
=} N, 1 1 
eingegangen, in dem wir einige Beispiele durchgerechnet haben. Für m = — bzw.n= — 


7 4 
haben wirin Tabelle 3 sämtliche Endformeln noch einmal zusammengestellt. 


Tabelle 3. Zusammenstellung der Formeln zur Berechnung der 


Grenzschichtgrößen für n-- 


Talli Ga 
& 
& V0,0262 + 1,0559 p® — 0,0147 p 
or 8 (1 +a90,3 
ver 0,0481 5.7342 9 + 0,31310)0,8 
5 2 
or Tr 
Van Slam 0,0481 & (1 + 02)0,3 (0,1342 p + 0,3131 0)0,2 
2 
or To 
re — 0,0481 (1 + 02)0,3 (0,1342 p + 0,3131 0)0,2 
E 2,1518 |/2° 
Kr ae ver) 


o* V0,0262 g? + 1,0559 — 0,0147 
5/5 
Eee (q? + 0**)0,3 
V Ha 0,0481 (5,1342 4 0,3131 +08 
VL 0,0481 &* (9? + &*2)0,3 (0,1342 + 0,3131 a*)0,2 
v o(ar)? ö EG rs ? San 
jr w_ 0,0481 q (9? + a*2)0,3 (0,1342 1 0,3131 &*)0,2 
Year N! se I 


5 /,. nD2 5 2 
& R? & ie ar —% 
— a. — 2,7318 ze . » on: 


5. Ergebnisse 


| ment Ze ee a ER 

Mit den Formeln der Tabelle 3 (n = 4) haben wir für die Werte p = 0; 25% 

und 1, sowie qg = I 0, ; E : : | und 1 die Grenzschichtgrößen ermittelt. Die Werte haben 
a 


wirin Tabelle 4angegebenundinden Abbildungen 2 bis 6graphisch dargestellt. 


a) Grenzschichtdicke (Bild 2) 


Hält man die Winkelgeschwindigkeit fest, dann nimmt die Grenzschichtdicke (obere 
Kurve) sehr schnell mit wachsender Anströmungsgeschwindigkeit ab. Der Ordinatenwert für 
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a/o = 0 (rotierende Scheibe in ruhender Flüssigkeit) ist zu 0,526 berechnet worden, während 
Th. von Kärmän [2] 0,462 angibt. Der letzte-Wert beruht aber auf einem Rechenfehler, 
wie man leicht zeigen kann. Der richtige Wert 
ist übrigens auch von H. Schlichting[3] 
schon angegeben worden. 

Hält man dagegen die Anströmungs- 
geschwindigkeit konstant, dann ändert sich 
die Grenzschichtdicke (untere Kurve) nur 
wenig mit . Der Ordinatenwert für &/a = 0 
(rotationssymmetrische Staupunktströmung) 
ist zu 0,092 berechnet worden, während man 
für die Lösung nach (35) 0,086 findet. Die 
Übereinstimmung ist. recht zufriedenstellend. 


b) Wandschubspannung 
(Bild3 und 4) 


Bei festgehaltener Winkelgeschwindig- 
keit nehmen beide Schubspannungskompo- 
nenten (untere Kurve in Bild 3 und obere 
Kurve in Bild 4) mit wachsender Anströmungs- 
geschwindigkeit beträchtlich zu. Ebenfalls 
wächst bei konstanter Anströmungsgeschwin- 
digkeit die Schubspannungskomponente in 
Umfangsrichtung (untere‘ Kurve in Bild 4) 
sehr schnell mit wachsender Drehgeschwindig- 
keit, während sich die Schubspannungskom- 

Bild 2. Grenzschichtdicke ö; e = Wert entsprechend (85) ponente in radialer Richtung (obere Kurve in 
Bild 3) nur wenig mit der Drehzahl ändert. 

Der Ordinatenwert für w/a — 0 (rotationssymmetrische Staupunktströmung) beträgt 0,0419, 
während sich nach der Lösung nach (35) 0,0414 ergibt. Die Übereinstimmung beider Wert 


T; 
9 gan?“ 


- 004 


-002 2 
je — a 


007 


0,2 0,4 06 08 10 


a,@ 2 

w’a a 

Bild 3. Wandschubspannung in radialer Richtung T, 
° = Wert entsprechend (35) 


;@ 
a 


Bild 4. Wandschubspannung in Umfangsrichtung Top 


ist sehr gut. Wir haben damit gezeigt, daß wir mit unserem Verfahren, den gesamten Bereich 


2 g l e b = 


2. . R e “= ö - 3 A 
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e), Drebimoment (Bild > und 6) 


5 2 


In Bild 5 haben wir über dem Parameter — den Drehmomentenbeiwert Cy 


der einseitig benetzten Scheibe aufgetragen. Man sieht, daß bei konstant gehaltener Winkel- 
geschwindigkeit das Drehmoment mit wach- 
sender Anströmungsgeschwindigkeit sehr 
schnell zunimmt. Für große Werte von 
a/w (»< a) errechnet sich die Änderung aus 
(43) zu 


0,113 en 
cu = Tor en (turbulent) (44). 
1 
Im laminaren Fall (vgl. [1] Gl.(49a)) hatten 
wir eine schwächere Änderung mit a/o gefun- 
den, nämlich 


> (laminar) . . (45). 


Schließen wir auch den Bereich kleiner Werte Birds. Drehmomentenbeiwert cyp — M/ - a ( 2) 
von a/win unsere Betrachtungen mit ein, dann = 
können wir in guter Näherung annehmen, daß sich das Drehmoment linear mit der Anströ- 
mungsgeschwindigkeit ändert, 


700c 
70 


1 1 7 Bil ea a 
HH HH 


OS 
[eo } 
I If 1 eier 


0* 2 De 10° 2 wer 2 ee 


oO 
Bild 6. Drehmomentenbeiwert der einseitig benetzten rotierenden Scheibe cy >= M ie w® R® 


Für die in ruhender Flüssigkeit rotierende Scheibe gilt in Übereinstimmung mit von 
Kärmaän [2] 


re ee (Fuck lenb) Me Bann Den 350] (AG)! 
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———— 
Tabelle 4: Zusammenstellungder Erge —— 

ee on 8 |tjar_ u | far TE 
a | Ag: | = vr | v 0 (wr)? E v ee 
2 ’ 
= 0 oo 0,1618 0,526 0,0043 0,0266 0,073 
ir 1/4 4 0,2999 0,256 0,0098 0,0326 0,089 
1/2 2 0,5314 0,166 0,0206 0,0388 0,106 
3/4 4/3 0,7765 0,130 0,0347 0,0446 0,122 
1 1 1,0255 0,112 0,0524 0,0511 0,140 
i “ ar ee re | SSR: 
TE a* ae Teen. er 2 TI M 
@- a er v o(ar) v en! v 
Sehe De en : > 
= 1 it 1,0255 0,112 0,0524 0,0511 0,140 
= 4/3 3/4 1,0199 0,104 0,0480 0,0351 0,162 
1,0161 0,098 0,0222 0,211 
1,0136 0,093 0,0105 0,348 
1,0129 0,092 ©?) 


zZ 
"2 
Bild 7. Geschwindigkeitsverteilungen Fall A, o>a O=2z=Ö Bild 8. Geschwindigkeitsverteilungen Fall B sa> wow O=2=Öö 
während sich im laminaren Fall 1.74 
9 . 
Cm = (laminar) BEN a Ka ee ee 47 
v 


nach unserer Näherung, vgl. [1], ergibt 3). 


2) Asymptote: 0,113 (a/w)9,8 
3) Der exakte WertnachCochran ist 1,935 statt 1,74, 
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Um die Ergebnisse des Bildes 5 für das Drehmoment in bequemer Weise auf die ange- 
strömte Scheibe vom Radius R anwenden zu können, wollen wir noch die Größe a/w gemäß (2) 
durch das Verhältnis aus der Anströmungsgeschwindigkeit W, zur Umfangsgeschwindigkeit 
V,= w Rausdrücken. Damit wird aus (44) für die asymptotische Lösung 


0,8 ’ 
le) a TR ee 


RW 
5 


In Bild 6 haben wir den Momentenbeiwert über der Reynoldsschen Zahl w R2/» mit 
W./V, als Parameter aufgetragen. Miteingezeichnet sind zum Vergleich die entsprechenden 


Kurven des laminaren Falles nach [1]. 
Die aufgetragenen Kurven gelten für die einseitig benetzte Scheibe. Wenn wir näherungs- 


weise annehmen, daß hinter der Scheibe infolge des Totwassers keine Geschwindigkeitskompo- 


—— radiale Richtung —— radiale Richtung 
| rl, 


Umfangsrichtung ——— 
Umfangsrichtung —— 


Bild 9. Polardiagramme der Geschwindigkeitsverteilungen; t = 2/ö 
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nente in Achsrichtung herrscht, dann können wir den Beiwert der beidseitig benetzten Scheibe 
(bezogen auf die einfache Scheibenfläche # = R?) etwa wie folgt berechnen 


ya = Cmı + emo - + (49). 


Es bedeuten cyr, den Beiwert der in ruhender Flüssigkeit rotierenden Scheibe, Cy ‚den Beiwert 
der einseitig benetzten Scheibe (Bild 5 und 6) und cyr, den Beiwert der beidseitig benetzten 
Scheibe. 

d) Geschwindigkeitsverteilung (Bild 7, 8 und 9) E 

Die Geschwindigkeitsverteilungen in radialer Richtung und Umfangsrichtung sind durc 

(9) und (8) gegeben. In Bild 7 haben wir die Verteilung für den Fall A (w>.«) und in Bild 8 
für den FallB (a>w) dargestellt. In allen Fällen zeigt sich die schnelle Anderung der Ge- 
schwindigkeitskomponenten mit wachsendem Wandabstand, die bei turbulenten Strömungen 
infolge der Mischbewegung entsteht. Eine anschauliche Vorstellung von der Geschwindig- 
keitsverteilung gibt auch das Polardiagramm in Bild 9. Hier ist ähnlich wie für den laminaren 
Fall in [1] von einem Pol aus der wandparallele Geschwindigkeitsvektor mit den Komponenten 
w und v für äquidistante Wandabstände t = z/ö dargestellt worden. 


6. Zusammenfassung 


Es wird die turbulente Strömungin der Nähe einer unendlich ausgedehnten rotierenden 
Scheibe berechnet, die gleichzeitig in Richtung der Drehachse angeströmt wird. Die Lösung 
wird mittels des Impulsverfahrens näherungsweise gefunden. An der Scheibe bildet sich eine 
Grenzschichtströmung mit Geschwindigkeitskomponenten in radialer Richtung und Umfangs- 
richtung aus. Die Dicke der durch die Drehbewegung mitgerissenen Schicht (Grenzschicht), 
die Wandschubspannung und das Drehmoment, sowie auch die Geschwindigkeitsverteilung 
in der Reibungsschicht sind stark abhängig vom Verhältnis der Anströmungsgeschwindigkeit 
zur Umfangsgeschwindigkeit der Scheibe. Während der Drehmomentenbeiwert im laminaren 
Fall umgekehrt proportional der Wurzel aus der Reynoldsschen Zahl w R?/» ist, ändert 
sich der Drehmomentenbeiwert im turbulenten Fall umgekehrt proportional der fünften 
Wurzelausder Reynoldsschen Zahl. Die Abhängigkeit vom Verhältnis der Anströmungs- 
geschwindigkeit zur Umfangsgeschwindigkeit ist im laminaren Fall proportional der Wurzel 
aus diesem Verhältnis, während sie sich im turbulenten Fall nahezu linear hiermit ändert. 
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Über die Verteilungsfunktionen von Vektorsummen 
Von E. A. Deuker in Vernon (Eure) 


In einem p-dimensionalen, euklidischen Raum seien n Vektoren von beliebiger zunächst fester Länge 
gegeben, deren Orientierungen statistisch mit gleicher Wahrscheinlichkeit verteilt sind. Der Betrag des 
Summenvektors ist eine statistische Größe, für deren Verteilung eine allgemeine, für alle positiven, ganz- 
zahligen p und n gültige Formel in Gestalt eines uneigentlichen Integrals über Besselsche Funktionen her- 
geleitet wird. "Diese Verteilungen gehen, unter bestimmten Bedingungen, für große Werte von n asymptotisch 
in bekannte Verteilungen der mathematischen Statistik über. 


The autor considers n vectors of any length in a p-dimensional Euclidean space. The directions of the 
vectors are statistically distributed with equal probability, so that the length of ihe sum vector is a random 
variable. The distribution function is stated, for all positive integer values of p and n, in the form of an 
improper integral of Bessel functions. Under certain conditions, the distributions tend, for large values ofn 
to well-known distribulions of mathematical statistics. ; 


Dans un espace Euclidien p-dimensionnel n vecteurs d’une longueur quelconque sont donnes. Leurs 
directions sont distribuees statistigquement avec probabilite &gale. Le module de la somme est une variable aleatoire 
pour la distribution de laquelle, on donne une formule generale, valable pour toutes les valeurs positives, entieres 
de pet n sousla forme d’une integrale infinie de fonctions de Bessel. Sous certaines conditions, ces distri- 


butions, pour de grandes valeurs de n, se changent asymptotiguement en distributions connues de la statistique 
mathematique. 


| 
j 
h 
4 
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llyctp B p-MepHoM EBKAUNOBOM IIPOCTPAHCTBe NAHO n BEKTOPOB JI000H, BHAYA,Ie, OIPe- 
AelleHHOU_ AIMHEI, OPMEHTAIHN KOTOPEIX CTATUCTHyeckH pacıpeneleHkl C PaABHOH BepoAT- 
HOocTbI. Monysp CyMMHpyIomeroO BEKTOPA ABJIAETCH CTATUCTUYECKOH BEIHYHHON, IA Pacıpe- 
AeneHnA KOTOPOH BBIBOAHTCA O6bImaa dOpMyıa, MeÄicTBuTeNBHaA MIA BCeX IIOJIOKHTEJIBHRIX, 
HEeJIOYHCJHIEHHBIX P U n, B BU] HECOÖCTBEHHOTO HUHTELPana mo PyHurunam Beccena. Ilpn onpe- 
MeJIeHHBIX YCIOBHAX TH PACIpeneteHHA aCHMITOTHYEeCKH MepexonAT AA 6OJIBIIHX BHAYeHMH n 
B XOPOMIO HUSBeCTHBIeE P&CHPpefelleHHA MATEMATUYeCKONM CTATUCTUKH. 


1. Einige einfache Beispiele für p = 2 


Wir gehen aus von dem Sonderfall 9 = 2, dessen vollständige Erledigung wir an den 
Anfang stellen, und wählen als Beispiele n = 2, 3 u. 4. Die zwei, bzw. drei u. vier Vektoren 
sollen außerdem gleiche Länge haben, die wir gleich 1 setzen können. Eine Vorstellung von den 
zu betrachtenden Wahrscheinlichkeiten bekommt man etwa, wenn man zwei bzw. drei oder 
vier Streichhölzer viele Male auf eine Tischplatte wirft, sie, entsprechend der durch den Kopf 
gekennzeichneten Richtungen, aneinanderfügt und das Verhältnis der Anzahl mit der Längen 
zwischen —1/2 di +1 und 1 + 1/2 dl aufgetreten sind zu der Gesamtzahl der Würfe über die 
Länge I aufträgt. 

Da alle Richtungen der Vektoren gleiche Wahrscheinlichkeiten haben sollen, sind die 
entstehenden Konfigurationen für jedes n nur bis auf eine Drehung in der Ebene bestimmt, 
so daß man einen beliebigen Vektor als Bezugsrichtung herausgreifen kann und sich auf die 
relative Lage der anderen Vektoren zu ihm beschränken kann. Im Falle n = 2 verbleibt dann 
nur noch ein Freiheitsgrad (s. Bild 1) und die Wahrscheinlichkeit, daß der zweite Vektor 
innerhalb des Winkelraums dg fällt, ist: 


Wun=,. — nr<gq<sae bzw. AR. — I dy Degen. 
Mit 
I=2.c0s Zr ie re a TE N) 
folgt: 
: 1 dl 
du.=— = “ 0 0. Aue) (1.3) 
PR 
also 
Bar 0<I=<2 
Rum. 2 we 
| 1 ı>2 


mit der Frequenzfunktion 


Bild L 
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So einfach die Herleitung und das Resultat in diesem Falle ist, so ist es doch der Ausgangspunkt 
für die Ableitung der allgemeinen Formel, wie wir später sehen werden. a 

Beim zweiten Beispiel n = 3 ist die direkte Berechnung der Verteilungsfunktion nicht 
mehr trivial, sondern bereits einigermaßen mühsam. Man hat: 


Fu | | Qua. TEE RR 
mit i 
B=3 +2[csg Feosg +cos(1 —9)) - - +: (EA 
Durch die Substitution: 
ı=Uu—v; Gem u... ei ae Bene (6 


gehen (1.6) u. (1.7) über in: 
Fu-za)| du dv ; 2— 1 +4cosv(cosW-+cose) - - - - (1.9) 


und nach Ausführung einer Integration folgt: 
( I+1 
alc C0oS | — u 


Pe, 
= are cos|1 — — cos u] dv NE! 
TC An: 4 cos v 
Re ae . 5:(1:10% 
alt COS —— 
E 3 
Ne | arecos|. 2 eos o| dv 1.2155 
| eo; 4cos v 


Einen Überblick über die jeweiligen Integrationsbereiche ] = const gibt Bild 2. 

Wir übergehen hier die Einzelheiten der weiteren Rechnung, die auf vollständige ellip- 
tische Integrale erster Gattung führt und schreiben nur das Ergebnis an. Man erhält als 
Frequenzfunktion: 


5. 
VER. K) 0<I1<1 


NE. re (A, 
| 7 
(0 >38 


wobei Ä(k) das vollständige elliptische Integral erster Gattung ist, und der Modul k als Funk- 
tion von l, sich wie folgt, berechnet: 


161 
; 8 zır Stel 
"OT |e-na+n era rel) 
Speziell gilt noch: 
1 
F.l)=Z: 
In ähnlicher Weise kann man im Falle n = 4 bei der Aufgabe: 
1 
Fa — u dgı ne dgs “ dgs m Kalle 2a a s . . . (1.13) 


?—=4+2[cos q, + cos ga + c05 95, + cos (9, — 95) + cos (g3 — 95) + cos (9, — 91)] 
vorgehen. Eine längere Rechnung führt zu dem Ergebnis: 
1 


f k(a)- Kik(D)] de = Klik) Pe 
0 Yı+21s+8 1—a V a 


1+2/2 +2 y1— a2 


2<sıs4 


} 
) 
a j Klik) da 
Rz 4 Yı+21c+% 1— 22 
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Es ist wieder K(k) das vollständige elliptische Integral erster Gattung und mit der Abkürzung 
y=!+21ogilt: 
N re a —1<ys0 
8+4y— y?+8fl+y 
syl +Yr+8+4y— 7? 
16y1 +y 


Wie man sieht, erhält man bereits für kleine n, recht verwickelte, aus höheren transzen- 
denten Funktionen aufgebaute Ausdrücke für die Frequenzfunktionen, im Gegensatz zum 
Fall p = 3, wie wir sehen werden, und mit n — 4 dürfte praktisch die Grenze für eine direkte 
Behandlung der Aufgabe: 


ee 
ns: [ie dga* +» dgn—ı | ER (1.16) 
?=n-+2[cos g, +C0s gg +++ COS n-ı + C08 (9ıı — 92) +: --] 


gezogen sein. 


ke — RT 


'0<ys8 


U. Der allgemeine Fall p = 2 


Unser erstes Beispiel gibt aber die Möglichkeit eine Rekursionsformel aufzustellen, die 
sich allgemein lösen läßt. Wir verallgemeinern zu diesem Zweck das erste Beispiel auf den Fall 
beliebiger Länge der beiden Vektoren, die # und a„+1 seien. Die Verteilungsfunktion ist 
dann, wie man leicht nachrechnet: 


[0 <a — 1] 


21 1 

ns — | — mul<I<u ra 241%. 
n Tara latmmn! Hi n+1 (2.1) 
0 I> at 

Es sei nun x die Länge des Summenvektors von n Vektoren mit der Frequenzfunktion fg,» (%), 


dann ist } die Länge des Summenvektors von n + 1 Vektoren und nach den Gesetzen der 
Wahrscheinlichkeitsrechnung gilt: 


Ja2,.(®%; Int» D) = 


ne d 
x) d® 
Sant) = m | V4 a? en. an— pp TR 
=@n+1l 
Macht man in (2.2) die Substitution: 
2 — 7.4 any = 2 anlı > cosO ne 2), 
so werden die Grenzen des Integrals 0 und x und man erhält: 
l 1 R 
aD 1 [eiiii er 
mit 
u ET 


Es gilt nun, eine Funktion zu finden, die ein’auf R bezügliches Additionstheorem besitzt und 
das ist die Be sse1sche Funktion Nullter Ordnung J, (R). In der Tat erfüllt die Funktion: 


on 

Al LH Jlat) SI LTHdeN. rr. .(2,6) 
0 k= 

die Rekursionsformel (2.4) und stellt damit die allgemeine Formel für die Frequenzfunktionen 

von Summen ebener Vektoren dar, wenn wir noch zeigen, daß sie, z.B. fürn =2, Gl. (2.1) 

liefert. Aus dem Additionstheorem der Besselschen Funktion J, (R) folgt nämlich: 


E13 x oo 


l n n+1 
| Mia IR Erd Aldıd..... an 
TU k=1 st 


(al) ö 
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ey 
und für n = 2 gilt speziell: 
Aula) Aldi... (2.8). 
6 


Dieses Integral ist aber nach der Formel von Sonine u. Dougall!) gleich: 
[0 Iı< la — @2| 


— au I — Dr En ar 

e& <i<a UM 6 2 ) 

J..(' z Vaatai (a? ai 1)3 | 1 2| 1 2 
& l>a, +0 


also identisch mit (2.1). 


Mi 


halk) ayraz:! 


Bild 1 


bulk) a 


EEE 
LE 


hslk) 


34792795°9,:9,=1 


3 


03 


02 


U 


hs) 


03 
ERERRERERERT] 


02 


01 


| 
In 


Bild 5 


Aus (2.6) folgt durch Integration, daß die Wahrscheinlichkeit, daß eine Länge zwischen 
O0 und J auftritt, der Ausdruck 


Fal=L] Hai) Ildd 0.2.2.0. 20... 010) 
- 0 k=1 


2)s.2..B.Magnus-OÖberhettinger, Formeln und Sätze für die speziellen Funktionen der matbe- 
matischen Physik. Berlin/Göttingen/Heidelberg 1949. 
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ist. Damit ist zunächst für den Fall 9 = 2eine allgemein gültige Darstellung für die Verteilung 
der Summe von n-Vektoren gegeben. Einige Beispiele spezieller Frequenzfunktionen f,, n(}) 
geben die Bilder 1—5. Die gestrichelten Kurven beziehen sich auf die entsprechenden 
asymptotischen Ausdrücke (4.20). 


3. Der Fall beliebiger Freiheitsgrade p 
Wir betrachten wieder die Summe von zwei Vektoren mit den Beträgen # und 4.1. 
Der Endpunkt der Strecke a„+1 liegt jetzt, statt auf der Peripherie des Kreises vom Radius 
@n+1, auf der Oberfläche einer p-dimensionalen Kugel vom Radius a„+1ı und für dF,,: (l) 
können wir schreiben: 


dO, 

STE DE en 
wenn O, die Oberfläche der p-dimensionalen Kugel bedeutet. Für 1? gilt die analoge Beziehung 
wie (1.2.) 

Pi a at + 28 000g ae Era en 


Unter Einführung p-dimensionaler Kugelkoordinaten findet man: 


re 


en N re N 
dF,,2 ET sin? g:dgq 
2 
KiEs E 
2, I[4 +18 — (m+ı + © — P)?] 2 A (3.3) 
also als Frequenzfunktion 
0 1> | — an+ıl 


je ö 
Bir = IA a0 — (a + BT 
a: I(p—1) Bier 

0 I> a + An+ı 


Faßt man wieder x als Länge des Summenvektors von n Vektoren mit der Frequenzfunktion 
Fo, (x) auf, so erhält man die zu Gl. (2.4) entsprechende Rekursionsformel, wenn man ebenfalls 
die Substitution (2.3) macht: 


re ee 


Fo. (l) = 


|® = +ıl He + An+i (3.4). 


lan lest, 
= ee sin Da BENDER 2310) 
ö 
Diese Rekursionsformel wird durch den folgenden Ausdruck erfüllt: 
2 D an 22 . n JP_, (a; 1) 
I» n(l) = 2 2 he r(3) . 12 | 22 er Ile Jp_, (l t) ıpl2 dt aan (3.6); 
\ 27. er N 2 
(a, 1)2 


denn nach dem Additionstheorem der Besselschen Funktionen ist: 


Jr» _,(Ri) » —, =; FED ai (An+ı ©) 
REDE 7-1): 2(&-14m) : Er 
Bei 2 m=0 2 z —i 
: #1 
J» 1rm(lt) ER 
x ; Fey (CDT Ay fa Bash se Fuel ale), 
Pe 


wobei die Funktionen C„ (cos 9) die Gegenbauerschen Polynome bedeuten. Da diese 
Polynome zueinander orthogonal sind entsprechend der Formel: 


n 0 m Zz m 
On (cos 9) - Om (cos 9) sin?’ ddd — - I'(2v + m) 
0 22 —1(y + m) m! [I (v)] 


Mm (3,08) 


a e B ° F h. 
168 . Deuker, Über die Verteilungsfunktionen von Vektorsummen p4.34 ar Marıası 


gilt speziell: 
0 mu 


| 62° (cos 9) sin? Bd Fremen a: 


so daß man für die rechte Seite der Gl. (3.5) mit dem Ansatz (3.6) erhält: 


Pf) ra-n 


2 
2 
aa N en ae) 
“ 2 pe 
u [92 (&) [> 1 r(2)) uf) (a4) 
x [ u — | Ja _ (Ai d= — ——— | I | — Ip _,(lt)tPi? dt 
k=1 2 Er k=1 SALE 2 
0 (art) 2 


(a, {2 ri 12 N) 


(3.10) 
also, wie behauptet wurde, usa ‚da Gl. (3.6) für n = 2 Gl. (3.4) liefert. 


Unter Einführung der ganzen Funktion: 
1 v 
A (» a 2) — (2) "Io MDF IE 3 en (3.11) 


schreibt sich noch Gl. (3.6) bequemer: 
Ip,» ee ar Dr (® 2 ak e) -Js_ (HPRd—.. . (8.12 
BEur (2) zeit ara 2" 
Die der Gl. (2.10) entsprechende Formel lautet: 
F,,n (l) en Ze p a ?_ı 
93-1 r(2) un e er at) Joe (ii)te "dr. . 23.15). 


-1 


Mit den Gl. (3.12) und (3.13) haben wir die einheitliche, für alle positiven ganzzahligen p und n 
gültige, Darstellung für die Verteilung des Betrages einer Vektorsumme erhalten. 

Mit diesen Formeln beherrscht man auch weitgehend die Verteilungen von Summen, 
wenn die einzelnen ‚„Summanden‘ selbst irgendwelchen Verteilungen mit den Frequenz- 
funktionen h,, (@,) gehorchen. Wir werden später dafür Beispiele geben, wollen zunächst aber 
einige allgemeine Sätze, einige Sonderfälle, sowie das asymptotische Verhalten der Primitiv- 
Verteilungen (feste a,) für große n untersuchen. 


IV. Einige allgemeine Sätze, Sonderfälle und asymptotisches Verhalten für große n 


n 
Da f»,n(l) für I> 2% offenbar identisch null ist, bestätigt sich z. B. auf diesem Wege 


eine in Formelsammlungen zu findende Beziehung zwischen Besselfunktionen, wie sie 
etwa bei Magnus-Oberhettinger in der Form: 


u - di 
| Ira) = 0 a4,>4, 4 Q; +. +0, ee (4.1) 
0 


angegeben ist. Weiter erscheint nach Gl. (3.12) f,n (l) as Hankel- Transformierte, und 


die Umkehrung der Hankel- Transformation liefert die Beziehung: 


oo 


p In: 1 
II FE, En a) = (rn: rl? zer) ai RE: 
‚ 0 


k=1 
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Die G1.(3.12) und (4.2) sind als die Erweiterung der in der mathematischen Statistikfür p = 1 
benutzten Fourier- Transformation anzusprechen, wobei die Fourier- Transformierte 
der Frequenzfunktion gewöhnlich als charakteristische Funktion bezeichnet wird. Im Sonder- 
fall p = 1 lauten nämlich die entsprechenden Formeln: 


oo 


| IF cos agt- cos It dt a Te eh 
7) kei 


sin bt 


N 
Pu) =&| Ifcos aut: N 
g = 


Den Streckenspektren im Falle p>1,n>1 steht also im Falle p =1 ein diskretes Spektrum 
gegenüber. Es ist angebracht, hier die Beschränkung auf positive und a, aufzugeben. Es ist 
dann zu schreiben: 


=] Te oe a SE NEE 
0 


k= 


- 


RE SEE 


Die charakteristische Funktion für die Summe von rn positiven oder negativen Zahlen a; ist 
also das Produkt: 


N 
rl) lheesapt.. 2. er 
k=1 
Im allgemeinen Fall sind die charakteristischen Funktionen die Produkte: 


. p 
mn) Hor, Re 1 


und an die Stelle der Fourier-TIransformation tritt die Hanke 1-Transformation. 


Alle Momente gerader Ordnung sind aus Gl. (4.2) durch Vergleich gleicher Potenzen von t 
abzulesen. Insbesondere gilt für das zweite Moment die Beziehung: 


ep) = at NETTER LT: 


Die Momente ungerader Ordnung dagegen, insbesondere also auch das arithmetische Mittel: 


ee a In Er (4.8), 
0 


sind nicht so ohne Weiteres anzugeben, bis auf die durch die Gl. (4.3a) und (4.4a) definierte 
Verteilung für p = 1, bei der alle ungeraden Momente verschwinden. 


Die allgemeine Formel für die Frequenzfunktionen weist ferner auf einen prinzipiellen 
Unterschied zwischen geraden und ungeraden Dimensionen hin. Für ungerades p hat man es 
mit Besselschen Funktionen halbzahliger Ordnung, also elementartranszendenten Funk- 
tionen zu tun, während für gerades p Besselsche Funktionen ganzzahliger Ordnung ein- 
gehen. So erhält man z.B. fürp =3: 


oo r 2 f 
fs.n() = r | fg ee 


SIDE DO EEE rer. (4,9% 
gie 0% 


Im Gegensatz zu p = 2sind alle Frequenzfunktionen rationale Funktionen, nämlich Polynome 
in I, von denen wir zum Vergleich mit p = 2einige für Vektoren gleicher Länge 1 angeben. So 
gilt: 

12 
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ki 
— 0 
Js. () .. 2 
0 we 
Br Ve 
Js,s (l zn 
an i2122 
Bi 
8-3) 0<I=2 
Js,. (I 17 
116° Kae ah 
1 . (4.10). 
re 1. 
mo 2° 15P 230: u u Hei 
e en 3<I<5 
Er 2, (06 - MP +57) ed 
RE 
fd = | g.95 (RL HAM —5R—-M0)  Aeicd 


Br (67 | 1<1<6 


Das asymptotische Verhalten der Verteilungen für großes n kann aus der Darstellung (3.12) 
unter gewissen Voraussetzungen leicht erschlossen werden. Zunächst stellen wir fest, daß der 
Betrag des Faktors in der charakteristischen Funktion (4.6) für 2>0 kleiner als 1 ist. Die 
charakteristische Funktion strebt also für großes n für alle#> 0 gegen Null, sofern die a; nicht 
selbst gegen Null streben. 

Wir schreiben (3.12) nun in der Form: 


& » 1 : 
ae Im n(3:- Zee) 


Jon(l) = j BE 35(kt) erlitt. nA: 


> : 2 1 
En = 2 »die Nullstellen von „1, (2; FT 2): 


Es seien &» 
2 


Dann gilt für alle z die Produktdarstellung: 


1 co 22 
A ee = ae N ech 


und man hat: 


Ei = “ at? 
Ipl2 me Br 2 ei 
Jrn(l) = gm mid, (di... (4.13), 
2 


Entwickelt man jeden Summanden in eine Potenzreihe, so ist der Koeffizient von {?? gegeben 
durch: 


il; n co an’ 
oh ee . Su eg . . . . . . . D . . . . . 
a: a 
=—1 m 
Setzt man noch zur Abkürzung: 
N 
Ber = Urn en Rn Re 


Z. angew. Math. Mech. a P . E 
BESUM, 4/5 April[/Mai1gpa Deuker, Über die Verteilungsfunktionen von Vektorsummen 171 
ee en ee N RE EEE NET ES Ne EN 


so hat man für die charakteristische Funktion 9,, (t) die Re; 


N 1 2 A 
= are ler teren Lie) 
k=1 2 4 
Wir machen nun die Substitution: 
a ee nn 


und erhalten: 


a 7=,6) va)" [m GEEZER Ver) Eden. (kis) 


Unter der Voraussetzung, daß die Bedingung 


imd-0.... HN EN EIg) 


erfüllt ist, wie es z. B. für den Fall gleicher Strecken der Fall ist, gilt dann: 


oo 


p/2 m\? le er 
Jnl)7 : Var. an ( Bir)an2ar 
22-1r(2) N en 2 
2 . 


Ar ee 


Diese asymptotische Verteilung ist aber identisch mit der in der theoretischen Statistik als 
4°-Verteilung bezeichneten, und entspricht z.B. für p=3 der Boltzmannschen Ge- 
schwindigkeitsverteilung in einem idealen Gas. Für 9=1 erhält man speziell die halbe 
Gau Bß-Kurve, oder wenn man von Gl. (4.3a) ausgeht, die vollständige Normalverteilung. 
Ist die Bedingung (4.19) nicht erfüllt, so werden die Verteilungen im Allgemeinen auch nicht 
das betreffende asymptotische Verhalten zeigen, wofür wir ein Beispiel geben wollen. Wir 
gehen aus von Gl. (4.3a), betrachten also die Frequenzfunktion 


© 


fin -- Tronkane cos. Edi Zee 


U gakeı 


und wählen für die Zahlen a, die geometrische Folge: 


. (4.22) 


n 
Die charakteristische Funktion, das Produkt //cos a,t strebt jetzt für n— ooeinem end- 


lichen Grenzwert zu.-Es gilt nämlich: k-1 
o t sin t 
au er (2) Se a en (4.23) 


und die Frequenzfunktion wirdeinDirichlet scher Faktor, nämlich: 


[ 1 

— —1<iI<1l 

re pet. cosltdt= 1 2 z en 9), 
0 lt 

Aus diesem Beispiel folgt der Satz, daß die Summe der Zahlenfolge a 5% , wenn das Vorzeichen 


jedes Summanden gleiche Wahrscheinlichkeit hat, jeden Wert zwischen —1 und +1 gleich 


häufig annimmt. 
12* 
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V. Verallgemeinerung der Verteilungsfunktionen auf den Fall vorgegebener Verteilungen der a, 


Wie bereits erwähnt, lassen sich die abgeleiteten Verteilungen ohne Weiteres auf den Fall 
erweitern, daß die a, selbst irgendwelchen vorgegebenen Verteilungen genügen. Es seien: 


klar) EN ee 


die Frequenzfunktionen der Variabeln a,, dannn bestimmt sich die Frequenzfunktion der Sum- 
menverteilung nach der Formel: 


/2 n 1 
Dlleir were FAN Eve EC (5.2). 


Wir wollen nun an einigen prägnanten Beispielen aus der Literatur, die Nützlichkeit der Formel 
(5.2) zeigen. 

Aus Gl. (5.2) folgt als erstes Beispiel der Satz, daß die „Summen‘-Verteilung eine xX--Ver- 
teilung ist, wenn die a; selbst einer solchen Verteilung genügen. Es seien also die Werte a; 
entsprechend den Frequenzfunktionen: 


P—1 p/2 a a: 
ut a ee 


fl b 
ar ö 


= 


verteilt. Jedes einzelne Integral 


ii 1 
1,0) = | Ann (a) orl2; —atr)aa . . - et 
0 $ 


ist ein Exponentialfaktor: 


Br 
Tg Abpl) Fe) an u ee Se 
so daß die Frequenzfunktion für die Verteilung der ‚‚ Summe“ sich aus der Gl. 
?]2 2 15a n 
ee a (derd d=-Nu:..- .(56) 
ONE r(®) ; Fa re 
9,2 


berechnet, was für die Verteilung von } wieder auf eine Frequenzfunktion nach Gl. (5.3), näm- 
lich auf: 


2 
2 


Ir—1 p p2 1 Pp 
Jp.n (1) = 2) u Eee sy re 
—1 ei ) 
ee 

führt. In diesem Ergebnis ist als Spezialfall für p = 1 ein Satz der mathematischen Statistik 
enthalten, der besagt, daß eine Summe wieder normal um Null verteilt ist, wenn die einzelnen 
Summanden normal um Null verteilt sind. 

Als weiteres Beispiel dafür, daß die ‚Summe‘ einer gleichen Verteilung wie die der Sum- 


manden genügt, wird in Lehrbüchern der theoretischen Statistik ?2) für p = 1 vielfach die so- 
genannte Cauch ysche Verteilung mit der Frequenzfunktion: 


a 1 


ee —o< RT 0 u ne (5.8) 
angeführt. 
Eine damit eng verwandte ist die von R. A. Fisher eingeführte Verteilung mit der 
Frequenzfunktion 
re) ER, 
2 p* 
le z a = er, (> 0). er 
rar) are 


Wir wollen nun zeigen, daß die Verteilung von Fisher für g=1 und für die Variable x? 


statt x, das Analogon zu (5.8) bei p Freiheitsgraden ist, und wie die Cauch ysche Verteilung 
die Eigenschaft besitzt, sich für die Summenverteilung zu reproduzieren. 


2) z.B. H.Cramö6r, Mathematical Methods of Statistics, Princeton 1946 und M. G. Kendall ‚ The 
Advanced Theory of Statistics, Vol. I, London 1948. 


Z. angew. Math. Mech. : . R ; 
Ba. 94 Nr. 4/5 April/Mai 1054 Deuker, Über die Verteilungsfunktionen von Vektorsummen 173 


Wir nehmen also an, daß die a; einer Verteilung mit den Frequenzfunktionen 


nie! 4 
han (= er — ea 
re) rm 


genügen, in denen (5.8) für die halbseitige Verteilung &> 0 für p — 1 enthalten ist. 
Jedes Integral liefert wieder, wie im ersten Beispiel, eine Exponentialfunktion, nämlich: 
RT ee ee re nlier: Mn 


so daß die neue charakteristische Funktion als n-faches Produkt sofort weh ist, wie auch 
die Frequenzfunktion als Hanke l-Transformierte der Exponentialfunktion e-®!; 


2/2 F N 
Jr (l) = Fer J2 (li) ir? di b= 3b RE 

er): | 

2 
bekannt ist und die Funktion: 
pP + n 
2b. \ = P—1 

Jr. »(l) = . (5.13) 


Im r(®) wi 


liefert. Die Verteilung für I hat also die gleiche Frequenzfunktion, wie die einzelnen az. 

Bei der allgemeinen Verteilung von R. A. Fisher ist zwar die Frequenzfunktion für I 
nicht mehr explizit angebbar, aber für die charakteristischen Funktionen besteht ein inter- 
essanter Zusammenhang. 

Es sei 


Arldr a EM Ar ee: 
“ DON (dk + ar) ? 
Man findet dann: 
28 a gl: r(2) * Io_‚(ayt) apl? day 
In(bi) = ER r | z — (5.15) 
rat 
Tg, (bt) = m r Ka (bxt) 


wobei K,(z) die modifizierte Hanke l-Funktion bedeutet. Die neue charakteristische Funk- 
tion ist also: 


E> 1 n 4 ; 
ae e) 0 Kaplan) > > > > + 6.10), 
n\—— q n = 
een 
der die Frequenzfunktion: 
p]2 
Dee : | 7 ] ä | oe? Ras(drt) I» _ (It) WR di (5.17) 
u 


zugeordnet ist. 
Zum Abschluß betrachten wir noch zwei Verteilungen der mathematischen Statistik, die 
in einem interessanten Zusammenhang zu unseren ‚‚Summen‘“-Verteilungen stehen. 
Als erstes nehmen wir die sognannte hypergeometrische Verteilung 
1 


—, . m—1 „gn—1 2 RE RE 
dF TER (1 — ar 1 gr de, zus ( ) 
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die wir etwas allgemeiner und mit der Variabeln & = a; in der Form schreiben: 


nt 


hy,x(0) = TE (bay ..... 019). 
rg: r(2) DER 
Es folgt dann: 
br. 
I u a} Da pn day; 
p,% (bt) E 2,% (4x) 021 IR 4 k k 
Ypl2. r(# zu ) d% 
En 2° | Fl OE—adrı zapl? da, 


Pi re, Fe 


Bu 
9p12. es nr ) a (5.20) 
= — u ik 4 (b,t sin 9) - sin?/?2 9 cos219 dd 
ER Si 
road) ° 
er in ri a ) 
= en RS (b = m 174: 3): 2 
er a Sr a Re 
und die neue Frequenzfunktion lautet: 
]p/2 f n 'p Y b F 
ern | DRS In zb) 7, _ (yw8dı.. .. (5.21). 
rs MEr 2 


Die charakteristische Funktion einer Vektorsumme bei p-Freiheitsgraden und Primitivvertei- 
lung (konstante a;) geht also bei einer hypergeometrischen Verteilung der a, nach Gl. (5.19) in 
die entsprechende charakteristische Funktion von 9 + 2 q Freiheitsgraden über. 

Als letztes Beispiel nehmen wir die Verteilung 


LEN 
hp, (4) — a eyege u . 
In diesem Fall findet man: 
p oo 
98 =; r 2) r 
2). DR 1 
ee N Are | Ten, (ta) a2 day =——— 2 (6.28) 
ATI a+oj® 
und 
a ‚ Jr, pi? dt 
N a 
2? 22) s (1+1) 2 
a’ 0 
jp—1 eng an 
ge u ne; I 2 
a r(3) vn Zi 2 ) ee #7 ns = m 


wobei wieder K,(z) die modifizierte Hanke l-Funktion bedeutet, - 


Eingegangen am 4. Mai 1953, 
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Zusammenfassender Bericht 


Zolotareffsche rationale Funktionen 
Von H. Piloty in München 


Für mancherlei technische Probleme, vornehmlich solche aus der Synthese elektrischer 
Wechselstrom-Schaltungen — Auffinden solcher Schaltungen mit vorgeschriebenen Eigen- 
schaften — spielen gewisse, von dem Mathematiker E. Zolotareff 1877 erstmalig unter- 
suchte (gebrochene) rationale Funktionen eine wesentliche Rolle !). Sie lassen sich in gewissem 
Sinne als eine Verallgemeinerung der Tschebyscheffschen Polynome auffassen, da 
sich die letzteren durch einen Grenzübergang aus jenen gewinnen lassen. Ähnlich wie die 
 Tschebyscheffschen Polynome lösen sie die Aufgabe, den Wert 0 in einem reellen 
Intervall der unabhängigen Variabeln gleichmäßig zu approximieren. Sie unterscheiden sich 
jedoch von ihnen dadurch, daß sie im komplementären Intervall den Wert co gleich- 
mäßig approximieren, was die Tschebyscheffschen Polynome nicht tun. Das 
Null-Approximations-Intervall ist der reelle offene Bereich |z|<1, das komplementäre der 
reelle offene Bereich |2|> 1. Je näher die Intervallgrenzen an + 1 heranrücken, desto schlechter 
wird bei gegebenem Grad die Approximation. 

Funktionen dieser Art — ich nenne sie zu Ehren ihres ersten Bearbeiters Zolota- 
reffsche Funktionen (kurz: Z-Funktionen) — werden in diesem Bericht beschrieben. Sie 
sollen hierdurch demjenigen, der sie benutzen will, in allen wesentlichen Eigenschaften. bekannt 
gemacht werden. 

In erster Linie interessiert der Zusammenhang zwischen Güte der Approximation (maxi- 
male Abweichung A vom Wert 0 im Intervall — Yk <z < + Yk) und Intervallgrenze Yk bei 
gegebenem Grad. Er ist durch Gl. (9) gegeben und in Bild 2 für die ersten 12 Grade 
dargegestellt. Weiterhin braucht man die Konstanten (Nullstellen und Pole), welche die 
Funktionen bestimmen. Sie sind durch GI. (2) und (3) gegeben. Alle sonstigen Eigenschaften 
der Z-Funktionen folgen am einfachsten aus ihrer Parameter-Darstellung (Gl. (5) bis (8)). 
Dies ist der Inhalt von Teil I. Der Teil II beschreibt das Verhalten der Z-Funktionen im ein- 
zelnen und zwar genauer für reelle Werte der unabhängigen Variablen z, in großen Zügen auch 
im Komplexen. Der Teil III schließlich befaßt sich mit einer Aufgabe, die bei den technischen 
Anwendungen der Z-Funktionen immer wieder vorkommt, nämlich: explizite Lösungen der 
algebraischen Gleichung F(z)—=const anzugeben, was für beliebigen Grad n infolge der speziellen 
Eigenschaften der Z-Funktionen möglich ist. Von praktischer Bedeutung sind insbesondere 
die Teilaufgaben, bei denen const entweder eine reelle, oder eine rein imaginäre, oder schließlich 
eine komplexe Konstante vom Betrag 1 ist. 

Die mathematischen Hilfsmittel entstammen durchweg der Theorie der (Jacobi- 
schen) elliptischen Funktionen. Obwohl sie in der Literatur leicht zugänglich sind — vgl. z.B. 
Jahnke-Emde, Tafeln höherer Funktionen, 4 Auflage, Teubner 1948 — habe ich 
doch die benutzten Formeln in Anhang I der Bequemlichkeit halber noch einmal unbewiesen 
kurz zusammengestellt. Eine Ausnahme machen diejenigen Regeln, die aus der Theorie der 
Modul-Transformationen elliptischer Funktionen stammen. Ihr Beweis ist in den Anhängen II 
und III angedeutet. Die Anwendung der Modul-Transformationen auf die Z-Funktionen im 
Sinn von Anhang II (Parameter-Darstellung der Z-Funktionen) stammt von amerikanischen 
Elektro-Ingenieuren (E. L. Norton u. S. Darlington). Diejenige der speziellen 
Landenschen Modul-Transformation auf die Berechnung des Verhaltens der Z-Funktionen 
auf dem Einheitskreis (Anhang III) habe ich nirgends finden können. 

Im übrigen verfolgt die Darstellung die Tendenz, möglichst nur zu beschreiben, die 
Beweise nur anzudeuten und sie dort, wo auch die Andeutung zu viel ablenken würde, in An- 
hänge zu setzen. 


= 


I. Definitionen und Parameter-Darstellung 


Die zu beschreibenden Funktionen F(z) gehören zu einer Klasse gebrochen rationaler 
Funktionen mit folgenden Eigenschaften: (Der Grad sei rn) 


(a) F (z) =eell für reelle 2 
(b) F (1/2) Y/R @) 


(c) Die n Nullstellen a; liegen auf der reellen Achse im Interval_— 1< ag; <-+1. 


1) vgl. die historische Note von W.CauerinZ.angew. Math. Mech. 20, (1940) S. 358. 


m 
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Wegen (b) ist F2(1)= 1. Der Einfachheit halber soll noch normierend 
(d) Fi)=+1 


festgesetzt werden. j 
Die allgemeinste, diesen Eigenschaften genügende Funktion hat offenbar die Darstellung: 


ki N 2 
Fo=- IT az mie dr Lee ae ns 


Die Z-Funktionen haben darüber hinaus noch die Symmetrieeigenschaft, 
entweder: (e) ungerade[F (— 2) = — F(z)] oder gerade [F (— 2) = F(2)] 


zu sein. 
Aus der allgemeinen Darstellung (1) wird dann für u 
ungerade Funktionen gerade Funktionen 
Fe 
=: a, | lies Re a 3) 
mit o<ier jet 
Gradn=2m-+1 Gradn=2m 


Zolotareff hat nun die Aufgabe gelöst, mit Funktionen dieser Gestalt den Wert 0 in 
einem reellen Intervall mit den Grenzen 2) — Yk und Yk, 0<k<1gleichmäßigim Tsche- 
byscheffschen Sinne zu approximieren, ebenso wie dies etwa die wohlbekannten 
Tschebyscheffschen Polynome im Intervall mit den Grenzen — 1 und +1 tun. Im 
Gegensatz aber zu diesen Polynomen approximieren die Z-Funktionen im reellen Intervall 


l2| > re wegen der Eigenschaft (b) den Wert ® gleichmäßig. Daraus und aus Eigenschaft (d) 


erhellt auch, warum sie den Wert O nicht im ganzen, abgeschlossenen Intervall — 1 bis +1 
approximieren können. 


Die Lösung der Aufgabe besteht darin, den Parametern ag, bzw. a@3,_ı der Darstellung (2) 
die Werte 


— 2» ee 
= Vhon| 5 Ei | = Yken| ae K;k Au 


nal)... Mm 
zu geben. Das Zeichen sn steht für die Jacobische elliptische Funktion sinus amplitudinis, 
k ist ihr Modul, K die reelle Viertelperiode, welche mit dem Modul durch das vollständige 
elliptische ‚Integral erster Art 


n/2 
: d 
IX = ame Alec Berg Bu, ae a & ae me en ee EEE (4) 
Yi — KR? sin? 
0 
verknüpft ist. 


Rationale Funktionen mit einer Darstellung nach Gl. (2) und Parametern gemäß den 
Gl. (3) und (4) sind Gegenstand unserer weiteren Diskussion. Obwohl natürlich das Verhalten 
der Z-Funktionen durch die explizite Darstellung mittels Gl. (2) und (3) völlig bestimmt ist, 
kann man es aus einer Parameter-Darstellung viel leichter erkennen. 
Wir setzen unter Einführung einer Hilfsvariabeln v 
z(v)—=Yksn[voK; iii ee 
und erhalten dann für w = F (z) im ungeraden bzw. im geraden Fall: 
(0) = (— 1)"Yhısn[noK;k) | wW)=(- 1" Yasn[(nv-+1)K,;kl. - (6). 


Hierin ist der neue Modul k, gegeben durch 


u a ‚2r—1 uk, 7 
li II” Enns K;h Ne wi 3 ee (7) 


2) Die Wahl der Bezeichnung + /k für die Intervallgrenzen wird sich gleich rechtfertigen, 


Zul 
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und die neue reelle Viertelperiode K, entweder durch Gl. (4) nach Ersatz von % durch k, oder 
auch durch 


3 ars h 
RE r vl 
ee an ER 


PAR ee 
IT| - RK; 


Der Beweis für diese Behauptungen findet sich in Anhang II. 


II. Eigenschaften der Z-Funktionen 


1. Verhalten für reelle z 2 
1.1. Das Verhalten im 0-Approximations-Intervall —yYk< DZ Y% läßt sich unmittelbar aus 
der Parameter-Darstellung (Gl. (5) und (6)) ablesen. Durchläuft die Hilfsvariable v das Inter- 


vall—1bis +1, so oszilliert wzwischen den Werten — Yk,und + /k,, und zwar wird der Werte- 
bereich zwischen diesen Werten genau » ma] durchlaufen, wie es die Prinzip-Skizzen ‚von 
Bild 1 zeigen. 


! 


Bild 1. Grundsätzlicher Verlauf der Z-Funktion für die Graden =1bisn =5 
Die maximale Abweichung vom Wert O.wird also durch die Zahl 
== Mm 2»— 1, 
A= Yk, — (Yk WALSn: re Be RE ee 
vl : \ 
gegeben. ‘Gl. (9) verknüpft die Approximationsgenauigkeit A mit 
dem durch Yk gegebenen Intervall. Bild 2 stellt diesen Zusammenhang 
dar. Aufgetragen ist 1/A abhängig von Yk für die Grade n=1 bis 12. Man sieht, wie 
rasch die Güte der Approximation mit wachsendem n zunimmt. Ist ein Wertepaar A und Yk 
als erwünscht vorgeschrieben, so kann man Bild 2 den Mindestgrad entnehmen, mit dem 
IF (2)| < A im Intervall z| < Yk erreichbar ist. 

Die neinfachen Nullstellen sind diejenigen im Intervall liegenden Stellen 
an denen nv bzw. nv +1 eine gerade Zahl ist, also, wie es wegen Gl. (2) sein muß, die 


Stellen 
z=0- und +@%, z= I My : 


nach Gl. (3); vl... iM. 


Die n—1 Extrema sind diejenigen im Intervall liegenden Stellen, an denen nv 
bzw.nv +leine ungerade Zahl ist, also mit 


= _— 1 — 2v—2 
Gi Yk sn & = Re: ı | Agy_a =Yksn Ei K: |. Se 10) 
die Stellen 
ee, Wei...m) | a=liund ar Ver). 


Die Extremalwerte + A werden auch an den beiden Intervallgrenzen + Jk angenom- 


men und zwar der Wert + Abeiz— + Ykund der Wert (— 1)" Abeiz—= — Vk. Jedoch sind 
diese Stellen keine Extrema, 
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Mittels der Parametera lassen sich dieApproximationskonstante4 
und die reelle Viertelperiode K, (der Gl.(6)) so ausdrücken: 
m 


— YkHa},ı MT a ae (11). 
y=1 y—i 
Mm m 
1 ITaz, 1 IT a v—2 

—1 r vV= 

K,=— KT Bear ie ee ee ee (12) 
Ila; v—1 IT 0} v—1l _ 
v-1 y= 

Gl. (11) folgt unmittelbar aus den Gl. (3), (9) und (10), Gl. (12) aus den Gl. (3), (8) und (10). 


10* 


0 0102.08 0042 000 08 WOLLE 


Ik — 


Bild 2. Aproximationsgüte 1/4 als Funktion von Vr Parameter: Grad n der Funktion (Fz) 


a 


; e 1 
1.2. Das Verhalten im ©-Approximations-Intervallz < — me > + Tr ergibt sich einfach 


aus der Eigenschaft (b), welche besagt, daß, wenn F(z) bei z=z, den Wert c annimmt, F(z) 
an der Stelle z= 1/2, den Wert 1/c annimmt. 
Insbesondere liegen die m einfachen Pole bei 


z—= wo und + 1/as, = + YYay,_ı 
vl.) (VL ER.M) 
und zwischen ihnen die n—1 Extrema 
z—= + 1/ay,_ı z= @ und + 1/as_a 
(A | (lahm), 


wobei reziprok zu einem Minimum jeweils ein Maximum liegt und umgekehrt. 
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Der Wert der Extrema ist + 1/4 und es ist im Intervall stets |F(z)| > 1/4. 
Die Werte + 1/A werden auch an den beiden Intervallgrenzen + 1/Yk angenommen, und zwar 


der Wert + 1/A beiz = -+ 1/Ykund der Wert (— 1)" 1/4 bei z= — 1/fk. Auch diese Stellen 
sind keine Extrema. 


y% 
ebenfalls aus der Parameter-Darstellung der Gl. (5) und (6) entnommen werden. Wegen der 
Eigenschaft (d) (die Funktionen sind ungerade oder gerade) genügt es, nur das positive Inter- 
vall zu betrachten. Man erhält das Gewünschte, wenn man das Argument vK vom Punkt K 
aus in imaginärer Richtung laufen läßt. Die so gewonnene Darstellung, welche die reellen 
Werte von zund w durch eine komplexe unabhängige Variable aus drückt, kann aber durch 
eine einfache, der Theorie der elliptischen Funktionen entnommene Transformation reell 
gemacht werden. Das Ergebnis lautet mit der Hilfsvariabeln % 


1. 3. Das Verhalten in den Ükergangs-Intervallen —. a 2 Ss— Yk und-+ Yk Ss2st = 7 kann 


rent 
Te 


Yk, 


Hierin bedeuten % — 1 — k2 und == 1 — k? die zu k und k, komplementären Moduln, 
K’ und Kj die zugehörigen reellen Weerrelperioden, dn das Zeichen für die Jacobische 
elliptische Funktion delta amplitudinis. 


_ Die Funktionen 2 (y) und w(y) oszillieren zwischen den Werten Yk und 1/Ykbezw.yk, und 
1/Yk,. Sie haben die Symmetrie-Eigenschaft 


Be, dere eh. 


in Übereinstimmung mit der Eigenschaft F(1/2) = I/F(z). Durchläuft y das Intervall 0 über 
1/2 nach 1, so durchläuft z monoton das positive Übergangsintervall Yk über 1 nach 1/yk und w 


monoton die Werte /k, — A über 1 nach 1/4. Gl. (13) und (14) sind einfache Konsequenzen 
der Formel (9) des Anhangs I, die G]. (15) solche der Formel (12) Anhang I. 


2. Verhalten für komplexe z 


Einige allgemeine Regeln für das Verhalten im Komplexen lassen sich sofort 
aufstellen: 


(«) Außerhalb der reellen Achse liegen weder Pole noch Nullstellen, 


da nämlich alle n möglichen Pole und Nullstellen auf der reellen Achse (einschl. 2= 0 und © 
bei den Ben. win liegen. Aus den Grundeigenschaften (a) und (b) folgt: 


B) IF@l21ı für klsı. 
Weiter ist leicht einzusehen, daß außer den Werten O0 und & auch: 


(y) die reellen Werte w mit Beträgen |w| < A und |w| = 1/A nur auf der reellen Achse 
angenommen werden. 


Schließlich: 


(6) Durchläuft (— 1)” z die imaginäre Achse von — soo bis 0 und dann von Obis + vo, 
so durchläuft w: 


beiungeraden F(@) beigeraden F() 
monoton ansteigend die imaginären Werte monoton absteigend die positiven Werte 
—t oo bis 0 und dann +1/Abis +A und dann 
O bis +io. monoton ansteigend die Werte 
+Abis +1/A. 


Alles das kann aus den Formeln (2) und (6) ohne Mühe entnommen werden. 


a . h. 
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III. Lösungen der Gleichungen F(z) = const. 
Bei den technischen Anwendungen der Z-Funktionen kommt häufig die Aufgabe vor, die 
jeweils » Wurzeln der Gleichungen: 
1) F (x) =reelle Konstante 
2) F (z) =imaginäre Konstante 
3) F (z) = komplexe Konstante vom Betrag 1 


zu bestimmen. Die nächsten Abschnitte EREFS die Lösung dieser drei Aufgaben. 
1. F(2)=A(A reell) 

Hier muß man unterscheiden, ob }A| im ee A<|A|<1/A oder außerhalb Has 
Im letzteren Fall sind die » Wurzeln alle reell. 


1.1 jA]|<A=/k.: 
Die Substitution De 
i A pt 
führt die GIn. (6) sowohl im ungeraden, wie im geraden Fall unabhängig von » über in: 
w= +Ykhsn[n&K,;k]. 
Ist x die auf das Intervall O< «< 1/n beschränkte Lösung der Gleichung 
Yk, sn[nz K,;k]= |Al, 


also der Wert: 1 
u ar . “ 4 
ee int °) | 


ee 2 (16), 
op = arcsin (|Al/Yk,) (0<4<#) | 
so sind die n gesuchten Wurzeln 2, =1...n)der Gleichung F(z) = A gegeben durch 
® Irma 
„= Em |(1- 2 1 ya) R;r| Se a 


oberes 


> 0 
Vorzeichen für A 
unteres 


<0, 
Sie liegen im Intervall — Yk<2,<-+yk. In den beiden Grenzfälen A=0 bzw. 


Al= Yk,: — A kommt man — in anderer Numerierung — auf die Nullstellen nach GI. (3) 
bzw. die Extrema nach Gl. (10) zurück. Diese letzteren sind doppelte + A- oder — 4-Stellen. 


Außerdem treten dann noch die Stellen + Yk als einfache Wurzeln auf, und zwar für 


ungerade Fe); n=2m-+1 | gerade F(z);n=2m 

die Stelle 2, — Yk als Lösung die Stellen 2 —yk und „= Pk lu: 
vonF@)=+4, als Lösung von F(?)= +4. I; = 

die Stelle 2, = — Yk als Lösung 


von F()=— A. | 
1. 2 JA|>1/4=1/Yk.- 


Wegen der Reziprozitäts-Eigenschaft (b) von F(z) haben wir hier an Stelle von Gl. (16) 
und (17) einfach: 


it . “ ) 7 
> Er F(k;g,) mit 'o =aresın pr: (0 SsoS S ee) 
und f 
IB 2» —1 } 
VE |(1 re m o)K; ER 20r 
oberes > 0 
n , x 
uüteren | Vorzeichen für A |< 0° 


REIN: = : 
Die — liegen im Intervall — Yk < = ZA Yk. 
2y 24 


3) F(k; o) ist hier und im Fol i fe a MT A € 
= N De Fe) im Folgenden eine a für das unvollständige elliptische Integral 1. Art (nicht 


F(k; 9) | en 
ok — k? sin? 


ENG. \W. 


In den Grenzfällen A—=oo bzw. |A| = 1/yk, — 1/4 kommen wir auf die Pole und die 
Extrema nach II. Abs. 1.2 zurück. Letztere sind dann Doppelwurzeln, zu denen gegebenenfalls 


sinngemäß zu Regel (8) noch die einfachen Wurzeln + 1/Yk treten. 
1.3. 41<JjAa|l<ıj4. 


9 ’ 
Die Substitution v=y—=1— ai +iy ne führt?) die Gln. (6) sowohl im ungeraden, 
wie im geraden Fall über in: - 
4 Er a dk 
w— (—- 1) Yk, nIK, +iyKil= N ee AL 


dn[yKı; ki] 
und zwar für beide Vorzeichen des imaginären Bestandteils von v,. 
Ist daher y die zwischen O und + 1 liegende Lösung der Gleichung 


dn[yKı;kıl = I NR RER u 
deren rechte Seite eine Zahl zwischen 1 und A, ist, also die Zahl 
1 E . 1 —— Y 
ISIN 6 en ae ( en 
u y (k;9) mit) g=arecsin K Yı —.h,/A4 Dirt ee 5 (23) 
so sind die 2 » Lösungen der Gl. F2?(z) — A? gegeben durch: 
zy = Yk sn ) Ktiyk;k| v0, er HER IE (24). 
Derrersteund letzte Index v lielert nur je,eine Lösung, und zwar 
eine reelle, nämlich: 
2 . 
hen ee (24a), 
Zn) \dn [yK’; kW] 
alle übrigen jedoch je zwei konjugiert komplexe Lösungen. 
Die zu Zee Indices » gehörenden Lösungen sind wegen Gl. (21) Lösungen der 
ungeraden | 
A 
Gleichung F(z) = Fi “ Ä 
Die komplexen Zahlen 2, mitv=1...n — 1 lassen sich reell in der Form 
2 =, im - 4 Va ME fan ko! Varel ke ae . (25) 
schreiben. Es ist 
— s, de 
Er — Yk 


2 + (ksıs2)? 
cd s2C2 


+ (ks, 2)? 


Nv — + yk c2 
2 


mit den Abkürzungen ?) 


(beide Vorzeichen) 


Hm 
29 
a on (1 - )R; a]; Tr A OR NN 
29 
a= ml -P)R:a; | 
N 
s—= sn[yK';k]; ca= en [yK';k]; dge= dn[yE‘;k]. ./. \. (28) 
Fry= 04, e1l;4=0)wmr=eneg, =—1; &=0) ergeben sich wieder die 


Ausdrücke (24a). Ist der Grad n=2m gerade, so folgen für v = m zwei rein imagi- 
näre Lösungen, nämlich (wegen ,=0; a =d,=|]) 


En =0; 7m = + YVkse[yK’;k] 2 ee Re re 1268), 


4) unter Berücksichtigung der in Anhang II Gl. (7) abgeleiteten Beziehung: K]/K, = n K’IK. * 
5) vgl. die Formel 9 des Anhang 1. 
6) In den Grenzfällen A?=k, und A?=1/k, wird 
9 =0 undgp =n/2 
y=0undy =|l. 
?) Außer den bereits benützten Funktionen sn u und dn u = Yı — k2sn?u tritt hier auch die Funktion 


ee a F EN U ler 
enu=/1— sn? u (cosinus amplitudinis) auf, weiter unten auch noch die Funktion sc u = om „(tangens amplitudinis). 
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ImVs pezialtfal ei 


1 ee Ei 
ergeben die Gl. (26) (wegen nn + = Ya ; 


die Lösung der Gleichung 
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= yr) 


= (— I 
durch; a 
8 C i ? 
&,=(l+% 1 = ng »„=+} ne (beide Vorzeichen) . . . . (26b) 


für [Al 1 


P=1...%—1) 


wobei & +7, —=1, wie es sein muß, da diese Punkte auf dem Einheitskreis liegen. Für v — 0 


ergibt sich natürlich &, =1; 
2. F(z2)=iB(B reell) 


n=0 und für ah Fer =D: 


Die Imaginärstellen von F(z) findet man durch die Substitution von 


A — 


1 R 


a 1 ea a KR 


in die Gleichungen (6). Man erhält: 


w— Yh, sn iyKi;kl=iyksclyKi;ki]l - - - - 297 
Ist daher y die zwischen — oo und + oo liegende Lösung der Gleichung 
se[yKi; kı] = B/Yh, (30); 
deren rechte Seite eine beliebige positive oder negative reelle Zahl ist, also die Zahl ®) 
1 , D az 
y-zgrfli.n mit p=artgByi; (-5 <p<+f) . (31), 
so ergeben sich die n Lösungen der Gleichung F(z) = iB aus: 
en —1 
= Em |1 ee JR DPtiyk;h) W=1:-..n). . . (82). 


+i n=1 -1 


n=2 
m=0 m=1 


F(z) positiv reell 
F(z) negativ reell 
F(z) positiv imaginär 
F(z) negativ imaginär 


o° Nullstellen 
x Pole 
Bild 3a. Verlauf der Kurven in der komplexen z-Ebene, auf denen 


F(z) rein reelle und rein imaginäre Werte annimmt, für die Grade 
r=1bisn = 6im 1; Quadranten des Einheitskreises 


von F(z) 


Für B=0;y=0 ergeben sich wieder die 
reellen Nullstellen, für B=»; y=-+1 die 
reellen Pole der Funktion F(z). Für end- 
liches und von Null verschiedenes B hingegen 
sind alle z, komplex. 

Die komplexen Zahlen z, lassen sich 
wiederum reell in der Form 


»=& in = Teer 
schreiben. Es ist 
= s,da } 
ee ne, 
Se % + (ksıs) | ... (34) 
Br ıyp-ı /n Cd, s2c2 
ir 2 + (ks,s2)? 
mit den Abkürzungen °): 
pl 
ea (( nen N RK; |; 
N 
a (i Ay ') K: | 
= ||1-———)K;k); u} (35), 
= 
d, = dn (( 5 u) K; x 
N 
2 = snlyK5 ki en yRE El \ (36 
dga — dn[yK'; k'] J ) 


8) Ist B negativ, so rechnet man zweckmäßig mit 
|B| und kehrt dann das Zeichen von y um. 

®) Bei den Gl. (36) ist zu beachten, daß y 
positiv oder negativ sein kann und daß sn (— u) = 
— snu; en (—u) = cnu; dn (— u) = dnu. 
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Ist der Grad n=2m-1 ungerade, so liefert der mittlere Indxv»—=m-+1 

eine rein imaginäre Lösung. Dann ist nämlich , = 0; c, =d, =1, und es ergibt sich: 
Em+1 —) und Nm+1 = (— 2% SC mK'; k] Mr er (34a). 
Bild 3 zeigt grundsätzlich den Verlauf der Kurven in der z-Ebene, auf denen F(z) rein reelle 
und rein imaginäre Werte annimmt, für die Grade n = 1 bis » = 6 nur im 1. Quadranten des 
Einheitskreises, für n = 5 und n = 6 dann noch ausführlicher in der ganzen komplexen Ebene. 


. 


— 


Sa 
> —— m 


— 


Bild 8b. Verlauf der Kurven in der komplexen z-Ebene, auf denen F(z) rein reelle und rein imaginäre Werte annimmt, für die 
Graden=5u.n=6 Zeichenerklärung siehe Bild 3a 


128 == 
Im Spezialfall B=-+1ergeben die c1. 84) (wegen v3; 4= +4: 
ae 
2] See 
s ned 
N a Pr ED. 


wobei wieder & +; —=1 ist, die Abkürzungen nach Gl. (35) zu benutzen sind und in der 
Gleichung für n, das obere (untere) Zeichen zu nehmen ist, wenn B= + 1(— 1) ist. 


A 2 th. h. 
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3. F(z2) =eip ” 


Die auf dem Einheitskreis gelegenen n Stellen ,= ed» =0,1-:-n—1),andenen F(z) 
den vorgeschriebenen Wert €? annimmt, findet man durch die Suhztiflffon von 


% 1 „Ari nl 
e PREUZ wur 2K 
in die Gleichungen (5) und (6). Man erhält zunächst RT im ungeraden wie im geraden Fall): 


er rkn|1--")R+s Re; DE 
er = w— Ykısn ar + 5 Kl = ,. (38). 


Die Gl. (37) und (38) können mittels einer Modul-Transformation reell geschrieben werden (Be- 
weis vgl. Anhang III) nämlich: 


TC 4v 
el +2 1)8.;%) 2 ee a ee 
9= 5 + am ma—1)Ky;k] EEE ne 0 ee © . (40). 


Die hierin auftretenden neuen Moduln %, und k, sind gegeben durch 
Ay und k,— 2yk 


118 === 1er kı 13; Te De ee er © (41). 
Die zugehörigen reellen Viertelperioden K, und K, durch 
R-=- (1 %)X und R,=ll1 FE Der (42). 


Das Funktionszeichen am bedeutet die Ja co bissche (elliptische) Amplitude. 
Wählt man p aus dem Intervall -—na<g < +, so wird durch Gl. (40) diesem Wert ein 


im Intervall — - <es+ - liegender Wert von & eindeutig zugeordnet. Man erhält ihn 
aus der Auflösung von Gl. (40) zu | 
1 1 7 
»- 1 gr — ) ne 5. . (43), 


worin n 


h dr 
En )= — ,,..2. 0.0.0.0... 0. (44). 
BR ie Yi — k3sin2r (4) 
das unvollständige elliptische Integral 1. Art zum Modul %, ist 1%). In den Gl. (43) und (44) 


kommen Winkel zwischen — 3 2 und + z vor. Das ist unbequem. Man benützt daher 


besser für die vier Quadranten verschiedene (aus den Eigenschaften der Funktion am leicht 
beweisbare) Formeln: 


a RT Se 
et Bu 
= |! = F(k,; y) Sa + Eh ) 
ul — e —o® mit y=o— = 
& (45) 
7 7 
3. Quadrant:- n<ps— 4. Quadrat: - Tape 
5 1 1 Re 1 1 
Zt "Fk o) == Plee)) 
mit y=—9— I mit y=I+9. 
Der inihnen vorkommende Winkel yliegt in allen Fällen zwischen 0 und Eu 
10) Die F 1 (43) 1äß 
) Die Formel (43) läßt sich mit F, (k,;Q) = (= re — (44a) auch Inn Fu (k,; 9) (43a) 


schreiben. 


| 
| 
| 
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Ist # aus einer der vier Gleichungen (45) bestimmt, so folgen die gesuchten n Werte von 
9, aus der Gl. (39), wenn man » n aufeinanderfolgende positive oder negative ganze Zahlen 
durchlaufen läßt. Zweckmäßig zählt man von v» = 0 bis an — 1, wenn >0 und von» =: 1 bis 
n, wenn z<0. Dann liegen alle n Werte 9, im Intervall0O <9,<2%r 

Für die numerische Rechnung ist noch anzumerken, daß am(— u) = — am (u); 
. am (u+2K)=am u-+rist. Damit kann man alle Berechnungen auf solche des 1. Quadran- 
ten zurückführen. 


"Spezialfälle. 
Von besonderer Bedeutung sind die Spezialfälle: 


a a. here lo en a N 


und 


P=+75 a ee Be a ih 


Sie sind in den Abschnitten 3.1 und 3.2 bereits behandelt. Dort sind die zugehörigen, auf dem 
Einheitskreis liegenden z-Punkte in rechtwinkligen Koordinaten berechnet worden. Hier er- 
halten wir die Lösung in Form des arcus. Fassen wir die Lösungen der jeweils zwei, unter (&) 
und (ß) formulierten Aufgaben zusammen, so liefern die Gl. (39) und (45) mit leicht geänderter 
Zählweise das Ergebnis: 


"Aufgabe (oa) 
DE ER EL TEN TEN 15) 
mit 
= +am | 1) Rh, a N 


Es genügt wegen der Symmetrie der Lösungen m Zahlenwerte zu berechnen (n—= 2m -+1lim 
ungeraden, n—=2 mim geraden Fall), nämlich die im 1. Quadranten liegenden Winkel 9 * * * Öy. 


Es ist 
=! 
TC 2» 
eye am I-#)2:%| (em) 
Day -- DD, | ER RO (48). 
Di „=1::-"m) 


Inn Hd 
Rt = an Dr Di 

Im ungeraden Fall erhält man durch Gl. (48) genau 2n verschiedene Werte, im geraden 

Fall zwei Werte doppelt, nämlich m = On = u und Ontm = Den m = 3 = ‚ also wieder 


genau 2n verschiedene Werte, wie es sein muß. 
Aufgabe (ß) 
mit 
2v—1 = 
= - + am je + ) R,; A BEN In) 2. ee 150): 


Auch hier genügt es, m Zahlenwerte zu berechnen, nämlich die im 1. Quadranten liegenden 
Winkel 9, Bu. 


Es ist 
2v—1 
%, = - — am (i Bay 2) R;: | (a Tr: m) 
Gay m 0, ] 
Et | a a er 1 BF 
Don yı1 —=2n — © | 


16) 
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Gl. (51) liefern 4 m verschiedene Lösungen. Das sind im geraden Fall alle. Im ungeraden Fall 
fehlen noch zwei, nämlich 97 +1 und Ögm4+2. Diese sind 


7 
Öm+i = By 


(ungerader Fall) solar 


1, 
Oamı2 = 3 6 
Anmerkung. Die numerische Auswertung der in dieser Arbeit enthaltenen -For- 
meln wird sehr erleichtert durch die Tafel von E. Glowatzki!). Seine Tafel enthält die 
Jacobische elliptische Amplitude als Funktion rationaler Bruchteile der reellen Viertel- 
periode K und des Modulwinkels © = arc sin k, also die Funktion 


p= am E K; sin 0) 


und zwar in dezimal unterteilten Graden der 90°-Skala mit 3 Dezimalstellen fürn —= 2,3: ::12; 
-m=1,2::-n— 1 und für ganzzahlige Werte des Modulwinkels © von 0 bis 90°. 

Die elliptischen Funktionen sn u= sin und cn u—= cos o können dann ausden Glo- 
watzkischen Tafelwerten mit Hilfe passender trigonometrischer Tafeln, z. B. denjenigen 
von J. Peters2), unschwer berechnet werden. 


Anhang I 
Einige im Hauptteil benützte Formeln aus der Theorie der Jacobischen elliptischen 
Funktionen ®) 


1.Moduln und Perioden-Konstanten 
Modul k; komplementärer Modul k — Yl — k? 


2 d 2 d 
T T 
ee ee ER Nr. 
I Van: n 
ö ö 
2.Die drei Funktionen sn,endn undihre Umkehrungen 
Es sei 
€ d 
J Yı— Resin®r 2). 
ö 
Dann ist: 
amu—p, u—=F(k;o) mit: } 


= SNH—=SInO P=sAaTC sing, 


L = CNU==COSP O—IATC COS T, 


[| _ 0 ——_ 


%=dnu= yı — k2 sin? p; p=arcsin en en 
sn u 
a p=arctg x, 
3.Spezielle Werte: 
(0) 0; en(0) =L; dn (0) — 1 
K Sen N K Fr 
sn Der- V en STR | DLR dn Eas Yk en ae (4), 
snK=]; MA 0: din K=Kk | 
WK' 0 iK' 1% yo a | 
SM —— — —; EN. = er N — — Kk 
2 Yk =») \ jr ’ dn 9 S= yı + k 5 (9): 
Er sn iK’ — oo; en iK’ —= oo; dniK’ —=o ) 


m 
= x ‚ Abhandlungen 


ae Akademie der Wissenschaften, Mathematisch-naturwissenschaftliche Klasse, Neue Folge 
12) J. Peters (siehe Literaturverzeichnis bei 11). 
13) Alle diese Formeln stehen entweder direkt in „Jahnke-Emde ‚ Tafeln höherer Funktionen, 4. Auf- 
lage, Teubner 1948°“ oder können aus dort stehenden durch leichte Spezialisierung gewonnen werden. a 


U)E.Glowatzki, Tafel dr Jacobischen elliptischen Funktion p=am 


2. R . b | 
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4. Periodizität und Symmetrie-Eigenschaften der Funktion »n 
sn(-uW)=— snu, 
n2 KL W)=F snu, 
sn(4K+u=-+tsnu, 


ER 1 a 
sn(K Lu)= et | | 
sank Luy=+tsnu 


)) 


5.Die Funktion sn für komplexe Argumente 


Allgemein: 
= .\  S1da + Ücıdı s2c5 
sn (2 +iy) = a N) 
mit 
Seen (z;K), sen (U, k), | 
a=en(s;k), @=cn.(9;%); le a a Tal 
dh=.dn(x;k), d; = dn (y; k') | 
Speziell: 
a) ==; (= RER ENT A REF Ha NUS), 
c2 
b) SR; SE re 
2 
K' r | 
c) y=5; sn (e+ 3:1) 7, sp + isn A rat ai): 
mit 
cosy—=(1+k) sin y —= re! (10a) 
RT ee ITTris DR mL 
6.Zwei Additions-Theoreme 
sn” U — SM? uU, 
or (rt een en IR 
K K ‚ 
ano \an(S +2) 3 ea Mr (12), 
7. Ableitung der Funktion sn 
smu=onu:dnu EN EEE RE TERERT LS 
8. Eigenschaften der Funktion am 
am (— u) = — am u, am2K tur tamu. »..00. (14), 
Anhang II 
Beweis der Parameter-Darstellung von F(2). 
Setzt man 
2— Yksn [WK; k] a A 
und 
n—=2m +1 (ungerader Fall) n—=2m (gerader Fall) 


et — [2v—1 
Ag,„—Yk sn ® K; ı W=1:...m). . (2a) | a,_Yı= Yksn en K; (=1::: m) (2b), 


so erhält man für 


Fe) re, en ? re 0,1 
—= =? te er Tg — —— 7 SEE 
= £ I 1— @,2 II 1 - a v—1*? 


13* 
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mittels Formel (11) des Anhangs I 
Sl 
N; w— ki? I“ (v En a) %: r 5 
= KR II le +Z 22) 8; ı 
a | on (or; 4 
oder mit Abänderung der Numerierung in beiden Fällen: 
n e ae] Er | 
DEE e-1 ee, ET N EN, 


Das Ergebnis (1) und (3) schreibt sich etwas einfacher mit v=v — 1zu: 


eben (ul) Rh 2 0 ae u el 
n = A Ey 
w=hk > Toll + 7) K; 1 a ee ee 


4 
w als Funktion von w ist nach Gl. (5) eine elliptische Funktion mit der reellen Periode = und 
der imaginären Periode 2:2 (x = K’/K) und den in einem Periodenparallelogramm liegenden 


1 3 5 3 
Nullstellen I und > und Polen — +:K und — +1iK. 
N N N N 


Dasselbe trifft zu von der Funktion 


wlan nu TV) Ra). 2 ne nee 
in welcher der Modul k, bestimmt ist durch das Periodenverhältnis 
Kı K' 
er! Kun en 


und K, das zugehörige vollständige elliptische Integral ist. Die Funktionen w, und w können 
sich als elliptische Funktionen mit den gleichen Perioden und den gleichen Nullstellen und Polen 
nur durch einen konstanten Faktor unterscheiden, sind also identisch, wenn man in Gl. (6) den 


Faktor C passend bestimmt. Setzt man in Gl. (5) und (6)u=i = ‚ so erhält man 
v1 und w= (/yh. 


Ersteres, weilin dem Produkt jeder der Faktoren nach Gl. (10) Anh. Iden Betrag 1 hat und nur 
paarweise konjugiert komplexe Faktoren (im ungeraden Fall auch einer vom Wert +1) auf- 


treten. Letzteres nach Gl. (10) Anh. IT unter Berücksichtigung von Gl. (7). Es ist also C—= Yk, 
zu setzen und damit 


w— Ykısn[(nu-+1)KR,;k] 2 ee SEES: 


Modul %k, und Periodenkonstante K, können auch direkt bestimmt werden. Nach Gl. (8) und 
(5) ist nämlich : 


w(0) = Yk, = — 18 Ion | — K; | 


oder: 


Yk, Pr Ton IR; |: EEE in, = 


Die Konstante K, ergibt sich aus dem Vergleich der Werte der Ableitung 2 an der (Null-) 
1 
Stelle — “ mittels der Formeln (13) und (4) aus Anhang I. Es ist 
Kl, 


1 I 
nach Gl. (5) w' = = =.K SEE, 


2v R z E 2v 
— K;kl=K kul2 sn? | — 


Are DakK 
nach Gl. (8) w - = — Yk, nK, 


De a a 
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und daher unter Verwendung des Ausdrucks (9) 


Ken |) & k 
K=7 ee . . . . a) Ze . . .. . (10). 
a—i sn? Erz T2C: k 


Ersetzt man in Gl. (8) dieses Anhangs wieder u durch v — 1, so folgen die Gl. (6) des Hauptteils. 
Die Gl. (9) und (10) dieses Anhangs sind die Gl. (7) und (8) des Hauptteils. 


Anhang II 
Verhalten der Funktion 2 = Yk sn [vK; k] 
auf der Geraden u— i5 «= KK). 
Nach Gl. (10) und (10a) des Anhangs 1 ist 
z—=Yksn (a+:%) x: Hl c0s 0 + isin d rg ee 
mit: 


sn [xK; k] } 
1+ksn[xK;k] 
en [&K; k]- dn[xK; k] | 4 
1+ksn?[xK;k] ) 


g9(2)=cos®—=(1+k) 
an, 


hi2)—=sin A 


Die Funktionen g(x) und h(z) sind elliptische Funktionen mit den Perioden 4 und ?x und mit 
denselben Perioden, Nullstellen und Polen wie die Funktionen sn [@K;; k;] bzw. en[xK;; kz], 
wobei k, der zum Periodenverhältnis 1*) 

Real ut 
RE EIERN RE (3) 


%; 


gehörige Modul und X, das zum Modul %k, gehörige vollständige elliptische Integral ist. Die 
Funktionen g(z) und sn [xK;; ks] einerseits, h(x) und cn[xK,;; k,] andererseits können sich da- 
her nur um einen konstanten Faktor unterscheiden. Der Faktor sowie Modul k, und Perioden- 
konstante K, ergeben sich durch Vergleich der Funktionswerte für —= 1; = 1/2 und der 
Werte der Ableitung für x = 0 bzw. 1. Das Ergebnis ist: 


9(z) = sn [xK;; k;] == sin am [x K;; k,] \ 


h(x) = en [xK;; k,]= cos am [xK;; k;] | { Zi 
mit = 
2yk 
ka ER EN EIER 
K 75 7 FR, + Dr, 
(3 — (1 Fi k) K ’ K; = ot K . . . . . . . . . (6). 
Vergleich der Gl. (2) und (4) liefert 
T_9=am Re EEE) 


2 
Wendet man diese Regel auf Gl. (37) und (38) des Hauptteils an, so erhält man die dortigen 
Gl. (39) bis (42). Die Wahl der dort angewendeten Substitution bezweckt, daß für = = 0 die 
Werte 9, = 0 und 9= 0 herauskommen. 
14) Dies ist die sogenannte „Landen-Transformation‘“. 


Eingegangen am 7. März 1953. 
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KLEINE MITTEILUNGEN 


Statistische Prüfverteilungen endlicher 
Spannweite. 


Zur Beurteilung von statistischen Verteilungen $(%) 
bei stetigem Merkmal mit annähernd symmetrischem 
Verlauf benutzt man üblicherweise die Ga uBßsche 
Normalverteilung M 

oo 


hm ta=aje Pdı=1 (), 


D®D(2 —o,h)= — 

x 

die man der gegebenen dadurch zuordnet, daß man & 

durch den arithmetischen Mittelwert festlegt und das 

Streuungsmaß % durch die Forderung bestimmt, daß 

die Streuungen beider Verteilungen übereinstimmen 

sollen. Die beiden wahrscheinlichkeitstheoretisch be- 
gründeten Bestimmungsgleichungen heißen dann 


5 p(a) (2 — 0) da —=0, 


+ © 
f ydı=ıl (2). 


+ 

S (a) 1-2 (2 —o))da=0, 
In Verallgemeinerung dieses Verfahrens erhält man eine 
bequeme Möglichkeit der Zuordnung einer beliebigen 
Prüfverteilung an die gegebene dadurch, daß man ver- 
langt, daß ebenso viele Momente beider Verteilungen 
übereinstimmen sollen, wie Parameter in der Prüfver- 
teilung auftreten. Auf diese Weise werden z.B. die 
Pearsonschen Verteilungen festgelegt. (R. v. 
Mises, Wahrscheinlichkeitsrechnung und ihre An- 
wendung. 1931, S. 269.) Dieses Bindungsprinzip für 
die beiden Verteilungen ist jedoch nicht das einzig mög- 
liche. Wir wollen hier ein anderes Prinzip verfolgen, das 
sich an eine bemerkenswerte Umdeutung der Bestim- 
mungsgleichungen (2) anschließt und auf ein Maximum- 
prinzip mit Hilfe der Sehwarzschen Ungleichung 
hinausläuft, das in seiner Geschlossenheit, wenn viel- 
leicht nicht für die Praxis, jedoch für die Theorie inter- 
essant erscheint. 

Man kann nämlich die Gleichungen (2) als Bestim- 

mungsgleichungen für das folgende Maximumproblem 
ansehen: 


va o(z)In d(z2—o,h)de=F(e,.h) =Max.. . (8). 


Indem man nun eine solche Auffassung verallgemeinert» 
kann man einer gegebenen Verteilung eine beliebige 
Prüfverteilung dadurch zuordnen, daß man ein ge- 
eignetes, (3) entsprechendes Prinzip angibt, mit dessen 
Hilfe die Parameter der Prüfverteilung bestimmt wer- 
den können. Man verzichtet dabei zwar auf einen wahr- 
scheinlichkeitstheoretischen Hintergrund, kann aber 
doch von einer wohldefinierten Anpassung der Prüfver- 
teilung an die gegebene Verteilung sprechen. 

Ein brauchbares Prinzip zum Vergleich zweier Ver- 
teilungen läßt sich aus der Schwarzschen Un- 
gleichung gewinnen. Schreiben wir nämlich mit nor- 
mierten Funktionen, die nicht negativ sind 


+» 2 +9 
( Ve yVP@ a, h) as) < [ va) de- 
e " BR 
wi Ya—o,h)d=1,) =| |, 


wobei vorausgesetzt wird, daß die hier auftretenden 
Integrale existieren sollen, so gilt das Gleichheitszei- 


chen nur für =. Für alle anderen Fälle läßt sich ' 


die Prüfverteilung Y/(x — a, k) inihren Parametern so 
wählen, daß 


wo yr (e —0,hde=1-(&,h= Mar.| 4, 


&& (a, h) = Min. 
eist dabei ein Maß für die Güte der Anpassung. 


Ein solches Anpassungsmaß läßt sich für die Prüfung 
durch eine Ga u 8 - Verteilung nach (3) nicht angeben, 
da das Maximum nicht eingeschrankt werden kann. 
Man benutzt stattdessen die g?-Methode. 


Nach dem neuen Prinzip lassen sich natürlich auch 
die Parameter einer Gau Bß- Verteilung bestimmen, 
nach (4) aus der Maximumaufgabe 


ERTEIT: _ (e— 0)? 
Vo (2) ZB dx = Max (0, h), 
 y% 


d.h. aus den Gleichungen 


.— ne (2—0)? 


we (2 — 0) de=0, 


= (A)? 


+» 
iR Vp(a) (1 — 2 (x —o)?) e de=0, 


die sich iterativ lösen lassen etwa unter Zuhilfenahme 
der Werte &, und h, aus (2). Die Durchführung ist je- 
doch unbequem. Wählen wir dagegen die folgende 
Prüfverteilung 


. (1 — h?(x — o)?)%, 


Y(x2—0,h = = 


1 . (5), 
——<rt<a +— | 
x iS 7 h 


die inden Pearsonschen Typen enthalten ist, hier 
jedoch anders zugeordnet wird, so ist 


1 
+, 


15h 
\ eV Te (1—-M(a—oP)de=1— 


Mn 


zum Maximum zu machen durch Bestimmung der Para- 
meter aus folgenden Gleichungen: 


. (6). 


+ 
| Vera - 5%(& — 02) de = 0 


Wir finden demnach für diese Prüfverteilung ähnlich 


einfache Bedingungen wie früher in (2), nur daß Vp 
anstatt pin den Integralen auftritt. 


Für die Rechnung ergeben sich infolge der a, h-ab- 
hängigen Grenzen keine Schwierigkeiten, wenn die 
Spannweite der Prüfverteilung größer oder gleich der 
Spannweite der gegebenen Verteilung ist. Man be- 


stimmt also für /p den Mittelwert %, und mit diesem 
das Streuungsmaß h,, ohne Rücksicht auf die Grenzen, 
für das Gesamtdefinitionsintervall von /p und hat die 
Rechnung a wenn die Verteilung im Intervall 


il 
<NM- 7 + e > enthalten ist. Andernfalls hat 
1 = 


man mit den gefundenen Werten iterativ weiterzu- 
rechnen, bis das Verfahren zum Stehen kommt, in der 
Praxis nach etwa zwei Schritten. 


a en a nn un 
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Mit diesen berechneten Werten & und % ist die Prüf- 
verteilung festgelegt und der Anpassungsfehler be- 
stimmbar. Es pr, 


frei 


u. 


— RR (x — 0)2) de 


und wegen der zweiten Bestimmungsgleichung (6) 


“+7 
15h zen] 
0=- VE N el (a j) 2 
E 
h 
ergibt sich 
1 
Ar 
3h nn 
u j Vode ..... (m. 
a! 
Ve 


Man hat demnach nur noch das Integral der gegebenen 


Verteilung Ve im Intervall der Prüfverteilung zu be- 
stimmen und erhält eals Gütemaß der Anpassung. 


Für folgende Beispiele findet man die Ergebnisse: 


1 
1. o(x) = 5. DES! 
0e >, 


$ 57/3 3 2 
Yz2—o,h),= I6 ee, ‚„e=0,44. 


3 
2. (a) Fa, el=1, 


0% «l>T, 
15 1/4 AN) 
V(e—a, h)= 16 5 we „e= 0,24, 
& (1 — 2], je[s1 
3. o(2) =, 2 SFR 
0 le! >1, 
ja. eure) 
Ya@—a,h-76 g 1 ri ‚„e= 0,28. 


4. Handelt es sich nicht um stetige Verteilungen, son- 
dern sind die relativen Häufigkeiten in einer Werte- 
tabelle diskreten Merkmalspunkten zugeordnet, so sind 
die in den vorliegenden Gleichungen auftretenden Inte- 
grale durch Summen zu ersetzen, die in der üblichen 
Weise ausgewertet werden: 


ale s 8 
x; |pi a Ir ea) ı I 1 Ir (0? 
= g S 
R E 
2|21|0,154| 0,392| 0,7834| —2 | 4 1,568 
3|31|0.231| 0480| 1440| —ı |ı | 0,480 
4,5 | 0,385| 0,620 | 2,480 010 0 
5|210,154| 0,392| 1,960 | +1 | ı | 0,392 
6 | 1 0,077\ 0,277 | 1,662| +2 | 4 1,108 
13 2,161| 8,326 |, =4| | 3,548 


. 1 2 Zn (ao)? 
Zn Zr; 
_ ZVri (ai — 20)? 


Zr; 


(2 8); 
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1 3,548 
— ee 2 2 — 
385g = 1,02 
1 
m —=2,839; 
15 1 © — 3,85\2\? 
w eh nat A 3 . 
, za5 |! | 2,85 )); 


| -y% Zr; = 0,0085; & = 0,30. 


Es ergibt sich in diesem Falle, wegen &, — N = 1,00; 


hı 

ai = = 6,70, daß die Spannweite der Prüfverteilung 
1 
die gegebene ganz enthält, so daß sich eine Iteration 
erübrigt. Läge dieser Fall nicht vor, so hätte man die 
herausfallenden Werte in der Tabelle fortzulassen und 
die Rechnung iterativ zu wiederholen. Wären dabei 
etwa zwei Verteilungen möglich, so wäre diejenige mit 
srößtem h auszuwählen. 

5. Nach dem vorliegenden Prinzip läßt sich eine 
Gaußverteilung ® auf eine Verteilung Y mit end- 
licher Spannweite umrechnen. 


Setzen wir 
Dee 
A 
Ye —o wer h? 212 
’ 7 16 > (2 —.0)?) E) 


so ergibt sich nach einiger Zwischenrechnung 


oe 


an &e>0,29, 


so daß sich die Streuung o der Ga u ß verteilung durch 
h ausdrücken läßt. Man findet 


3 1 
Y2H 


36= 


IR 


und erkennt, daß die 3o-Spanne der Ga u ß verteilung, 
innerhalb der 99,73% des Kollektivs vorhanden ist, un- 
gefäbr zusammenfällt mit der 1/h-Spanne der Vertei- 
lung mit endlicher Spannweite, die das gesamte Kollektiv 
einschließt. Für die Praxis der Anpassung einer Ga u ß- 
verteilung nach dem üblichen Verfahren bedeutet dieses 
Ergebnis eine Stützung der Auffassung, daß man die 
Anteile außerhalb der 30-Spanne als unwesentlich 
fortlassen kann. 

Wenn auch die Bestimmung der Parameterwerte 
nach dem genannten Verfahren nicht die gleiche Ein- 
fachheit erreicht wie für die Ga ußsche Normalver- 
teilung und Bedenken bestehen, das wahrscheinlich- 
keitstheoretische Prinzip durch ein Anpassungsprinzip 
zu ersetzen, so ist immerhin ein Verteil darin zu sehen, 
daß die zugrunde gelegte Prüfverteilung endliche Spann- 
weite besitzt und das Prinzip gleichzeitig noch ein sum- 
marisches Anpassungsmaß liefert. Es ist möglich, durch 
Aufnahme weiterer Parameter auch Schiefverteilungen 
mit endlicher Spannweite darzustellen, z. B. in der 
Form 


Ve —0o,h,k) 


k 
= ml! + h(2 — 0) (1 — k(x — o))*, 
1 1 
nr’ TI 


Die Bestimmungsgleichungen für die hier auftretenden 
Parameter &, h und k sind naturgemäß verwickelter als 
für symmetrische Verteilungen. 


Hannover. RK. Taeckel. 
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Ergänzung zu dem Aufsatz „Die mathema- 
tischen Grundlagen zu einer Meßmethode 
des Schubmoduls zäher Flüssigkeiten‘. 


Obiger Aufsatz ist in der „Z. angew. Math. Mech.“ 
31 (1951), Seite 83—92, erschienen. Die Beschreibung 
der Apparatur und des Arbeitens mit ihr findet man in 
der Arbeit vonH.Fromm und W. Hartung ‚„Mes- 
sung des Schubmoduls zäher Flüssigkeiten“, Die Tech- 
nik, 7. Jahrgang (1952), Seite 126—130. 


In meiner Arbeit sind die natürlichen Koordinaten 
(Zylinderkoordinaten) näherungsweise durch karte- 
sische Koordinaten ersetzt worden, um die Rechnung 
zu vereinfachen. Hier sollen nun die durch Anwendung 
der Zylinderkoordinaten sich ergebenden Verbesserun- 
gen der Gleichungen und Endformeln angegeben wer- 
den. Der Radius bis zur Innenseite des Spaltes von der 
Breite ö sei r,. Die Entfernung eines Punktes im Spalt 
vom Mittelpunkt der konzentrischen Zylinder sei 


r=ent+%® 


Um Anschluß an die Rechnung der Ausgangsarbeit zu 
finden, rechnen wir mit «statt r als Variablen. Die hier 
gewählten weiteren Bezeichnungen entsprechen denen 
der Hauptarbeit. Ebenso erhalten hier die Formeln die 


gleichen Nummern wie die entsprechenden der Ur- 


sprungsarbeit. Doch werden die Formeln dieser Arbeit 
zur Unterscheidung mit einem Stern versehen. 


Die den partiellen Differentialgleichungen (3) bis (5) 
entsprechenden Gleichungen entwickelt man aus den 
Formeln (5.3) und (8.2) — hier muß man T', gleich Null 
setzen — meiner im Ing.-Arch. 19 (1951), S. 271, er- 
schienenen Arbeit ‚‚Stoffgesetze und Wellen zähelasti- 
scher, isotroper Medien“. Dazu drückt.man die Ein- 
heitsvektoren im Zylinder-Koordinatensystem durch die 
Einheitsvektoren im kartesischen Koordinatensystem 
aus und berechnet die Tensordifferentialoperatoren der 
erstgenannten FEinheitsvektoren. Dann entwickelt 
man die partiellen Diff.-Glen. ähnlich der Entwicklung 
der Wellengleichung (9.5) der vorher erwähnten Arbeit 
im Ing.-Arch. Man erhält (Die Korrekturglieder werden 
in allen folgenden Formeln durch ein Rechteck um- 


schlossen): 
a le v n 
re e —) ll +7 er 
a „a ov v 
Aa fr p e ee TR (4*), 
TEN IB a o$ 
da’ nte Ta; 9), 
N) =u ) ee (6*), 
OL U) ro RW, er (7%). 


Das linearisierte Randwertproblem — Vernachlässi- 
gung von o/@ neben 1 — lautet dann: 


@r 1 On Gi Da 
ee Ta 9\de Ira + =08 Mta% 
(89), 
u [Or (0,1) 2 | 
Fr -( a 0). (9%), 


(6,0) 


—=opao? cosal (10*), 


2 
nr 


Durch den Ansatz 
= f(o)-sinat+ g(a)-cosat.. (11*) 


erhält man aus (8*) bis (10*) das Randwertproblem 
gewöhnlicher Differentialgleichungen 


5 1 ae 173 
er c +(@ a) 
+9) = 9 (12*), 
"a ++ om 
ee GE nr)? 
&0 ge: 
ie ö TR) >07. SERIE: (13*), 
1-&)o-4ro; 
rg 0 i 
— (149), 
2u [ 
1— | — =-—f0 
| ER a f(0) 
b A Es ) 
WER OR: 
ee \:- (15#), 
2 f 
2.0) en PA), 0 | 


Da r, stets sehr viel größer als & ist, läßt sich das 
Differentialgleichungssystem vereinfachen, indem man 
für 


1 1 z 1 1 x 
———_— Il — — R 
rg +% rg | 2) (r9 +)? 13 ( 22) 


setzt. 


Da das Glied x/r, bzw. 2 x/r, neben 1 vernachlässigt 
werden kann, weil x/r, und 2 z/r, in den eingekastelten, 
d.h. erst durch Einführung der Zylinderkrümmung sich 
ergebenden, Gliedern vorkommen, kann man die Diffe- 
rentialgleichungen sogar ohne Bedenken in lineare Diffe- 
rentialgleichungen mit konstanten Koeffi- 
zienten vereinfachen. In der weiteren Rechnung 
wird also in (12*) und (13*) & neben r, vernachlässigt. 
Die Lösung der linearen Differentialgleichungen mit 
konstanten Koeffizienten lautet dann.: 


1 
(- 2r ir ) R 
Er 
ae | (4#sinvz + A%cosv «) 
ni 
(- ET ut) 
+e (4$sinvc + A%cosv.) (16*), 
1 
(- > .) Basar: 
ga) =e (4% sinv x — At cosv «) 


n =): 
( 


— 4% sinvc + A*cosv «) (17, 


worin 
3 =%°p und vy= =) . (18%) 
p 
ist mit 
„= ET und VI Fat a, 
mit \ (19*) 
G 
— 1 N 
2 ( 4,2 Or? 02 | 
und 


nn 
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Die Konstanten A# bis Af ergeben sich aus den Rand- 
bedingungen. Man erhält das Gleichungssystem: 


(Ar +45) + (4% — AP) p? 


-41- 
v 


(42 — A) — (Ar +43) p? 


1,5 &) (At HAN. . (21*), 
rg 0 


-3(- 1) Ar) 2.22%) 
au are 2) -arlunrs 
0 
15 
+ [are + öl) + Ar|coro 
0 
— \ .@3*). 
men, 3106 
= *(o en + 4% |sin»ö 


Die 4. Gl. (24*), die aus (15*) sich ergibt, wird hier 
nicht angegeben, da sie zur Rechnung der ‚„Phasenver- 
schiebung 90°-Methode“ nicht gebraucht wird. 

Für die Phasenverschiebung 90° erhält man das Er- 


gebnis 
DW ya 
9 - >) 
P p 
s 1 1,5 (20,6 — 1,5) p2 
= oö°z oe 1 — - 
@ - Por, p) | vor, r3 (1+ q#) 
x Tg (dog. dp) 2)" 
Unter den in der Hauptarbeit angegebenen Bedin- 
gungen werden dann die GIn. (28a), (29a) und (30a) zu 


1,5 (z or, ö — 1,5) 9° 
Yorr 19 (l om p%) 


R HILFE... 
Mrz >. Bir 5 
2 1—-V), 
EEERTN | (282%), 
Ti 4,256, | 
0,302 | 
| 
OO OE\ 3 
k=0)] € a 5 (292 *), 
N 
De 29b*) 
nz (290*), 
ou? a? 
n=— 7 (1+04E+0,07R +...) 
N 
6 23,25 v 
— —— E en. *), 
(1 : ee )a ae j (30a *) 
Kassel 1. E.Pieruschka. 


Ein Näherungsverfahren für die Integration 
der laminaren, kompressiblen Grenzschicht- 
gleichungen. 

Wir haben die kürzlich erschienene Arbeit von W. 
Dienemann [1] mit Interesse gelesen und möchten 
hier über ein ähnliches von den Verfassern entwickeltes 
Berechnungsverfahren [2] berichten. Wie es bei dem 
Pohlhausenverfahren üblich ist, beginnt man mit 
Integration der x-Impuls- und Energiegrenzschicht- 
gleichungen in bezug auf die senkrechte Koordinate %, 
zwischen den Grenzen y = 0 und y = öy oder Ay, wobei 
öy und 4, die entsprechenden Grenzschichtdicken sind. 
Es ist bekannt (vgl. z. B. [3]), daß die integrierten 
Gleichungen eine besonders einfache Form annehmen 
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durch Einsetzen der Dorodnitsynschen Größe N 
welche durch die folgende Transformation bestimmt ist, 


m =3,|, „m, 
(0) =0, (Ay = A 


Hierbei bedeutet o die örtliche Dichte. Der Index ot 
bezieht sich auf die Staubedingungen der ungestörten 
Strömung. Für eine konstante Prandtlsche Zahl 
werden die integrierten Gleichungen!) 


dd d Aw .7; 
ru A ee @) 
d® d Bi q(&) \ 
dx EZ , Zr us Cp Oot Tot U 
wobei 
ou 
To f ah 
oT\ 
«= -(25,), 
ö 
u 1 
9-7 zan. 
0 
ö 
u ee) 
= (1 Ber a) dy, 
6 U 
A 


und T; die lokale Stautemperatur sind. Um Gl. (2) zu 
integrieren, müssen die sechs, in Gl. (3) hervorgehobe- 
nen Größen auf zwei reduziert werden. Dies wird beim 
Pohlhausenverfahren durch passende Darstellung 
der Grenzschichtprofile erreicht, wobei die betrachteten 
sechs Größen durch zwei Größen ausgedrückt werden, 
im wesentlichen ö und J. Es möge in diesem Zusam- 
menhang hervorgehoben werden, daß die Profildar- 
stellungen die Berechnung von Hautreibung und 
Wärmeübergang nur insofern berühren, als sie die 
erforderlichen Beziehungen für Reduktion der Vari- 
ablen liefern. 

Die gegenwärtige Theorie nimmt Geschwindigkeits- 
und Temperaturprofildarstellungen in der folgenden 
Form an, 


BEER) 


T — To 9 7 N 
T, —/too e ) ne Atı 9, en + Ata 9a ee J 


wobei /, dem gewöhnlichen Pohlhausen-Formpara- 
meter entspricht und A, dem Formparameter A; von 
Dienemann analog ist; A;, ist ein zusätzlicher 


Kompressibilitätsfaktor und /;, bedeutet die bekannte 
Funktion (To — Tx)/T&- Mit Hilfe der Grenzbe- 
dingungen an der Wand, welche man von den origina- 
len (nicht integrierten) Differentialgleichungen erhält, 
kann man alle unbekannten Formparameter als Funk- 
tionen der abhängigen Variablen ö und A ausdrücken. 
Diese, etwas langwierige Einzelrechnung soll hier nicht 
wiedergegeben werden. 

Der nächste Schritt, bevor man Gl. (2) integrieren 
kann, ersetzt die Größen 7,, 9, ö,, 41, q (x) und © durch 
die Formparameter und durch das Verhältnis 4/6. 
Damit wird die Zahl der abhängigen Größen auf zwei 
reduziert. Für die Integration des gegebenen Glei- 
chungssystems ist es relativ gleichgültig, ob die Größen 


ı) Es werden Dienemanns Symbole benutzt. In Ausnahme- 
fällen wird das verwendete Symbol ausdrücklich erklärt. 


Schließlich wird noch ein System von zwei Gleichun- 
_ gen benötigt, um die Anfangswerte der abhängigen Varia- 
'blen an der Plattenvorderkante zu bestimmen. Die erste 


ab ;e Variable gewählt w 
alle kannten Größen durch das ge 
gedrückt werden können. te 


Gleichung kann man leicht auf die gewöhnliche Art und 
Weise erhalten, denn es muß sein A,(0) =0 oder 
(dd/da)a=0# ®, je nachdem ob U(0)#0 oder 
U(0)=0. Die zweite Gleichung erhält man unter der 
Annahme, daß U 2 und U ®2an der Vorderkante ver- 
schwinden, was in jedem Falle gültig sein soll. Man 
‚erhält dann, unter Anwendung des Ta ylorschen 
Theorems, die Identität ; 
d d N Fr 
(Sr | a | 2) | 92 > ( 

ie u ) g 4 Sr | r z=0 m 
welche eine zusätzliche Beziehung zwischen den ab- 
hängigen Variablen liefert. Für den Grenzfall 7W=T% 
und M% (Machzahl der ungestörten Strömung) = 0, 
wird Gl. (5) unbestimmt. Aus der Annahme kontinuier- 
licher Anfangswerte in der Nachbarschaft dieses Singu- 
laritätspunktes ergibt sich eine neue Bedingung, die 
durch die Profildarstellung befriedigt werden muß und 
die für die gegenwärtige Theorie von Bedeutung ist. 
Der Bedingung wird dadurch Rechnung getragen, daß 
wir die Profile, (n/4) und 9, (n/4) in der folgenden 
Form ausdrücken, 


In (nd) = nı (mA) + Agn2 (n/A) | 
RG (2) re J 


wobei A als Eigenwert von Gl. (5) behandelt wird. 

Die Rechnung wurde zunächst für die ebene Platte 
durchgeführt, unter Annahme der korrekten Zähigkeit- 
Temperatur-Beziehung vonSutherland. Es wurde 
jedoch dabei festgestellt, daß die auf diese Weise ein- 
geführten nichtlinearen Glieder zu einer erheblichen 
Abweichung der Lösung von den korrekten Werten [4] 
führten. Dies deutet an, daß ein nichtlineares Zähig- 
keitsgesetz für eine Pohlhausen- Näherungsmethode 
nicht brauchbar ist. Die genaue Ursache hierfür wurde 
noch nicht festgestellt. Jedenfalls wird die Schwierig- 
keit durch Annahme eines linearen Zähigkeitsgesetzes 
vermieden, was gleichzeitig zu einer erheblichen Ver- 
einfachung ‘des Formalismus führt. Zur Darstellung 
der Geschwindigkeitsprofile wurden Polynome vierten 
wie auch sechsten Grades benutzt. Für die Temperatur- 
profile wurden ausschließlich Polynome sechsten Grades 
benutzt. 

Trotz ziemlich umständlicher Symbolik wird die 
Integration des allgemeinen Falles der Gl. (2) [welche 
die Berechnung der Schubspannung 7, und des Wärme- 
zustromes q (=) enthält] ohne besondere Schwierigkeiten 
durchgeführt und zwar mit Hilfe einer elektronischen 
Digitalrechenmaschine mit IBM-Kartenprogramm, bei 
der die allgemeine Verschlüsselung nur einmal am An- 
fang zu erfolgen braucht. Obwohl im allgemeinsten 
Fall keine Standardlösung existiert, so zeigt doch ein 
Vergleich der hier vorgebrachten Theorie mit genauen 
Lösungen von weniger allgemeinen Fällen, daß die 
erhaltene Genauigkeit für technische Zwecke hinrei- 
chend ist. 


Buffalo, New York. D. N. Morris 
J. W. Smith. 
Schrifttum 
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Br, Ram 
7 0,9827 « + 0,3110 b + 0,2867 Ee. 
Die Abweichungen, die die Ergebnisse jeder Näherungs- 
gleichung von den genauen Werten zeigen, sind für ver- 
schiedene Werte des Achsenverhältnisses b/a ver- 
schieden. Die Zahlentafel zeigt die Abweichungen A 
der drei oben angeführten Näherungsgleichungen in 
Prozenten des genauen Wertes in Abhängigkeit vom 
Achsenverhältnis bj/a. > 
Wenn man die maximalen Abweichungen der Zahlen- 
tafel betrachtet, sieht man, daß die Näherungs- 
gleichungen in derselben Reihenfolge an Genauigkeit 
zunehmen, in der sie an Einfachheit ihres Aufbaues 
verlieren. 
Im Grenzfalle 
Sue 
[77 Tr 
in dem die Ellipse in einen Kreis vom Radius a über- 
geht, wird der Grenzwert 


IN 2 

ET 
von den Gleichungen (1) und (2) genau erfüllt, im PR 
Grenzfalle 

a 

@ 


in dem die Ellipse in zwei zusammenfallende Strecken 
übergeht, wird der Grenzwert 
Dat 5 
+4 
von keiner der drei Gleichungen genau erfüllt. 

Es läßt sich nuneine verbesserte Näherungsgleichung 
angeben, die mindestens die Gleichungen (2) und (3) an 
Einfachheit übertrifft, aber allen drei Gleichungen an 
Genauigkeit überlegen ist. Sie lautet: 


U _nab +(a —b)® 
er ab 


Ihre Abweichungen sind in einer weiteren Spalte der 
Zahlentafel enthalten. Sie erfüllt die beiden Grenz- 
werte genau, ihre maximale Abweichung liegt erheblich 
unter den maximalen Abweichungen der anderen drei 
Näherungsgleichungen. Sie sollte daher in Zukunft an 
deren Stelle Verwendung finden. i 
Da die Abweichungen der neuen Gleichung nur auf 

der positiven Seite der Null-Linie liegen, kann man 
die maximale Abweichung noch auf die Hälfte herab- 
setzen, wenn man die Ergebnisse aus der Gleichung 
durch Multiplizieren mit dem Faktor 0,997 um soviel 
verkleinert, daß auf der positiven und negativen Seite 
etwa gleich große Abweichungen auftreten. Man be- 
nutzt dann die korrigierte Gleichung: 

172% nab + (a — b)? 

5 0,997 ne 2 (9); 


*) „Hütte‘‘, d. Ingenieurs Taschenbuch, Bd.I, 26. Aufl. Berlin 
1931, 8.119. 


. (4). 


Z. angew. Math. Mech. 
Bd. 34 Nr. 4/5 April/Mai1954 


Zahlentafel 


Achsen- 


verhältnis | ae khune a = D% , 100 der Näherungsgl. Nr. 
ea RR, a oe ),E 
1 0 1) +0,61 0 I 9,30 0 
0,7071 |-+ 0,01 0,00 —0,35 |+0,07 |—0,23 | —0,02 
0,5000 |-+ 0,06 1—0,09 —0,10 |+0,23 |—0,07 | —0,05 
0,2588 |+ 0,66 |—0,78 |+0,56 |+0,55 I|+0,25 | —0,07 
0,1736 |-+ 1,46 |—1,50 +0,50 +0,62 |+0,33 | —0,02 
0,0872 |+ 3,57 —2,79 —0,07 |+0,55 | +0,26 | +0,05 
0,0523 I+ 5,45 |—3,57 |—0,56 140,43 |-+0,13 | +0,07 
0,0175 |+ 9,40 |—4,58 —1,16 +0,19 |—0,25 | +0,05 
0 +17,81 |-919 —1,73| 0 —0,30 0 
Maxim. 
Abweich.|+17,81 )— 5,19 |—1,73 +0,62 | +0,33 | +0,07 


deren Abweichungen in einer weiteren Spalte der 
Zahlentafel eingetragen sind. 

Die Gleichung (4) wurde aus einer Reihenentwick- 
lung für die genaue Berechnung des Ellipsenumfanges 
hergeleitet, die an der gleichen Stelle der ‚Hütte‘ mit- 
geteilt wird. Sie lautet: 


Un ee 
tl) 


und setzt sich mit wachsenden geraden Potenzen von 


= nb t fort 
er unbegrenzt fort. 

Brieht man die unendliche Reihe bereits hinter dem 
zweiten Glied mit der zweiten Potenz ab, so muß man 
den Koeffizienten 1/4 dieses Gliedes durch einen etwas 
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höheren Zahlenwert ersetzen, um den Wegfall der 
übrigen Summenglieder auszugleichen. Bestimmt man 
diesen Koeffizienten aus der Bedingung, daß die 
Gleichung im Grenzfall 


den Grenzwert 


erfüllt, so nimmt der Koeffizient den Wert 


et, 
ı 
an. Damit ergibt sich die Gleichung (4). 

Man kann natürlich eine noch genauere Näherungs- 
gleichung gewinnen, wenn man die unendliche Reihe 
in obiger Gleichung erst nach dem dritten Glied ab- 
bricht und die beiden Koeffizienten der Potenzglieder 
so bestimmt, daß die Abweichungen vom genauen 
Hr möglichst klein werden. Eine solche Gleichung 
autet: 


U _ zab[2 (a? + b?) + 0,244 (a —b)?] + (a — b)? 
4 (a + 5)? Ba gene 


Ihre Abweichungen sind in der letzten Spalte der 
Zahlentafel enthalten. Diese Näherungsgleichung ist 
jedoch bereits so kompliziert geworden, daß sie ihren 
praktischen Wert einbüßt, während man andererseits 
ihre maximale Abweichung noch nicht auf 1/10 des 
Wertes der einfachen Gleichung (4) herabdrücken kann. 
Man wird daher praktisch wohl stets die Gleichung (4) 
oder (5) verwenden. 


Göttingen. 


(6). 


W.Weydanz. 


BUCHBESPRECHUNGEN 


Dr. Alfred Müller, Nomographie für die 
teehnische Praxis 2688. .m. 97 Abk. 
Leipzig 1952. Fachbuchverlag. Preis geb. 9,80 DM. 

Der erste Teil dieses Buches gibt eine ganz _ele- 
mentare, wenig mathematische Kenntnisse erfor- 
dernde, mit großem didaktischen Geschick geschriebene 
Einführung in die Nomographie. Wertvoll sind vor 
allem die im zweiten Teil zusammengefaßten Nomo- 
gramme aus den Gebieten der Maschinenelemente, der 
Werkzeugmaschinen, der spanlosen Formung, der 
Wärmewirtschaft, der Elektrotechnik, der Mechanik 
und Festigkeitslehre, die von auf diesen Gebieten 
arbeitenden Ingenieuren beigesteuert sind. Das Buch 
kann sowohl den in der Praxis stehenden Technikern 
eine Hilfe sein, wie auch einführenden Kursen zu- 
grundegelegt werden. 

Dresden. Willers 

Jahrbuch 1952 der Wissenschaftlichen 
Gesellschaft für Luftfahrt e.V. (WG). 
Herausgegeben von H..Blenk. 1488. mit zahl- 
reichen Abb. Braunschweig 1953. Friedr. Vieweg 
& Sohn. Preis geb. 24,— DM. 

Die Wissenschaftliche Gesellschaft für Luftfahrt, 
die 1952 neu gegründet wurde, sieht die vom 21. bis 
23. April in Braunschweig von der Technischen Hoch- 
schule, dem Bezirksverein des VDI und der Orts- 
gruppe der Arbeitsgemeinschaft für Luftfahrtforschung 
abgehaltene Tagung als WGL-Tagung an. In dem 
vorliegenden Jahrbuch werden nach einleitenden Auf- 
sätzen über die Geschichte der WGL, ihre Satzungen, 
ihrem Mitgliederbestand usw. die Vorträge der 
Braunschweiger Tagung wiedergegeben, die in über- 
sichtlicher Zusammenfassung über die Entwicklung 
der Luftfahrt und der Luftfahrtforschung seit Kriegs- 
ende berichten. 


Dresden, Willers 


Dr.-Ing. .W. Meyer zur Capellen, Leitfaden 
der Nomographie. IV +178S. m. 203 Abb. 
Berlin-Göttingen-Heidelberg 1953. Springer-Verlag 
Preis geb. 17,40 DM. 

Im ersten Teil des Buches werden die theoretischen 
Grundlagen besprochen, wobei sich schon hier wich- 
tige Hinweise für die Praxis finden, u.a. werden 
Genauigkeitsfragen bei verschiedenen Nomogramm- 
formen behandelt sowie die projektive Abbildung 
und Dualität, deren Kenntnis für die Beherrschung 
der Bereiche unerläßlich ist. 

Der zweite Teil, der fast die Hälfte des Buches aus- 
füllt, enthält Beispiele aus der Praxis, deren Ansatz 
erläutert wird. Dort, wo dieser auf eine zunächst 
noch unbefriedigende äußere Nomogrammform führt, 
weist der Autor auf Verbesserungen der Bereichs- 
darstellung hin, wobei die projektive Umformung 
wegen ihrer Schmiegsamkeit bevorzugt wird. Die 
exakten Zeichnungen verdienen besondere Erwäh- 
nung. Der Name des Bahnbrechers der Nomographie 
in Frankreich, d’Ocagne, ist im Buche leider durch- 
weg mit zwei c geschrieben. 

Das Buch, das einige mathematische Kenntnisse 
voraussetzt, wendet sich zwar in erster Linie an 
Ingenieure, Mathematiker und Naturwissenschaftler, 
es wird aber nicht nur diesen, sondern auch allen 
denen eine wertvolle Hilfe sein, die beruflich mit sich 
ständig wiederholenden Berechnungen zu tun haben. 


Berlin-Tempelhof. Joh. Fischer. 


W.J. Duncan D. Se., F.R.S. (Mech. Prof. in the 
University of Glasgow), The Principles of 
the Controlvand=Stabilıty oR Ayr- 
eraft. (Cambridge aeronautical, Series I). XVI-+ 
384 S. Cambridge 1952. University Press. Preis 40. 

Mit dem vorliegenden Buch wird das Erscheinen 
einer vom College of Aeronauties in Cranfield ge- 
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planten Reihe von Spezialwerken auf luftfahrttech- 
nischem Gebiet eingeleitet. Es soll damit dem auf 
diesem Gebiet bestehenden Mangel an lehrbuchartigen 
Darstellungen, die die neuesten Erkenntnisse berück- 
sichtigen und dem durch die Entwicklung der Luft- 
fahrt bedingten erweiterten Problemkreis Rechnung 
tragen, abgeholfen werden. Durch die Beschänkung 
auf Teilgebiete, wie hier auf die besonders wichtigen 
und komplizierten Fragen der Steuerung und Stabili- 
tät von Flugzeugen, ist die Möglichkeit gegeben, diese 
in erschöpfenderer Weise als sonst in den Aero- 
dynamiklehrbüchern üblich zu behandeln und zu- 
gleich die Bearbeitung einem Fachmann auf dem 
speziellen Gebiet zu übertragen (im vorliegenden Buch 
sind darüber hinaus noch zwei Kapitel von Prof. 
A.D. Young übernommen worden). 

Die im Hinblick auf die Fülle des Stoffs auch hier 
erforderliche Auswahl wurde so getroffen, daß der 
Leser mit allen wichtigen Problemen eingehend ver- 
traut gemacht wird. Dazu gehört neben den all- 
gemeinen Stabilitäts- und Steuerungsproblemen eine 
Untersuchung der zahlreichen speziellen Flügel-Ruder- 
Konstruktionen, der Einfluß der elastischen Verfor- 
mungen des Tragwerks sowie der Kompressibilität 
der Luft. Der Verf. hat mathematisch schwierigere 
theoretische Ableitungen im allgemeinen vermieden, 
indem er nur die entsprechenden Ergebnisse anführt 
oder aber den Fragekomplex in großen Zügen be- 
handelt. Gleichzeitig wird auch auf eine übermäßige 
Aufzählung experimenteller Daten trotz deren Be- 
deutung für das vorliegende Gebiet verzichtet, sondern 
vielmehr versucht, diese in den Rahmen einer mög- 
lichst elementaren Theorie einzubauen. 

Die Absicht des Verf., dem angehenden Flugzeug- 
aerodynamiker ein entsprechendes Lehrbuch in die 
Hand zu geben, kann als voll gelungen bezeichnet 
werden. Zugleich ist dieses Buch aber auch allgemein 
für den Aerodynamiker zur Orientierung über die 
modernen Probleme auf diesem Gebiet von Nutzen. 


Dresden. K. Krienes. 


H.Drescher, Untersuchungenaneinem 
symmetrischen Tragflügel mit spalt- 
los angeschlossenem Ruder bei ra- 
schen Anderungen des Ruderaus- 
schlags (ebene Strömung). (Mitteilungen 
aus dem Max-Planck-Institut für Strömungsfor- 
schung. Nr.6.) I-+ 71S. mit 45 Abb. u. 2 Tab. 
Göttingen 1952. Selbstverlag Max-Planck-Institut 
für Strömungsforschung. Preis brosch. 9,— DM. 

Von den umfangreichen Arbeiten aus dem Krieg 
ruht immer noch ein sehr großer Teil, der Allgemein- 
heit entzogen und vergraben. Da ist es besonders zu 
begrüßen, daß die Mitteilungen des Max-Planck- 
Instituts für Strömungsforschung sich auch in den 
Dienst der Hebung solcher Schätze gestellt haben. 
Die vorliegende Mitteilung bringt reiches experi- 
mentelles Material zu den Strömungserscheinungen 
am Tragflügel beim Übergang in einen anderen sta- 
tionären Zustand. Wenn die veröffentlichten Ergeb- 
nisse (24 Versuche) auch nur einen kleinen Teil des 
ganzen Materials (250 Versuche) darstellen, so geben 
sie doch einen guten Überblick über die Erscheinungen 
durch Druckmessungen an 21 Meßbohrungen und 
Photographien von Farbspuren, auch mit Film- 
kamera und Zeitlupe. Die Hauptergebnisse der mit 
Wasser durchgeführten Versuche sind folgende: 
l. Gruppe mit Anstellw. « = — 5°, Ruderwinkel 
ß = 0° > 15°, liefert durchweg anliegende Strö- 
mung; 2. Guppe mit =0° und —5°, = 
0°—>60° am Ruder abreißend; 3. Gruppe 
& = — 15° bis +12,5°, 6 = 0° — 60° an Ruder und 
Hauptflügel abreißend. Bei anliegender Strömung 
(kleine Winkel) zeigte sich gute Übereinstimmung mit 
der instationären Theorie für die geknickte ebene 
Platte. 


Weilburg. L. Schiller. 
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Francesco Giacomo Trieomi, Equazioni dif- 
ferenziali. 2. erweiterte Aufl. 353 S. m. 26 Abb. 
Torino 1953. Edizioni Seientifiche Einaudi. Preis 
geb. 4000 L. i 

Der Autor stellt sich die Aufgabe, einige Gebiete aus 
der Theorie der gewöhnlichen Differentialgleichungen, 
die in den üblichen Anfängervorlesungen gewöhnlich 
zu kurz kommen, mit der wünschenswerten Ausführ- 
lichkeit zu behandeln, wobei bis zu neueren For- 
schungsergebnissen, an denen der Verf. maßgebend be- 


teiligt ist, vorgestoßen wird. Unter bewußtem Ver- 


zicht auf elementare Dinge (‚integrable Typen“ usw.), 
die in allen guten Lehrbüchern der höheren Analysis 
nachgelesen werden können, beginnt das Buch in 
seinen ersten (einleitenden) Kapiteln mit dem fun- 
damentalen Existenzbeweis von Picard-Linde- 
1öf für Systeme von Differentialgleichungen erster 
Ordnung. Das zweite Kapitel befaßt sich mit 
geometrischen Eigenschaften der Lösungskurven, vor 
allem mit dem Verhalten in der Nachbarschaft von 
singulären Punkten und mit periodischen Lösungen. 
Im dritten und vierten Kapitel bebandelt der Verf. 
die linearen Differentialgleichungen, wobei er sich aus 
Raumgründen im allgemeinen auf diejenigen zweiter 
Ordnung beschränkt. Es finden sich dort insbesondere 
das Oszillationstheorem (unter Benutzung der auf 
H. Prüfer zurückgehenden Variablentransfor- 
mation), die Eigenwerttheorie, die Greensche 
Funktion und die Zusammenhänge mit linearen Inte- 
gralgleichungen, Betrachtungen über das asympto- 
tische Verhalten der Lösungen für mw (z B. 
Stabilitätsprobleme) und asymptotische Formeln für 
Eigenwerte und Eigenfunktionen für große Werte des 
Parameters nebst Anwendungen auf Laguerre- 
sche und Legendresche Polynome. Nach diesen, 
der reellen Theorie gewidmeten Abschnitten befaßt 
sich das fünfte und letzte Kapitel mit Untersuchungen 
im Komplexen. Neben dem fundamentalen Existenz- 
satz gelingt hier durch Verwendung funktionen- 
theoretischer Hilfsmittel eine abgerundete Betrach- 
tung der linearen Differentialgleichungen zweiter 
Ordnung der Fuchsschen Klasse, aber auch das 
Studium wesentlich singulärer Stellen nebst An- 
wendungen auf die konfluenten hypergeometrischen 
Funktionen und die Zylinderfunktionen. 

Das Buch wird jedem Leser, der über die Elemente 
der höheren Analysis und (für Kapitel V) der Funk- 
tionentheorie verfügt, einen tiefgehenden und durch 
seine klare und didaktisch geschickte Darstellung gut 
verständlichen Einblick in die genannten Gebiete 
gewähren. 


Dresden. Maruhn. 


A. Timpe (em. o. Prof. a. d. TU Berlin), Ein- 
führungkin-dierfinanz-undsWor:- 
schaftsmathematik. 2. verb. Aufl. VI+ 
217 S. m. 70 Abb. Wiesbaden 1953. Verlag für an- 
gewandte Wissenschaften. Preis geb. 17,80 DM. 

Die 1934 erschienene erste Auflage des vorliegenden 
Buches wurde in dieser Zeitschrift Bd.15 (1935), S. 181 
vonv.Mises eingehend besprochen. Es war längere 
Jahre vergriffen und erscheint jetzt in einer nicht 
wesentlich veränderten Form. In ihm wird ein um- 
fangreicherer Stoff behandelt, als man ihn üblicher- 
weise in derartigen Büchern findet. Dieser wird in 
acht Abschnitten dargeboten, die die Überschriften 
tragen: 1. Einfacher Geldverkehr, Wertpapiere, 
2. Zinseszins- und Rentenrechnung. Versicherungs- 
rechnung, 3. Tilgungs- und Abschreibungswesen, 
4. Kurs- und Rentabilität von Tilgungsanleihen, 
5. Tarifmathematik, 6. Erfolgsrechnung, 7. Wahr- 
scheinlichkeitsrechnung und mathematische Statistik 
in Anwendung auf Wirtschaftsprobleme, 8. Die Analyse 
von Zeitreihen. Neubearbeitet ist der Paragraph über 
Rentenrecbnung. Auch sonst sind Änderungen und 
Ergänzungen hinzugefügt. Gern gesehen hätte man, 
daß die erst in den letzten Jahren in Deutschland be- 
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kanntgewordenen, hierher gehörenden Arbeiten über 


die Anwendung der Methoden der mathematischen 


Statistik berücksichtigt wären. Aber auch in der vor- 
liegenden Form wird das klar geschriebene Buch, das 
an Mathematikkenntnissen nur die durch eine gute 
Schule vermittelten voraussetzt, vielen eine sehr 
brauchbare Einführung in die behandelten Gebiete 
sein. i 


Dresden. Willers. 


I. 6. Petrovskij, Vorlesungen über die 
Theorie der Integralgleichungen. 
Übersetzt von R. Herschel. 100 S. Würzburg 1953. 
Physica-Verlag. Preis 7.80 DM. 

Die vorliegende Einführung in die Theorie der 
linearen Integralgleichungen zweiter Art entspricht 
etwa dem, was in einer Vorlesung für Anfänger auf 
diesem Gebiete gebracht zu werden pflegt. Das erste 
Kapitel befaßt sich mit den Fredholmschen 
Sätzen, soweit sie determinantenfrei bewiesen werden 
können. Es werden also ausgeartete Kerne, die 
Neumannsche Reihe, die fastausgearteten Kerne 
und die (durch ausgeartete Kerne gleichmäßig 
approximierten) gleichmäßig stetigen Kerne behan- 
delt. Ausführlich wird auf diejenigen nichtbeschränk- 
ten Kerne eingegangen, auf diesich deFredholm- 
schen Sätze übertragen lassen. Nach einem kurzen 
zweiten Kapitel über Volterrasche Integral- 
gleichungen ist das dritte Kapitel den Gleichungen 
mit reellem symmetrischen Kern gewidmet. Es 
bringt eine für den Anfänger recht anschauliche Ein- 
führung des Funktionenraumes, die wichtigsten Sätze 
über Eigenwerte und Eigenfunktionen (insbesondere 
die Existenz eines Eigenwertes nach Gelfand) 
und die verschiedenen Entwicklungssätze. Allen 
diesen Betrachtungen liegt wie üblich der Rie- 
mannsche Integralbegriff zu Grunde. Ein Aus- 
blick auf Kerne, die im Lebesgueschen Sinne 
quadratisch integrabel sind, bildet den Abschluß des 
Buches. 

Das Buch ist klar und mit didaktischem Geschick 
abgefaßt. Obwohl es nicht wesentlich aus der Reihe 
der vielen Bücher mit gleichem Thema tritt, ist es 
doch jedem, der sich in diese Materie einarbeiten 
will, durchaus zu empfehlen. 


Dresden. Maruhn. 


Dr. Günther Pickert, Analytische Geo- 
metrie. (Mathematik und ihre Anwendungen in 
Physik und Technik, Reihe A, Bd. 24.) X + 398 S. 
mit 76 Abb. Leipzig 1953. Akademische Verlags- 
gesellschaft Geest und Portig KG. Preis 26.— DM. 

Verfasser stellt sich die Aufgabe, die Analytische 
Geometrie von Grund auf, ausgehend von einem 
Axiomensystem aufzubauen, das mit den Grund- 
begriffen ‚Punkt‘ und ‚Vektor‘ arbeitet und die 
reellen Zahlen verwendet. Er berücksichtigt, daß man 
für viele Herleitungen statt aller für die reellen Zahlen 
gültigen Gesetze nur die Eigenschaften eines Schief- 
körpers benötigt (in dem das kommutative Gesetz der 
Multiplikation nicht gilt). Ganz allgemein legt Ver- 
fasser besonderen Wert darauf, den Studierenden von 
vornherein an die Begriffsbildungen und Gedanken- 
gänge der modernen Mathematik heranzuführen, und 
verzichtet dafür lieber auf eine gewisse Breite des 
geometrischen Stoffes. Rund 270 Aufgaben dienen 
zur Einübung des Vorgetragenen und zu dessen Er- 
gänzung. Ein Literaturverzeichnis enthält Bücher 
und Abhandlungen, die den Leser über das in dem 
Werk Gebotene hinausführen können. 

Das Buch enthält drei Hauptteile: I. Affine 
Geometrie. Ein n-dimensionaler reeller affiner 
Raum A(r) wird als eine Menge von zwei Arten von 
Dingen erklärt, die als Punkte und Vektoren bezeich- 
net werden und für welche eine Reihe von’ Axiomen 
gelten. Nach Definition der Begriffe ‚‚Vektorschar‘ 
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und ‚linearer Unterraum‘‘ wird die lineare Algebra 


des affinen Vektorraumes abgeleitet, wobei von 
Matrizen und Determinanten Gebrauch gemacht wird. 

II. Metrische Geometrie. Es werden 
zwei neue Grundbegriffe, die „Länge eines Vektors‘“ 
und das „Senkrechtsein von Vektoren‘ und einige 
neue Axiome eingeführt. Ein A{r) mit n > 2, in wel- 
chem die zusätzlichen Axiomen gelten, ist ein n- 
dimensionaler metrischer Raum M(r). In diesem wird 
eine symmetrische Fundamentalform zweiten Grades 
eingeführt. Die weitere Entwicklung der metrischen 
Geometrie geht über die Begriffe Pythagoräischer 
Raum, Winkel zweier Vektoren, archimedische An- 
ordnung, Bewegung, Inhalt und orientierter Inhalt 
eines Simplex zu den Hyperflächen zweiter Ordnung. 

II. Projektive Geometrie. Der n- 
dimensionale projektive Raum //® wird aus einem 
affinen n-dimensionalen Raum durch Hinzunahme 
der uneigentlichen linearen Unterräume gewonnen. 
Als Kollineationen zweier projektiven Räume werden 
die eindeutigen Punktabbildungen definiert, bei wel- 
chen drei Punkte einer Geraden stets wieder drei 
Punkten einer Geraden entsprechen. Projektiv heißt 
eine Kollineation, wenn sie die Doppelverhältnisse in- 
variant läßt. Die Projektionen führen zu ausgearteten 
Kollineationen. Über Dualitäten, insbesondere Polari- 
täten und Nullsysteme im 7% gelangt man zu den 
Quadriken. Die affine Geometrie und die Metrik 
werden projektiv erklärt. Durch Auszeichnung einer 
(n — 1)-dimensionalen Hyperfläche zweiter Ordnung 
wird die nichteuklidische Geometrie und die zu- 
gehörige Metrik eingeführt. Die projektive Linien- 
geometrie des //(3) wird mit Hilfe einer Abbildung der 
Geraden des 1/(3) auf Punkte einer besonderen Qua- 
drik des //(5) behandelt. — Ein Anhang bringt 
als Ergänzung zum systematischen Aufbau der analy- 
tischen Geometrie dieEinführung der trigonometrischen 
Funktionen ohne Benützung des Winkelbegriffs, die 
Kennzeichnung der projektiven Ebene ohne Ver- 
wendung von Vektoren und Skalaren und die 
Automorphismen des Körpers der reellen Zahlen. 

Das Werk ist insbesondere solchen Studierenden zu 
empfehlen, die nicht nur die praktischen Methoden 
der analytischen Geometrie kennenlernen wollen, 
nach welchen konkrete geometrische Aufgaben mit 
Hilfe von Rechnungen gelöst werden können, sondern 
denen es vor allem darum zu tun ist, die Grundlagen 
der analytischen Geometrie, ihren begrifflichen Auf- 
bau und die tieferen Zusammenhänge zwischen ihren 
großen Zweigen, die durch die Gedankengänge der 
modernen Mathematik vermittelt werden, kennen- 
zulernen. 


Freiberg (Sa.). W. Schmid. 

Dr. P. Funk (o. Prof. a. d. T. H. Wien), Dr. Hans 
Sagan und Franz Selig, Die Laplace-Trans- 
formation und ihre Anwendungen. 
VII + 106 S. m. 13 Abb. Wien 1953. Verlag Franz 
Deuticke. Preis kart. 10.— DM. 

Die Laplace-Transformation wird heute schon viel- 
fach, insbesondere bei Behandlung von Einschalt- 
vorgängen verwendet. Auf der Hochschule ist jedoch 
meist keine Zeit, darauf in den mathematischen Vor- 
lesungen einzugehen. Da will das vorliegende Bänd- 
chen einspringen, dem Vorlesungen des ersten Ver- 
fassers zu Grunde liegen. Es ist in der Hauptsache 
für Physiker und Ingenieure bestimmt und verzichtet 
daher auf mathematische Einzelheiten, wie etwa auf 
die genaue Diskussion der Gültigkeitsbereiche der 
Operationen und Korrespondenzen, verweist vielmehr 
da auf entsprechende Lehrbücher. Nur das Heaviside- 
sche Theorem über asymptotische Entwicklung wird, 
wenn auch skizzenhaft, doch ausführlicher dargestellt, 
um zu zeigen, daß dieser die rein mathematischen 
Schwierigkeiten doch beträchtlich unterschätzt hat. 
Um eine Übersicht über Inhalt und Umfang des Bänd- 
chens zu geben, seien hier die Überschriften der einzel- 
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nen Abschnitte angeführt: Lineare Funktionaltrans- 


formationen, Erklärung und Eigenschaft der Laplace- 
Transformation. Die Abbildung der fundamentalen 
Funktionen. Der Heavisidsche Entwicklungssatz. 
Systeme von Differentialgleichungen mit konstanten 
Koeffizienten. Faltung im Objektraum, der Trans- 
formator. Die Abelsche Integralgleichung. An- 
wendung der Laplace-Transformation auf die Wärme- 
leitungsgleichung. Thomson-Kabel. Das Fouriersche 
Integraltheorem. Die allgemeine Umkehrformel. 
Asymptotische Entwicklungen. Übergangsfunktion 
und Frequenzgang. Regelvorgänge. Eine kurze 
historische Übersicht bildet den Abschluß. 

Das Bändchen setzt nur die von der Hochschule 
vermittelten mathematischen Kenntnisse voraus und 
ist durchaus verständlich geschrieben. Es sei daher 
nieht nur den Studierenden, sondern auch den in der 
Praxis stehenden Elektroingenieuren und Physikern, 
die sich mit Theorie und Anwendung der Laplace- 
Transformation vertraut machen wollen, angelegent- 
lichst empfohlen. 


Dresden. Willers. 


Prof. Dr. Friedrich Böhm, Versicherungs- 
mathematiklIl. (Sammlung Göschen 917/917a.) 
Zweite vermehrte und verbesserte Auflage. 205 S. 
Berlin 1953. W. de Gruyter & Co. Preis 4.80 DM. 

Der vorliegende Band II des Werkes ergänzt die im 
Band I, 180 der Sammlung, dargestellte elementare 
Versicherungsrechnung nach verschiedenen Seiten hin 
und zerfällt in die beiden Teile Lebensversicherungs- 
mathematik und Einführung in die technischen Grund- 
lagen der Sozialversicherung. Der erste Teil gliedert 
sich in die 8 Kapitel: Verzinsung und Zeitrente, 
Sterblichkeit und Lebenserwartung, Erlebensfall- 
versicherung und Leibrente, Todesfallversicherung, 
Versicherungen mit veränderlichem Betrage, Berech- 
nung der Beiträge und Deckungsrücklagen, Beitrags- 
rückgewähr, Versicherung verbundener Leben und 
bringt die analytische Theorie des Zinses und der 
Lebensversicherung sowie die Behandlung ver- 
wiekelterer Versicherungsformen. In den 3 Kapiteln 
des zweiten Teils werden die technischen Grundlagen 
der Invaliditäts-Versicherungsrechnung, die Berech- 
nung der Barwerte und Anwartschaften sowie die Bei- 
tragsberechnung in der Sozialversicherung dargelegt. 

Die leichtverständliche Darstellungsweise, ein an- 
gemessener Stoff statistischer und technischer Tafeln 
und viele geschickt gewählten Zahlenbeispiele machen 
das Buch zum Selbststudium und zum praktischen 
Gebrauche sehr geeignet; ein reichhaltiges Literatur- 
verzeichnis weist dem Leser den Weg zu tieferem Ein- 
dringen in die Versicherungsmathematik. Das Buch 
kann daher warm empfohlen werden. 


Dresden. P. E. Böhmer. 


Kurt Reidemeister (o. Prof. a. d. Univ. Marburg). 
Topologie der Polyeder und kombi- 
natorischeTopologie der Komplexe. 
(Mathematik und ihre Anwendungen in Physik und 
Technik. Reihe A, Bd. 17.) 2. Aufl. XI + 196 S. m. 
69 Abb. Leipzig 1953. Akademische Verlagsgesell- 
schaft. Geest u. Portig KG. Preis geb. 14.— DM. 

Die vorliegende neue Auflage enthält gegenüber der 
ersten Auflage keine wesentlichen Änderungen. Als 
Polyeder bezeichnet der Verfasser jede Punktmenge 
eines linearen n-dimensionalen Raumes, die sich durch 
Durchschnittbildung und Summenbildung aus end- 
lich vielen Halbräumen und linearen Unterräumen er- 
zeugen läßt. Durch diese Definition ist der elementare 
Charakter der hier entwickelten Theorie gekennzeich- 
net, die nur Inzidenzaxiome und Anordnungsaxiome 
der linearen Geometrie, aber keine Stetigkeitsaxiome, 
benutzt. In Kap.I werden die randtreuen Zer- 
legungen eines Polyeders (in konvexe Raumstücke) 
betrachtet. Invarianten eines Polyeders im Sinne der 
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kombinatorischen Topologie sind dann diejenigen 
Eigenschaften, die allen randtreuen Zerlegungen des 
Polyeders gemeinsam sind. In Kap. II wird u. a. der 
Jordansche Satz auf (n — 1)-dimensionale Poly- 
eder im n-dimensionalen Raum übertragen. Eine 
Polyederaddition und ein Verbindungsprodukt von 
Polyedern wird erklärt. Die endlichen Polyeder bilden 
mit diesen beiden Verknüpfungen einen Ring, dessen 
Struktur sich genau übersehen läßt. In Kap. III sind 


u.a. auch die Beziehungen zwischen Stetigkeits-_ 
topologie und kombinatorischer Topologie beschrieben. 


Es folgen Sätze über isolierte Transformationen. 

Kap. IV behandelt topologische Eigenschaften der 

Polyeder. 
Dresden. M. Landsberg. 


1° 


S. €. Dampier, Kurze Geschichte der 
Wissenschaften inihrer Beziehung zu Philo- 
sophie und Religion. XI + 288S. m. 14 Abb. u. 
9 Tafeln. Zürich 1946. Rascher-Verlag. Preis 15,— DM. 

Verf. faßte seine umfangreiche Geschichte der 
Wissenschaften in einen kürzeren Band zusammen, 
der 1943 erschien. Das vorliegende Buch ist eine 
Übersetzung dieses Bandes. Es wird darin die Ent- 
wicklung der Physik, Chemie, Biologie, Medizin, 
Psychologie und Geographie von ihren ersten Anfän- 
gen aninihrer Beziehung zu Religion und Philosophie, 
die besonders in den früheren Zeiten von Bedeutung 
ist, dargestellt. Der geringe Umfang des Buches be- 
dingt an manchen Stellen eine etwas skizzenhafte 
Darstellung. Die Entwicklung der letzten zehn Jahre 
fehlt natürlich. 

Die Übersetzung des ausgezeichneten Buches be- 
dürfte hier und da einer leichten Überarbeitung. Da- 
bei sollte etwa in Anmerkungen noch mehr auf den 
Beitrag Deutschlands zu dieser Entwicklung einge- 
gangen werden, da manche deutsche Forscher, deren 
Werk etwas ausführlicher berücksichtigt werden sollte, 
nur kurz erwähnt oder überhaupt nicht genannt wer- 
den. Konstanten sollten in dem auf dem Kontinent 
üblichen Maßsystem angegeben werden. 


Dresden. Willers. 


Prof. Dr.-Ing. habil. Ulrich Graf und Dr. phil. Hans- 
Joachim Henning. Formeln und Tabellen 
derMathematischen Statistik VI+ 
+ 1028. m. 9 Abb. und 12 Kurvenblättern. Berlin/ 
Göttingen/Heidelberg 1953. Springer-Verlag. Preis 
brosch. 9,— DM. 

Der Benutzer einer mathematischen Formelsamm- 
lung üblicher Art hat sich die hierfür nötigen Grund- 
lagen während seiner Ausbildung erworben. Die Mehr- 
zahl der Ingenieure, die mit vorliegender statistischen 
Formelsammlung arbeiten, wird sich nicht in dieser 
Lage befinden. Darin liegt natürlich eine gewisse 
Gefahr, denn nichts ist für die Verbreitung statisti- 
scher Methoden schädlicher als folgenschwere Fehl- 
schlüsse auf Grund von unsachgemäßer Anwendung. 
Für diejenigen jedoch, die sich die für ihre Arbeit nöti- 
gen statistischen Grundkenntnisse im Selbststudium 
erworben haben, ist das vorliegende Bändchen sehr zu 
begrüßen. Es enthält in übersichtlicher Anordnung 
alles Wesentliche, was der Praktiker für die tägliche 
Arbeit braucht. 

Auf den ersten 30 Seiten findet man eine Formel- 
sammlung, die inihrem Umfang etwa dem entspricht, 
was in dem Buch ‚,‚Statistische-Methoden bei textilen 
Untersuchungen“ der gleichen Verfasser gebracht 
wird. Auch die Tabellen und Kurvenblätter sind 
dem erwähnten Werk entnommen (Normalvertei- 
lung, Integralgrenzen der t-, F- und x?-Verteilung 
zu 95%, 99% und 99,9% Sicherheit, Hilfstafeln 
für Kontrollkarten, Existenzprüfung von Korrelations- 
koeffizienten usw.) und wurden noch um zwei Tafeln 


für 1g N! (bis N —= 999) und für n? und Vr (bis 
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er: 


rn — 2000) ergänzt. Bei Kurvenblatt H fehlt leider 
eine Genauigkeitsangabe. Von den 20 Beispielen, 


die den Gebrauch der Formeln und Tabellen er- 


läutern sollen, sei besonders auf Beispiel 17 (Prüfung 
von Ereignishäufigkeiten im Binomialpapier nach 
Mosteller-Tuke y) hingewiesen; dieses graphi- 
sche Verfahren, das für beliebige Stichprobenumfänge 
verwendbar ist, verdient wegen seiner Bequemlich- 
keit in der Handhabung große Beachtung. Dem Cha- 
rakter einer Formelsammlung entsprechend, wird auf 
theoretische Ausführungen völlig verzichtet. Dafür 
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sind im Literaturverzeichnis knapp 40 Bücher von 
a unterschiedlichem Schwierigkeitsgrad aufge- 
ührt. 

Weil das Bändchen im Gegensatz zu den meisten 
deutschsprachigen statistischen Büchern nicht nur 
für ein begrenztes Anwendungsgebiet geschrieben ist, 
ist es geeignet, der weitesten Verbreitung statistischer 
Methoden zu dienen, und es wird dieser Aufgabe durch 


geschickte Stoffauswahl und Übersichtlichkeit in 


besonderem Maße gerecht. 


Dresden. G.. O’putizi 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten). 


Richard von Mises (Prof. a. d. Harvard Universität 
Cambridge (Mass.)) und Theodore von Kärmän (Prof. 
a. d. Columbia Universität, New York). Advances 
in applied Mechanics Bd. II. X+ 
324 S. m. zahlr. Figuren. New York 1953. Academics 
Press Inc. Preis geb. 9.— $. 


Ph. Furtwängler, Allgemeine Theorie 
deralgebraischen Zahlen. Überarbeitet 
von H. Haase — W. Jehne, Hamburg. (Enzyklo- 
pädie der Mathematischen Wissenschaften mit Ein- 
schlußihrer Anwendungen Bd I 2, Heft 3, Teil II,19.) 
50 8. Leipzig 1953. B.G. Teubner Verlagsgesell- 
schaft. Preis geb. 6,— DM. 


Dr.-Ing. A. Nasvytis. Die Gesetzmäßig- 
keiten kombinatorischer Technik. 
(Wissenschaftliche Normung. Schriftenreihe, heraus- 
gegeben in Verbindung mit dem Seminar für Tech- 
nische Normung an der Technischen Hochschule Han- 
nover. Heft3.) VIII +103S. m. 29 Abb. Berlin/ 
Göttingen/Heidelberg 1953. Springer-Verlag. Preis 
geh. 9.— DM. 


Dr. Bernhard Sticker (Prof. a. d. Univ. Bonn). 
Fünfstellige Tafel der Trigonome- 
trischen Funktionen. Ausg. A für Altgrad- 
teilung. 52 S. m. 1 Interpolationstafel. Bonn 1954. 
Ferd. Dümmlers Verlag. Preis geb. 7,60 DM. 


Dr.-Ing. W. Grossmann (Prof. a.d. T. H. Hannover). 
Grundzüge der Ausgleichungsrech- 
nung nach der Methode der kleinsten Quadrate 
nebst Anwendungen in der Geodäsie. VIII + 261 S. 
Berlin/Göttingen/Heidelberg 1953. Springer-Verlag. 
Preis geb. 19,80 DM. 


Uno Olsson, Non-Circular Cylindrical 
Gears. XI-+216S. Stockholm 1953. Acta Poly- 
technica. Mechanical Engineering Series Vol. 2 Nr. 10. 
Preis 18,— Sw. Kr. 


Colloque sur les fonctions de plusieur 
variables. Tenu & Bruxelles du 11. au 14. mars 
1953. 159 S. Liege-Paris 1953. Verlag Georges Thone, 
Masson et Cie. Preis geh. 250 francs belges. 1800 
francs frangais. 


H. Prüfer, Projektive Geometris. Aus 
dem Nachlaß herausgegeben von G. Fleddermann und 
G. Köthe. 2. unveränd. Aufl. VII+314S. m. 
254 Abb. Leipzig 1953. Akademische Verlagsgesell- 
schaft Geest & Portig K.G. Preis geb. 9,— DM. 


Dr.-Ing. R. Zurmühl, Praktische Mathe- 
matik für Ingenieure und Physiker. 
XI +481S. m. 114 Abb. Berlin/Göttingen/Heidel- 
berg. Springer-Verlag 1953. Preis geb. 28,50 DM. 


Dr. phil. habil. Rudolf Beyer, Kinematische 
Getriebesyntbese. Grundlagen einer quanti- 
tativen Getriebelehre ebener Getriebe. Für den Kon- 
strukteur, für die Vorlesung und das Selbststudium. 
VII + 217 S. m. 258 Abb. Berlin/Göttingen/Heidel- 
berg 1953. Springer-Verlag. Preis geb. 36,— DM. 


Dr. Fritz Stüssi (o. Prof. an der E. T. H. Zürich), 
Vorlesungen über Baustatik. BandI: 
Statisch bestimmte Systeme — Spannungsberech- 
nung — Elastische Formänderungen — Stabilitäts- 
probleme — Seile. 2. Aufl. 370S. m. 336 Abb. 
Basel/Stuttgart 1953. Verlag Birkhäuser. Preis geb. 
40,05 Sfrs, brosch. 35,90 Sfrs. 


L. Howarth, H.B. Squire and C.N.H. Lock, Mo- 
dern Developments in Fluid Dyna- 
mies. High Speed Flow. (The Oxford 
Engineering. Science Series.) 2 Bände. XVI+ VIII + 
8758. m. 288 Abb. u. 26 Tafeln. Oxford 1953. Cla- 
rendon Press. Preis geb. 84.— s. 


Dr. H. Blasius (Oberstudienrat i. R. der Ingenieur- 
schule Hamburg), Höhere Mathematik. 
Mathematische Grundlagen vom technischen Stand- 
punkt. 294 S. m. 168 Fig. und 78 Aufgaben. Ham- 
burg 1954. Boysen & Maasch Verlag. Preis kart. 
8,— DM, geb. 10,— DM. 


Dr. E.Trost (Prof. am Technikum Winterthur). 
Primzahlen. (Elemente der Mathematik vom 
höheren Standpunkt aus, Bd. II.) 95 S. Basel/Stutt- 
gart 1953. Verlag Birkhäuser. Preis brosch. 13,50 SFr. 


A. Blanc-Lapierre (Prof. & la Facult6 des Sciences 
d’Alger) und R. Fortet (Prof. & la Facult6 des Sciences 
de Paris, Th6orie des fontions alöa- 
toires. Applications a divers phe- 
nomenes de fluetuaticen. (Collection 
d’ouvrages de mathematique & l’usage des physiciens 
publi6e sous la direction de G. Darmois),. XVI-+ 
693 8. m. zahlr. Abb. Paris 1953. Masson & Cie. 
Preis brosch. 6000 frs., geb. 6500 frs. 


N. I. Muskhelishvili, Some basic problems 
ofthe mathematicaltheory ofelasti- 
eity. 3. verb. und vermehrte Auflage. (Aus dem 
Russischen übersetzt von J. R. M. Radok.) XXXI + 
704 S. m. 66 Abb. Grooningen (Holl.) 1953. Verlag 
P. Noordhoff Ltd. Preis geb. 38,— fl. 


N. I. Muskhelishvili (Math. Institut Tiflis), Sing u- 
larintegralequations. 2. Aufl. (Aus dem 
Russischen übersetzt von J. R. M. Radok.) 4478. 
m. 16 Abb. Grooningen (Holl.) 1953. Verlag P. 
Noordhoff. N. V. Preis geb. 28,50 fl. 
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Dr.-Ing. Rudolf Richter (Em. Prof. a. d. T.H. 
Karlsruhe) mit Beiträgen von Dr.-Ing. Robert Brüder- 
link (o. Prof. a. d. T.H. Aachen). Elektrische 
Maschinen. Band2: Synchronmaschi- 
nen und Einankerumformer. Zweite, 
verb. Aufl. XXVIII + 707 S. m. 519 Abb. Basel/ 
Stuttgart 1953. Verlag Birkhäuser. Preis geb. 
46.80 fr. 

Winfried Oppelt, Kleines Handbuch Tech- 
nischer Regelvorgänge. 4488. m. 353 Abb. 
u. 110 Tafeln. Weinheim/Bergstr. 1954. Verlag 
Chemie GmbH. Preis geb. 31,80 DM. 


Review of InputandOutputEquipment 
Used in Computing Systems (Joint AIEE- 
IRE-ACM Computer Conference). 1428. m. zahlr. 
Abb. New York 1953. American Institute of Elec- 
trical Engineers. Preis geh. 4,— $. 


Dr.H.von Sanden (o. Prof. an der T. H. Hannover). 
Praktische Mathematik (Teubners Mathe- 
matische Leitfäden, Bd. 44), 3. erw. und umgearbeitete 
Aufl. 128 S. m. 25 Abb. Leipzig 1953. Verlag_B. G. 
Teubner. Preis kart. 3.20 DM. ut 


NACHRICHTEN 


Tagung der Deutschen Mathematiker-Vereinigung 
in Mainz vom 20.—26. September 1953 


Die Tagung war mit rund 200 Teilnehmern, unter 
denen sich wiein den Vorjahren eine Reihe namhafter 
ausländischer Gäste befand, sebr gut besucht. Das 
reichhaltige Programm wies gegen 100 Vorträge und 
Berichte auf, darunter die folgenden mit vierzig bis 
sechzig Minuten Dauer: 

Kam ke-Tübingen, 
Rechenmaschinen. 

Grün-Würzburg, Dimensionsgruppen, 
Ring und Hausdorffsche Formel. 

Kneser-Tübingen, Wertfunktion und Versiche- 
rung. Ein Beitrag zu den Grundlagen der 
Ökonometrie. 

Schäfke-Mainz, Entwicklung analytischer Funk- 
tionen nach Biorthogonalsystemen spezieller 
Funktionen und Anwendungen. 

Brelot-Grenoble, Neuere Formulierungen für 
das Dirichletsche Problem und Prinzip. 

Reiter-Chieago, Methoden und Ergebnisse der 
modernen harmonischen Analyse. 

Den y-Strasbourg, La recherche des noyaux pour 
lesquels il existe un theor&me du balayage. 

Hamburger-Köln, Zur Tbeorie des Residual- 
spektrums linearer Transformationen im Hilbert- 
schen Raum. 

Reidemeister-Marburg, Rauminhalt von 
Polyedern. 

Hornich-Graz, Lineare partielle Differential- 
gleichungen, deren Koeffizienten Polynome sind. 

Hlawka-Wien, Bericht über Geometrie der 
Zahlen. 

Hasse-Hamburg, Die Zetafunktionen arithme- 

tischer Funktionenkörper I. 

Deuring-Göttingen, die Zetafunktionen arith- 
metischer Funktionenkörper II. 

van der Waerden-Zürich, Jacobische Man- 
nigfaltigkeiten. 

Die übrigen Vorträge behandelten Fragen aus fol- 
genden Gebieten: Algebra, Gruppentheorie, Topo- 
logische Algebra, Zahlentheorie, Mengenlehre, Logik, 
Limitierungstheorie, Komplexe Funktionentheorie, 
Partielle Differentialgleichungen mit Rand- und 
Eigenwertproblemen, Funktionalanalysis, Geometrie 
einschließlich Grundlagen der Geometrie, Algebraische 
Geometrie, Topologie, Differentialgeometrie, Prak- 
tische Analysis, sowie zahlreiche andere Einzelfragen. 

Anschließend an die Mitgliederversammlung be- 
richtete Herr Behn ke-Münster über die Arbeit des 
deutschen Unterausschusses der IMUK. Den Ab- 
schluß der Tagung bildete eine Führung durch das 
Institut für Praktische Mathematik der Technischen 
Hochschule Darmstadt. F.Lösch. 

Stuttgart. 


Prinzipien elektronischer 


Liescher 


Am 4. September 1953 verstarb im Alter von 
70 Jahren Präsident a. D. Prof. Dr. Karl Willy 
Wagner an einem Herzschlag. 


R.v. Mises fr 
Am 20. Juli 1953 verstarb in Boston (Mass.) im 
Alter von 70 Jahren der ordentliche Professor der 
Aerodynamik und der angewandten Mathematik an 
der Harvard-Universität 


Dr.Dr.h.c. Richard von Mises. 


Es gibt wohl kaum ein Gebiet der angewandten 
Mathematik, das dieser vielseitige Wissenschaftler 
nicht durch wertvolle Arbeiten bereichert hat. Im 
Jahre 1920 — damals Professor der angewandten 
Mathematik an der Universität Berlin — gründete 
von Mises die ‚Zeitschrift für angewandte 
Mathematik und Mechanik“, die er bis zu seiner Be- 
rufung nach Istanbul 1934 leitete. Die Ziele, die von 
Mises der Zeitschrift in einem einführenden Auf- 
satz des ersten Heftes weitschauend gesteckt hat, 
verfolgt sie bis heute. Zu seinem 50. und seinem 
70. Geburtstage widmete sie ihm Sonderhefte. Schrift- 
leitung und Verlag werdenRichard von Mises 
ein ehrendes Andenken bewahren. 


Schriftleitung: Willers Verlag. 


München: Am 29. September 1953 verschied 
der Dozent an der Technischen Hochschule München 
Dr. Fritz Wenzl im 40. Lebensjahr. 


Berlin-Lichtenberg: Am 6. Januar 1954 ver- 
schied unerwartet Dr.-Ing. KarlBorkmannim 
Alter von 49 Jahren. 


Internationaler Mathematiker-Kongress 


Der Internationale Mathematiker-Kongress 1954 
findet in der Zeit von Donnerstag, dem 2. September 
bis einschl. Donnerstag, dem 9. September in Amster- 
dam statt. Die Anschrift des Sekretariats dieses Kon- 
gresses ist: 


2e Boerhaavestraat 49, Amsterdam, The Netherlands. 
Der Tagungsbeitrag beträgt 50 Gld. und ist möglichst 
bald an die ‚„Amsterdamsche Bank N.V., Filiaal 
Leiden, The Netherlands‘ einzusenden. Es finden in 
7 Sektionen 


(Sektion I Algebra und Zahlentheorie, 
II Analysis, 
III Geometrie und Topologie, 
IV Wahrscheinlichkeit und Statistik, 
V Mathematische Physik und angew. 
Mathematik, 
VI Logik und Grundlagen-Theorie, 
VII Philosophie, Geschichte und Unterricht) 


45 ganz- und halbstündige Vorträge statt. Außerdem 
können kürzere Vorträge beim Sekretariat des Kon- 
gresses angemeldet werden. Dieser Anmeldung sollte 
eine kurze Inhaltsangabe von nicht mehr als 400 Wor- 
ten beigefügt werden. 
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